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Abstract 



' The outstanding problem in cosmology today is undoubtedly the origin and evolution 

, of large scale structure. In this context, no model has proved as successful as the 

' standard Cold Dark Matter (CDM) model, based on a flat (scale-invariant) spectrum of 

, density fluctuations growing under gravitational instability. Nevertheless, pressure has 

recently been exerted on the model both from analyses of large scale clustering in the 



galaxy distribution and from the measurement of the level of microwave background 
anisotropics by the Cosmic Background Explorer (COBE) satellite. 
O , In this report we aim to give a unified view of the Cold Dark Matter model, beginning 

' with the creation of perturbations during an inflationary epoch and pursuing it right 

\ through to comparison with a host of observations. We discuss in detail the evolution of 

density perturbations in Fricdmann universes, utilising the fluid flow approach pioneered 
by Hawking, and provide a simple derivation of the Sachs-Wolfe effect giving large 
' angle microwave background anisotropies. We illustrate the means by which inflation 

^ \ provides an initial spectrum of inhomogeneities, the spectrum having a shape which 

can be readily calculated in a given inflationary model. We also include a discussion 
of the generation of long wavelength gravitational waves, which have recently been 
recognised as having a potentially important role with regard to microwave background 
anisotropies. 

Although the standard CDM model is based on a scale-invariant spectrum, the 
generic prediction of simple inflationary models is for a power-law spectrum, tilted 
away from scale-invariance to provide extra large-scale power. For many models such 
as chaotic inflation, this tilting is rather modest. However, in several models, such as 
power-law inflation, extended inflation and natural inflation, the tilting of the spectrum 
can be more dramatic, and potentially useful in the light of observations indicating a 
deficit of large-scale power in the galaxy distribution. In the former two of these, there 
is the interesting extra of a substantial production of gravitational waves. 



^To appear, Physics Reports. Original version entitled 'The Spectral Slope in the Cold Dark Matter 
Cosmogony'. 



We then discuss observational constraints on cold dark matter cosmogonies based on 
power-law spectra. We examine a range of phenomena, including large angle microwave 
background fluctuations, clustering in the galaxy distribution, peculiar velocity flows, 
the formation of high redshift quasars and the epoch of structure formation. One of 
our aims is to compute the current constraints on both the shape of the spectrum as 
defined by the spectral index n, and its normalisation as defined by the usual quantity 
as = 1/^8. 

We end by discussing briefly some variants on the CDM model, such as the incorpo- 
ration of a hot dark matter component, or the introduction of a cosmological constant 
term. We do not however investigate them in the depth that we do the tilted CDM 
models, though the techniques illustrated throughout the paper provide the background 
required for such an undertaking. 

E-mail addresses: arl @ uk.ac.sussex.starlink ; lyth @ uk.ac.lancs.ph.vl 
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1 Introduction 



For many (and probably most) cosmologists, the outstanding problem at present is to 
develop an understanding of the origin and evolution of large scale structure in the universe, 
with the ultimate goal of explaining such phenomena as the epoch of galaxy formation, the 
clustering in the galaxy distribution and the amplitude and form of anisotropies in the 
microwave background. The present time is a particularly exciting one, as experiments 
through the last decade have presented great leaps in our knowledge of the universe around 
us, and experiments proposed for the years to come promise to further revolutionise our 
view. Attempts to model the formation of large scale structure already face an impressive 
array of constraints they must obey. 

Gravity being the only universally attractive long range force, the formation of structure 
must begin with a small perturbation in the otherwise uniform distribution of matter in 
the early universe. Two alternative origins have been proposed for this perturbation. The 
first is that it originated as a quantum fluctuation during inflation which became 'frozcn-in' 
after horizon exit; such a mechanism typically yields a Gaussian, adiabatic perturbation 
with a more or less scale-invariant spectrum, of the kind which has long been favoured on 
grounds of simplicity (Harrison 1970; Zcl'dovich 1970). The second is that it was caused by 
topological defects like cosmic strings (Kibble 1976; Zel'dovich 1980; Vilenkin 1981) leading 
to a more complicated type of perturbation. 

Partly because of its simplicity, the first possibility has been investigated in far more 
detail than the second, and it forms the subject of the present report. A crucial qTicstion for 
this hypothesis is the nature and amount of the dark matter in the universe. Soon after the 
need for such matter came to be widely accepted in the early 1980's, it became clear that 
the hypothesis fails completely if the dark matter consists of massive neutrinos, because 
their thermal motion wipes out small scale structure (White, Prenk & Davis 1983). Given 
the failure of this hot dark matter (HDM) model, attention turned to the other extreme, 
of matter which has by definition negligible random motion. In its standard form this cold 
dark matter (CDM) model assumes that the universe has fiat spatial geometry, critical 
matter density and a spectrum which is precisely scale invariant. Then it is defined by the 
present Hubble parameter Hq, the baryonic matter density and the normalisation of 
the spectrum. For the Hubble parameter, a value Hq = 50kms~^Mpc~^ is usually taken, 
both in order to have an acceptably old universe and to make the model itself reasonably 
viable. The baryonic content is taken as accurately determined by the standard big bang 
nucleosynthesis picture. 

For the first few years of its existence, this standard form of the CDM model was very 
successful in describing the observed features of structure in the universe, from the scale 
of galaxies upwards (Peebles 1982b; Blumenthal et al 1984; Bardeen et al 1986; White 
et al 1987; Frenk et al 1988; Efstathiou 1990). More recently, though, the model has run 
into difficulties and is widely felt to be in need of modification. One possibility is to 'tilt' 
the spectrum away from its scale invariant form. Another is to allow the cold dark matter 
density to be less than the critical value, making up the difference with either a cosmological 
constant or hot dark matter. Of course there are other possibilities, but these are the ones 
that have received the most attention at the present time. 
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1.1 Generalities 



The standard hot big bang model (see eg Kolb & Turner 1990), which provides the frame- 
work in which gravitational instability scenarios are set, consists of a collection of parameters 
which are more or less well known. The rate of expansion of the universe is described by 
the Hubble parameter H = a/a {a being the scale factor) , whose present value is commonly 
parametrised as 

i?o = 100/ikms^^Mpc"^ , /i G [0.4, 1] (1.1) 

With c = 1, this can be written accurately as -f^o = /i/3000 Mpc~^. The energy density is 
usually written in terms of the critical density 

p -^juMl (12) 

as Q = p/ pc- It is convenient to regard any cosmological constant as a time independent 
contribution to p, and then the case $7 = 1 corresponds to flat spatial sections. 

Different direct measurements of the Hubble parameter give h in the range .4 to 1. The 
present value of is even more uncertain than h, with any value between roughly 0.1 and 1 
permitted by various experiments. The theorist's favourite, for fairly well founded reasons 
we shall examine later, is the = 1 case which describes a universe with flat spatial sections. 
Perhaps surprisingly, the density in ordinary matter is far better known, provided one takes 
the reasonable step of believing the phenomenally successful nucleosynthesis calculations 
(see eg Walker et al 1991), as lying in the rangeQ 

0.010 < nsh'^ < 0.016 (1.3) 

It is immediately clear that for the theorist's 0, = 1 option to be correct, the bulk of the 
universe must be in some as-yet-undetected form of dark matter. The nature of this dark 
matter turns out to be a crucial ingredient in models of structure formation. 

The age of the universe has also been estimated, most reliably from the ages of globular 
clusters as upwards of 13 Gyr. In a flat (0 = 1) universe, the age is simply t = 2i?o"V3 = 
6.5//1 Gyr, and so one requires /i < .5 for consistency. Combined with the result h > A 
from direct observation, a value /i ~ .5 is therefore preferred over higher values. Should 
measurements giving a high Hubble parameter prove correct, one would be forced either to 
let be significantly less than 1 or to introduce unusual dynamics such as a cosmological 
constant in order to reconcile them. 

In itself, the hot big bang model has nothing to say on structure formation, remaining 
homogeneous for all time. However, prompted by certain conceptual problems with the hot 
big bang — why should the initial conditions be so homogeneous, why should the universe be 
so close to the critical density, what prevented an overabundance of magnetic monopoles, etc 
— cosmologists devised the paradigm of inflation, whereby the universe in its earliest stages 
undergoes a period of accelerated expansion a > (Guth 1981; Kolb &: Turner 1990; Linde 
1990). It was later realised that inflation, by a mechanism we shall shortly discuss, could 
generate small density irregularities on large (superhorizon) scales, and that these would 

^In the context of cosmology ordinary matter is usually referred to as baryonic matter since the baryons 
(protons and neutrons) far outweigh the electrons. Thus fis is referred to as the baryon density. 
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naturally lead to the gaussian, adiabatic, scale invariant density perturbation which had 
already been proposed by Harrison and Zel'dovich, and adopted on grounds of simplicity by 
most large scale structure workers. With inflation to provide an initial spectrum of density 
inhomogeneities, the spotlight is transferred to developing an understanding of whether or 
not such a spectrum can evolve so as to reproduce large scale structure observations. The 
evolution of the spectrum depends sensitively on the assumptions (or observations) one 
makes as to the matter content of the universe. 

In attempting to explain the increasing range of observations, no model has provided as 
dramatic an advance in understanding as that which has come to be known as the standard 
Cold Dark Matter (CDM) model (White et al 1987; Prenk et al 1988; Efstathiou 1990; 
Ostriker 1993). As implied by its name, the dark matter is assumed to be cold, which for 
most purposes can be taken to mean non-relativistic. By definition, dark matter does not 
interact significantly with more conventional forms of matter by any means other than grav- 
itationally, and in particular is beneficial for structure formation in that it is not subject 
to pressure forces from interaction with radiation which prevent baryonic density inhomo- 
geneities on scales smaller than superclusters from collapsing before radiation decouples 
from matter. Structure can thus start to form earlier in the dark matter providing initial 
gravitational wells to kick-start structure formation in baryonic matter after decoupling. 

Over the last decade, the cold dark matter model has been extensively explored, most 
notably by Davis, Efstathiou, Frcnk and White in a series of papers. The assumptions 
above provide a scenario in which calculations and predictions are readily framed. As 
long as they remain small, the inhomogeneities can be split up into each comoving scale 
evolving separately, enabling analytic progress to be made. On scales much larger than 8h~^ 
Mpc, such calculations arc still applicable today. On smaller scales the perturbations reach 
the nonlinear regime and become mode-coupled, and more complex methods, primarily 
large iV-body computer simulations, must be employed to produce predictions. In the 
late eighties, the CDM model could boast impressive successes in reproduction of galaxy 
clustering statistics, structure formation epochs and peculiar velocity flows, while keeping 
fluctuations in the microwave background at a level safely below the observational upper 
limits. The only addition to the model which observation required is that which remains 
the most controversial to the present day — the question of Massing. 

The concept of biassing arose from the realisation that in the theory one calculates the 
distribution of matter in the universe, but in observations one measures the distribution 
of galaxies. In general these need not be the same; the question is whether or not 'light 
traces mass'. The unbiased CDM models appeared unable to produce enough clustering 
in the galaxies, and the standard assumption became that light does not trace mass. The 
clustering of galaxies, as measured by the two-point correlation function, was assumed to be 
greater than that of the mass by a multiplicative constant, the bias parameter b, originally 
thought to lie anywhere in the range 1 to perhaps 2.5. With the incorporation of the extra 
free parameter, CDM was able to reproduce galaxy distributions, at least out to Wh^^ Mpc. 

More recently, standard CDM model has come under fire from two directions, and 
though neither is as yet completely conclusive these are certainly troubled times for CDM. 
The first problem is that on scales 10 to 100 Mpc the standard CDM model does not give 
the correct scale-dependence for the galaxy correlation function; it seems to have too little 
power on the large scales in this range, relative to the power on smaller scales. Perhaps 
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the most prominent data leading to this conclusion is the APM survey (Maddox et al 1990, 
1991) giving the galaxy angular correlation function, but also important are the QDOT 
survey (Saunders et al 1991), observations of X-ray galaxy clusters (Lahav et al 1989) and 
of radio galaxies (Peacock 1991; Peacock & Nicholson 1991), and optical galaxy redshift 
surveys such as CfA (Vogeley et al 1992) and the Southern Sky Redshift Survey (Park, 
Gott & da Costa 1992). 

The second problem for standard CDM is that there seems to be a conflict between 
the normalisation of the spectrum of the perturbation, as specified for instance by the 
usual quantity as = ^/bg, which is required by different types of observation. Perhaps 
the most reliable normalisation comes from the first positive measurement of microwave 
background anisotropics, reported by the COBE DMR collaboration in April 1992 (Smoot 
et al 1992). It probes the density perturbation on scales of order 10^ Mpc, and requires 
cTg = .95±.2. Second, there are surveys of the galaxy distribution and bulk flow which probe 
the normalisation on scales of order 10 to 50Mpc, notably from the QDOT (Kaiser et al 
1991) and POTENT (Dekel et al 1992) groups; these seem to require erg in the range .7 to 1.1, 
which is compatible with the COBE normalisation. The problem comes with observations 
of the pairwise velocity dispersion of galaxies on scales < 3 Mpc, which according to most 
calculations seem to require ag < .5. Gelb, Gradwohl and Prieman (1993) have claimed 
recently that through non-linear effects these velocity dispersion observations probe the 
normalisation of the density perturbation on a wide range of scales, extending into the 
regime 10 to 50 Mpc probed by the bulk flow data. If that is so, the problem is of a 
different type from that represented by the galaxy correlation function data, and cannot 
cannot really be paraphrased by stating that there is a relative lack of power on large scales. 
Recent calculations of the predicted galaxy cluster abundance (White, Efstathiou &: Frenk 
1993) also indicate less small-scale power than indicated by the bulk flow data. 

Faced with troubling observations but a generally successful model, the impulse of the 
community is to re-examine, and if necessary relax, the underlying assumptions of the 
standard CDM model. Several means are available to achieve this, which we shall discuss 
later. One of the purposes of this paper is to examine the option which, while not necessarily 
the most useful, is to the particle cosmologist the most compelling. That is to examine the 
form of the primeval spectrum. 

Standard CDM is based on a scale-invariant spectrum. Whilst initially this was based 
entirely on aesthetics (it is for instance the only option which does not require a choice of 
normalisation scale), it was later justified on the basis of being a prediction of inflationary 
cosmology. Nevertheless, it has long been known, at least within the inflation community, 
that the predicted spectrum is only approximately flat, and that within any given model 
there are readily calculable deviations from scale invariance. An example is provided by 
the standard 'chaotic' inflation scenario, (Linde 1983, 1987, 1990), wherein a logarithmic 
correction (see eg Salopek, Bond &; Bardeen 1989; Mukhanov, Feldman &: Brandcnbcrgcr 
1992; Schaefer & Shafi 1992) adds some power at large scales, giving a density contrast at 
horizon crossing which is 5-10% greater at 1000/i~^ Mpc than at Wh~^ Mpc. Hitherto, when 
our knowledge of the primeval spectrum was restricted to studies of local clustering up to 
scales of perhaps WOh^^ Mpc, these corrections where rightly regarded as insignificant. The 
recent results (Smoot et al 1992) from the Cosmic Background Explorer satellite (COBE) 
have changed this picture, by providing for the first time good estimates of the amplitude of 
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the spectrum on scales of upwards of 1000/i~ Mpc. The time has come for these correction 
terms to be taken seriously, and possibly the real test of inflation lies not in how flat the 
primeval spectrum is, but in the type and size of deviations from it. 

Some inflationary models predict an even more dramatic change from the Harrison- 
Zel'dovich case. The simplest is power-law inflation (Abbott &: Wise 1984b; Lucchin & 
Matarrese 1985), which can be realised via a scalar field evolving in a potential of exponential 
form and which leads to an expansion of the universe which is not the conventional nearly 
exponential growth, but rather has the scale factor growing as a rapid power-law in time. 
A second option is provided by the so-called extended inflation scenarios (La &; Steinhardt 
1989; Kolb, Salopek k. Turner 1990; Kolb 1991), wherein modifications to Einsteinian 
gravity allow inflation to proceed via a first-order phase transition. Yet another option 
is provided by 'natural inflation' (Freese, Frieman & Olinto 1990; Adams et al 1993), a 
specific case of the more general situation of a scalar field evolving near the top of an 
inverted harmonic oscillator potential which also gives a power-law spectrum. And finally, 
it is worth emphasising that although technically logarithmically corrected, models such 
as chaotic inflation produce spectra which are excellently approximated across large scale 
structure scales by a power-law, though typically much closer to the scale-invariant case 
than the genuine power-law models allow. Only in exceptional cases (we discuss a two- 
scale inflation model) does inflation seem capable of providing spectra tilted in the opposite 
direction, giving extra small-scale power. 

Another potentially vital aspect of inflationary models is that they can generate long 
wavelength gravitational waves (Fabbri k, Pollock 1983; Abbott k, Wise 1985; Starobinsky 
1985). Recently, a host of papers (Krauss & White 1992; Salopek 1992a, 1992b; Davis et 
al 1992b; Liddle & Lyth 1992; ; Lidsey & Coles 1992; Lucchin, Matarrese k Mollerach 
1992; Souradeep k Sahni 1992) have served to remind the community that while this 
can be a very small effect in some inflationary models, it is not necessarily the case that 
this be so. We discuss the generation of these modes during inflation, and emphasise the 
inflationary models in which this effect is most important. Typically, a strong gravitational 
wave component is associated with tilt, though the converse — that tilt necessarily implies 
gravitational waves — need not be true. 

1.2 About this paper 

The primary aim of this report is to provide a unified view of the cold dark matter model 
based on inflationary cosmology, examining the assumptions made and testing the obser- 
vational viability. In the first half of the paper, we concentrate on developing the theory 
of the origin and evolution of perturbations in Friedmann universes. To do this, we utilise 
an approach, initiated by Hawking (1966), which is somewhat different from the usual met- 
ric perturbation approach of Lifshitz (1946). In the perfect fluid approximation it leads 
to a straightforward derivation of the evolution of the density perturbation and related 
quantities, using fluid flow equations which closely resemble the non-relativistic ones. They 
cover a variety of physical situations, ranging from pressure-free adiabatic perturbations to 
isocurvature perturbations in two component fluids. On large scales they are applicable in 
a very general setting, making no assumptions about the type of dark matter. We include 
detailed discussions of peculiar velocities, which we use to provide a simple derivation of 
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the Sachs- Wolfe effect giving the large scale cosmic microwave background radiation (cmb) 
anisotropy arising from the density perturbation. We also discuss the anisotropy which may 
arise from gravitational waves. Then we go on to discuss the origin of these perturbations 
in inflationary cosmology, emphasising that their spectra typically exhibit deviations from 
scale invariance. We calculate the extent of the spectral tilt for the most popular models, 
and also calculate the relative contribution of gravitational waves to the cmb anisotropy. 

In the second half of the paper we move on to the observational side. We concentrate on 
CDM models, generalising standard CDM by incorporating arbitrary power-law spectra as 
motivated by inflation. Parametrising the inflation-generated spectrum of density pertur- 
bations by its normalisation erg and by its degree of tilt as specified by the spectral index n, 
we first map out forbidden regions in the n-crg plane coming from different types of obser- 
vation. We examine a range of observations spanning the range of scales from the present 
horizon size (microwave background anisotropies) down to the scale of typical galaxies, and 
find that if the observations are all taken at face value no choice of the parameter pair n-crg 
can fit them all. We end by briefiy discussing the possibility that the cold dark matter 
makes up only part of the critical density, with the remainder either a cosmological con- 
stant or hot dark matter. The emphasis in this second half is on providing an illustration of 
the techniques used to obtain the constraints, and we caution the reader that the detailed 
results may in places be rapidly superceded by improved observations. 

1.3 Some basic facts and assumptions 

Before commencing, let us note some important cosmological parameters and scales we shall 
need. 

For ease of comparison with other work we often take h'^ Mpc for the distance unit. 
The units are taken to be such that c = 1 and (in Section |5|) ?i = 1, and in these units the 
Planck scale is defined by mpi = = 1.2 x 10^^ GeV. 

The total energy density is taken to be critical, $7 = 1, and except in Section |9| we 
discount the possibility of a cosmological contribution to it. The present Hubble parameter 
is taken to be given hy h = 0.5. The most important scales are the following (Hogan, 
Kaiser & Rees 1982). Expressions are quoted first with the relevant factors of h to show 
the scaling, and then for our choice of h = 0.5. The mean density at the present epoch is 
po = 3H^ /8ttG, which in convenient astrophysical units is 

po = 2.78h'^ X lO^^M0(/i-^Mpc)-3 (1.4) 
= 6.94 X lO^°M0Mpc"^ for h = 0.5 (1.5) 

where Mq is the solar mass. We normalise the scalar factor to be a = 1 at the present, 
so that it is related to redshift simply by o = (1 + z)~^. With this choice, comoving 
and physical scales coincide at the present, with physical scales growing proportional to 
a. Thus, comoving scales are specified by their physical size at the present. The present 
horizon size|^ is simply 2Hq^ = 6000/i~^Mpc = 12000 Mpc. The last scattering surface is 
centred at a redshift z = 1080, taking on roughly the form of a gaussian of width Az = 80 

^This is somewhat notional, being the distance light could travel were the standard hot big bang extrap- 
olated back to its origin. With inflation, the true horizon size can be vastly greater, while in practice the 
last scattering surface provides a limit to the distance we can see which is close to this value. 
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(Jones & Wyse 1985). In comoving units, this width is around 7h~^Mpc = 14 Mpc, and for 
microwave features of this size or smaller the width must be taken into account. Matter- 
radiation equality occurred at a redshift z = 24000/i^. 

We shall be defining horizon crossing at a given epoch as A; = aH. We can thus give the 
horizon-crossing wavenumbers for these different important epochs. They are (with h = .5) 

^hor = eOOOMpc (1.6) 
^dec = ISOMpc (1.7) 
k-^ = 79Mpc (1.8) 



2 The density perturbation 

2.1 Cosmological perturbation theory 

Departures from homogeneity and isotropy are small in the early universe, and even at 
the present epoch on sufficiently large scales. In this situation one can use cosmological 
perturbation theory, which develops linear equations for the perturbations which are valid 
to first order. This means that the equations become exact if each perturbation is multiplied 
by a common parameter which tends to zero. A perturbation can be either a small change 
in some quantity which is non-zero in the limit of homogeneity and isotropy, or the entire 
value of a quantity which vanishes in that limit. 

Two essentially different ways of handling the perturbations have been developed in the 
literature. The usual one, due to Lifshitz (1946), works with coordinates so that pertur- 
bations of the metric tensor components are involved. The other, due to Hawking (1966), 
is a coordinate-free approach. The former approach is described in many textbooks and 
reviews (Weinberg 1972; Peebles 1980; Mukhanov, Feldman Sz Brandenberger 1992) and 
wc shall use it now to explain how cosmological perturbations may be divided into three 
independent modes. Then we shall abandon it in favour of Hawking's approach, which is 
simpler for the mode which involves the density perturbation. 



The metric perturbation approach 

In discussing perturbations, one can employ any coordinate system which reduces to one of 
the standard choices in the limit where the perturbations vanish. In this limit the universe 
is homogeneous and isotropic, and its energy momentum tensor necessarily has the perfect 
fluid form 

T"" = pgi^" + {p + p)u^'vy (2.1) 

where p is the energy density, p is the pressure and is the fluid four-velocity. The usual 
choice of coordinates in this limit corresponds to the line element 

ds^ = - a'^{t)5ijdx'dx^ (2.2) 

Then p and p depend only on t, and = (1, 0, 0, 0). 

The perturbations in the metric tensor components may then be defined by 

ds^ = {l + hoo)dt^ + 2hoidtdx' - a^{t){5ij + hij)dx'dx^ (2.3) 
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Ignoring the phenomena of particle diffusion and free-streaming the energy-momentum ten- 
sor continues to have the perfect fluid form, its perturbation being defined by the pertur- 
bations 5p, 6p and Su^ = (0,t;*), where v'' is the three- velocity of the fluid in the chosen 
coordinate system. 

Einstein's field equation yields a set of linear partial differential equations involving the 
metric perturbations, 5p, 6p and v^. They were written down first by Lifshitz (1946). For 
the case of critical density which we are considering here, each perturbation 'lives' in flat 
space, and can be expanded as a Fourier series in a comoving box. One then has ordinary 
differential equations which do not couple separate Fourier coefficients. The following three 
modes propagate independently 

• The perturbations 5p and 5p, and the irrotational part of u*. 

• The rotational part of (vorticity) . 

• Gravitational waves, characterised by the traceless transverse part of hij. 

Each of the first two modes is also associated with a definite type of metric perturbation, 
which we have not bothered to specify. These modes are usually referred to as respectively 
scalar, vector and tensor modes because of the spatial transformation properties of the 
metric components, but more instructive labels are 'density', 'vorticity' and 'gravitational 
wave' modes. 

Vorticity decays with time, making the early universe less homogeneous than the present 
one, so it is presumably absent. Gravitational waves on cosmological scales might be gen- 
erated during inflation and we shall need to consider them later. For the moment, though, 
our interest is in the mode involving the density perturbation. This mode describes the 
perturbation in the density and motion of the matter and radiation in the universe, and can 
be generalised to include the effects of particle diffusion and free streaming (in principle it 
then couples to the other modes, but in practice the coupling is completely negligible). 

For a given perturbed space-time the perturbations in the metric and the energy- 
momentum tensor are uniquely defined only when a definite coordinate choice has been 
made. A coordinate choice, or more loosely a family of coordinate choices, is called a gauge. 
A coordinate choice can be thought of as a slicing of spacetime into spacelike hypersurfaces 
(with constant time coordinate) plus a threading of spacetime into timelike lines (with con- 
stant space coordinates). For the density mode the choice of gauge has a big influence on the 
appearance of the equations, and largely for this reason the metric perturbation formalism 
of Lifshitz was worked over by many subsequent authors as described in several text books 
and reviews (Weinberg 1972; Peebles 1980; Mukhanov, Feldman &; Brandenberger 1992). 

The most widely used gauge, at least until recently, was the synchronous gauge. In full 
generality the synchronous gauge employs arbitrarily chosen geodesies and the hypersurfaces 
orthogonal to them. Thus in Eq. ( |2.3| ) hoi = and hoo = 1. In practice, solutions of the 
synchronous gauge equations are selected which correspond to geodesies which become 
comoving in the limit of very early times (Lyth & Stewart 1990b), and from now on the 
term 'synchronous gauge' will denote this particular choice. 

In an influential paper, Bardeen (1980) considered several other gauges, among them the 
'comoving' gauge in which comoving worldlines are used instead of geodesies, along with the 
hypersurfaces orthogonal to them which are conveniently termed 'comoving hypersurfaces'. 
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The comoving gauge was shown by Lyth (1985) to lead, on each scale, to particularly simple 
equations well before horizon entry as will be explained below. Bardeen also developed a 
'gauge invariant' formalism which we shall not need. 

The fluid flow approach 

A description of perturbations which makes no mention of the metric perturbation was 
proposed by Hawking (1966). For the density mode of the perturbation it is simpler than 
the metric perturbation formalism, at least when particle diffusion and free streaming are 
negligible. For that mode, it boils down to evolving the quantities of interest along each 
comoving worldline, these quantities being the energy density p, the pressure p and the 
locally defined Hubble parameter H, all measured with respect to the worldline. The 
evolution equations are a pair of fluid flow equations, which hardly differ in form from their 
non-relativistic counterparts. 

To implement Hawking's approach for this type of perturbation, one has to introduce 
spacelike hyper surf aces, on which each quantity is given as the sum of an average plus a 
perturbation. As in the metric perturbation approach, one can derive differential equations 
for the Fourier components of the perturbations. The form of the equations is uniquely 
determined once the hypersurfaces have been chosen, and is of course the same as that 
obtained from the metric perturbation approach, using any gauge whose time coordinate is 
constant on each hypersurface. The big advantage of this approach for the density mode 
is that it eliminates all mention of the perturbation in the choice of space coordinates 
which is defined by threading timelike lines through the spacelike hypersurfaces, and of the 
associated metric perturbations hij and hiQ. One needs only the perturbation /iqo, specifying 
the proper time separation dr of adjacent hypersurfaces, which is obviously defined once 
the slicing into hypersurfaces has been defined. Because of this, we shall from now on use 
the term 'gauge' to denote simply a slicing of spacetime into hypersurfaces. 

Hawking's original implementation of his approach was valid only for the case of zero 
pressure gradient. It was carried through for nonzero pressure gradient by Olson (1976) 
in what amounted to the synchronous gauge, and by Lyth and Mukherjee (1988) in the 
comoving gauge. The approach has since been discussed in detail, and extended to the case 
of non-critical density and non-isotropic stress, by Ellis, Bruni and others, in a series of 
papers which can be traced from that of Bruni, Dunsby and Ellis (1992). 

However it is treated, the list of perturbations given above is incomplete because it does 
not mention the constituents of the universe. For our purpose these are nuclei, electrons 
and cold dark matter particles (matter), as well as photons and neutrinos (radiation). Until 
Section |9| the neutrinos are assumed to be massless, corresponding to the absence of hot dark 
matter. The usual assumption is that long before horizon entry the perturbations satisfy 
the adiabatic condition. This implies that there is a common radiation density contrast and 
a common matter density contrast, related by Spm/pm = \^Pr/ Pr- A less natural initial 
condition is the isocurvature condition, that there is negligible perturbation in the total 
energy density. The most general initial density perturbation is the sum of an adiabatic and 
an isocurvature one. For an adiabatic perturbation, the cosmological perturbation theory 
described above gives an unambiguous description of very large scale perturbations, which 
enter the horizon well after matter domination. In particular it allows one to describe their 
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effect on the cosmic microwave background radiation (cmb), which is their only observable 
signature, the necessary formalism being worked out first by Sachs and Wolfe (1969). 

On smaller scales additional considerations become necessary. At the very least, one has 
to take into account the separate evolution of the perturbations of the radiation and the 
dark matter, on the assumption that both are perfect fluids. The necessary equations in the 
synchronous gauge are given in, for example, the text of Peebles (1980). Simpler equivalent 
equations in the comoving gauge were derived by Kodama and Sasaki (1984, 1987) using the 
metric perturbation approach, and rcdcrivcd by Lyth and Stewart (1990b) using the fluid 
flow approach. To obtain accurate results one usually has to go beyond the perfect fluid 
approximation, to deal with free streaming and diffusion of the radiation. The necessary 
formalism is given in the text of Peebles (1980) and the review of Efstathiou (1990), in the 
synchronous gauge. 

In this report we give a treatment of the perturbations which is complete as far as 
large scales are concerned and is also rather simple. Using fluid flow equations and the 
comoving gauge, we use Hawking's approach to track an adiabatic density perturbation 
from a supposed inflationary origin to the present. A simple treatment of large scales is 
given, including the Sachs-Wolfe effect. Then the same thing is done for 'isocurvature' 
perturbations, using a 'two-fluid' approach. Finally, a possible contribution to the large 
scale cmb anisotropy from gravitational waves is described, starting with the supposed 
inflationary origin of such waves and evolving them to the present. 

On smaller scales our treatment is less complete because the effects of free streaming 
and diffusion are not described (save in a very simple approximation). Nevertheless we 
are able to deduce the qualitative behaviour of the transfer function relating the final, 
evolved density perturbation to the primeval one. Parametrisations of the adiabatic transfer 
function which various authors have obtained by taking into account free-streaming and 
(sometimes) diffusion are considered at this point, with an eye to the second half of the 
paper where theory and observation arc compared. The CDM prediction for the smaller 
scale cmb anisotropy is also displayed, again without detailed discussion. 

Although the cold dark matter scenario is the focus of the report as a whole, the part of 
it dealing with cosmological perturbation theory has much wider application. In particular, 
the treatment of perturbations remains valid before horizon entry even if the dark matter 
is not cold. Provided that there is no cosmological constant it remains valid to the present, 
for those very large scales which dominate the large scale cmb anisotropy. The effect of a 
cosmological constant is easy to evaluate using the formalism, though we have not given 
the details. 

2.2 Fluid flow equations 

We are supposing that the early universe may be regarded as a fluid, with a smoothly 
varying four- velocity field u^. At each point in space-time, a comoving observer is defined 
as one with this four-velocity. Relative to a comoving observer, the momentum density is 
zero. More formally, at each point in space-time the time-space components Toj of energy- 
momentum tensor vanish, in an orthonormal basis such that u'^ = (1,0,0,0). We shall 
call such a basis a comoving one. At least at sufficiently early and late times one can 
use the perfect fluid approximation, which means that in a comoving basis the space-space 
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components have the form Tij = —p5ij, where p is the pressure. The energy density Tqq in 
this basis is called p. 

Around a given point in space-time, it is convenient to use a locally inertial coordinate 
system in which the fluid is momentarily at rest. Then = (1, 0, 0, 0) at the point, and in 
a small region around it the space components define the three velocity of the fluid. The 
remaining component is given by 

u° = {1 - u'ui)^/^ (2.4) 
At the point, the Einstein field equation is equivalent to 

p = -3H{p + p) (2.5) 

n^ = (2.6) 
p + p 

Here, di (i=l,2,3) denotes differentiation with respect to the locally inertial space coordi- 
nates, and the dot denotes differentiation with respect to the locally inertial time coordinate, 
which may be identified with proper time along the comoving worldline passing through the 
point. The first equation is equivalent to the energy conservation equation dE = —pdV for 
the energy ii^ in a comoving volume. The second is the familiar 'acceleration equals force 
divided by mass', except that the mass density is replaced hy p + p (not just by the energy 
density p). 

The spatial derivatives of the three- velocity may be split uniquely into an antisymmetric 
vorticity LOij, a symmetric traceless shear crjj, and a locally defined Hubble parameter 

diUj = LOij + aij + 5ijH (2.7) 

From Eq. (|2.4| ), the derivatives of the time component vanish at the point where the fiuid 
is at rest, 

iio = diUQ = (2.8) 

These equations hold at each space-time point, in a locally inertial coordinate system in 
which the fluid is instantaneously at rest. 

If the vorticity vanishes throughout some region of space-time, there exist hypersurfaces 
orthogonal to the flow lines, which we shall call comoving hypersurfaces. They are the 
surfaces of simultaneity in the neighbourhood of each comoving observer. 

The unperturbed universe 

In a perfectly homogeneous and isotropic universe, the vorticity and shear are zero, and H, 
p and P are constant on each comoving hypersurface. Between any two comoving hyper- 
surfaces, the proper time interval along a comoving worldline is independent of position, 
and may be used as a time coordinate labelling the hypersurfaces. From the equation of 
motion Eq. ( |2.6| ) one may derive the Raychaudhuri equation (Raychaudhuri 1955, 1979; 
Ehlers 1961) 

H = -H^-^{p + 3p) (2.9) 
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Combining it with the energy conservation equation Eq. (^]^) leads to the Priedmann equa- 
tion 

H' = ^p-^ (2.10) 



where K is a constant. As usual a denotes the scale factor of the universe, satisfying 
H = a/a and normalised so that a = 1 at the present epoch. We are assuming critical 
density, corresponding to if = 0, so 

= -^p (2.11) 

As explained in Appendix A, the quantity K is related to the intrinsic curvature scalar of 
comoving hypersurfaces by 

= 6^ (2.12) 

For critical density i?^^-* vanishes and the comoving hypersurfaces are flat (the curvature 
scalar determines their curvature entirely, because they are homogeneous). As a result, the 
space coordinates can be taken to be Cartesian. It is convenient to use comoving coordinates 
x*, which are related to Cartesian coordinates by y* = ax*, and are constant along each 
comoving worldline. 

Prom Eqs. ( |2.9[ ) and ( 2.11| ), it follows that during any era when w = p/p is constant, 

H = — - — (2.13) 

3 + 3w ^ ^ 

and therefore a oc t2/{3+3u>)^ -^y-g ^^.^ interested in matter domination, w = 0, and radiation 
domination, w = 1/3. 



2.3 The perturbed fluid flow equations 

Now we consider perturbations away from homogeneity and isotropy. The vorticity is 
expected to be negligible, and can in any case be eliminated by redefining the fluid four- 
velocity as indicated in Section This having being done, comoving hypersurfaces exist. 
On each hypersurfaces one can define average plus a perturbation, 

p(x,t) = />(t) + Mx,t) (2.14) 
p(x,t) = p(t)+5p(x,t) (2.15) 
//(x, t) = H{t) + 5/7 (x, t) (2.16) 

Here t is the time coordinate labelling the hypersurfaces, and x = {x^,x'^,x^) are space 
coordinates. We would like to choose them to be comoving coordinates, related to Cartesian 
coordinates by = arc*, with a the average scale given hy a/a = H (we shall not have 
occasion to define a perturbed scale factor). This cannot be done exactly, because the 
comoving hypersurfaces are not flat in the presence of perturbations. However the departure 
from flatness is of first order, and can therefore be ignored when consider perturbations 
which are themselves of first order. In other words all perturbations 'live ' in flat space. 
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Independent scales 

Each perturbation / can be written as a Fourier series, defined in a comoving box much 
bigger than the observable universe 

5/(x,t) = 5]5/k(t)e*^-^ (2.17) 

k 

The beauty of this expansion is that each Fourier mode propagates independently, as long 
as the cosmological perturbation theory that we are developing here is valid. The inverse 
wavenumber a/k is said to define a scale, which is specified by giving its present value . 
The term scale is, of course, appropriate because a feature — such as a density enhancement 
that will later become a galaxy — with size r is dominated by wavenumbers of order 1/r. 



The evolution of the density perturbation 

Now we derive differential equations for the perturbations. In doing so we have to remember 
that the comoving worldlines are not in general geodesies, because of the pressure gradient. 
As a result, the proper time interval dr between a pair of comoving hypersurfaces is position 
dependent. Its average may be identified with the coordinate time, and then its variation 
with position is given by 

dr f 6p 



dt \ p + pj 

A simple proof of this important formula is given in Appendix A. When considering pertur- 
bations the variation is a second order effect and can be ignored, so that coordinate time 
is identified with proper time. On the other hand, the variation is a first order effect when 
considering the unperturbed (average) value of a quantity, and has to be included. This 
fact was not recognised in Hawking's pioneering development of the fluid flow approach. 

The energy conservation condition Eq. ( ^.5| ) is not affected by the perturbations, but 
the Raychaudhuri equation Eq. (|2.9|) becomes to first order (Lyth & Stewart 1990b) 



^ = + (2.19) 

3 o p + p 

Here is the Laplacian on a comoving hypersurface, given in terms of comoving coordi- 
nates by 

= a-^6'^didj (2.20) 

The energy conservation equation Eq. ( p. 51) and the Raychaudhuri equation Eq. (|2.19D 
determine the evolution of the energy density and the Hubble parameter along each world- 
line, including first order perturbations away homogeneity and isotropy. Averaging them 
over a comoving hypersurface, one verifies that p and H satisfy the unperturbed equations 
Eqs. ( |2.5| ) and (^). Subtracting the averages and using Eq. ( |2.18| ) gives 

{6p)- = -3{p + p)6H - 3Hdp (2.21) 
, , AttG lV^6p 

(6H)- = 2H6H 6p (2.22) 

^ ' 3 3p + p ^ ' 

The dots denote differentiation with respect to time. 
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Remarkably, the pressure perturbation appears in these equations only through its spa- 
tial gradient, because of the relation Eq. ( |2.18| ) between proper time and coordinate time 
on comoving hypersurfaces. For any other choice, such as a synchronous one, there would 
have been additional terms which involve 5P itself, not just its gradient. 



For the Fourier components Eq. (|2.22| ) becomes 



iSH.y = 2H5H, - '-^5p, + V ^ (2-23) 



3 3 \a J p + p 



Eliminating SHk with Eq. (|2.21|) gives a linear second order differential equation for each 
Fourier component of the density perturbation. Often it is convenient to use the density 
contrast 

6 = ^ (2.24) 

P 

and in terms of this quantity the equation is 

H-^^ + [2 - 3i2w - cl)]H-H,, - 1(1 - 6c2 + ^w- 3w')5], = - f 4?^ — (2-25) 

2 \arl J p 

where w = p/ p and = p/ p. This equation was first given (with a different choice of 
variables) by Bardeen (1980). 

Horizon entry 



The right hand side of Eq. ( 2.25| ) involves the ratio aH/k of the scale a/k to the Hubble 



distance H~^. As long as gravity is attractive, so that a < 0, it is a decreasing function 
of time. Each scale is said to enter the horizon at the epoch aH/k = 1. Before horizon 
entry, when aH/k > 1 it is said to be outside the horizon, and afterwards to be inside the 
horizon. As we shall see, the factor aH/k is ubiquitous in cosmological perturbation theory, 
and quite different physics operates in the eras before and after horizon entry. 

2.4 The case of zero pressure gradient 



The right hand side of Eq. (2.25), which involves the pressure gradient, will be negligible well 
before horizon entry unless 5p is far bigger than 5p. As we shall see this seemingly bizarre 
circumstance does occur if the initial conditions are what are called isocurvature, but it 
does not occur with the standard adiabatic initial condition which we shall take for granted 
unless otherwise stated. The pressure gradient is negligible also after matter domination. 



simply because there is no significant source of pressure. In these circumstances Eq. ( |2.25 ) 



can be reduced to a first order equation, which has a very simple interpretation (Lyth 1985; 
Lyth & Mukherjee 1987). 

To derive the solution and its interpretation, it is best to go back to the original fluid 
flow equations Eqs. ( |2.5| ) and ( 2.19] ). They are valid in the presence of first order perturba- 



tions, and we see that if the pressure gradient is negligible they reduce to the unperturbed 
equations. Each comoving region of the universe evolves as if there were no perturbation, 
when the pressure gradient is negligible. As a result, the unperturbed Friedmann equation 
is valid for every comoving region, with a time independent value of K which is still related 
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to the spatial curvature scalar by Eq. (2.12). We are assuming that the average of K on 
each hypersurface vanishes (ie., the average density is critical with respect to the average 
Hubble parameter), so that the Friedmann equation becomes 



-P 



6K 



(2.26) 



This equation is valid to first order in the perturbations, and to that order we can ignore 
the perturbation of a in the last term. The perturbations therefore satisfy 



2H6Hk 



SvrG 



(2.27) 



3 

Combining Eqs. (^ ) and (^^ one finds a first order differential equation for the 



density contrast, 



2H-^ d 
5 + Sw dt 



+ 



2 + 2w 
5 + 3w;' 



(2.28) 



where w = p/p. 

The time-independent quantity T^k on the right hand side of this equation is defined by 



^ 3 5Kk 



(2.29) 



The second equality uses Eq. (|1D, which as we show in Appendix A remains valid in the 
presence of the perturbation. According to this equality 7?.k is a measure of the perturbation 
in the spatial curvature of comoving hypersurfaces. Although TZ is extremely useful in 
cosmological perturbation there is unfortunately no standard notation for it. The first 
author to introduce it seems to have been Bardeen (1980), who called it (j)m- He noted that 
it is constant outside the horizon, but did not relate it to the comoving density perturbation. 
Kodama and Sasaki (1984) called it TZm and noted (in somewhat obscure language) its 
simple relation Eq. ( ^.41]) with the inflaton field perturbation. On scales far outside the 
horizon it coincides with (/3 of Bardeen, Steinhardt & Turner (1983) and the ^ of Salopek, 
Bond & Bardeen (1989). 

During any era when w is constant, Eq. (2.28|) has a solution 



k J ^ 5 + 3w ^ 
The most general solution is this one, plus a solution of the homogeneous equation 

p^^-[(5+3«,)/(3+3«;)] 



(2.30) 



(2.31) 



The latter decays, so that a few Hubble times after the onset of the epoch in question the 
density contrast will be given by Eq. ( 2.30| ). 

Eq. ( 2.31 ) is valid only when the pressure gradient is negligible. Thus it is valid after 
matter domination. 



-7?-k (final) 
5 



(2.32) 
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and it is also valid during an initial era, well before horizon entry and well before matter 
domination, giving 

'ai7\2, 4 , , , 

—j <^k = g^k (initial) (2.33) 

2.5 Relation to the synchronous gauge 

So far we have worked exclusively in the comoving gauge. Many authors work in the 
synchronous gauge described in Section As we now explain, the relation between the 
gauges is very simple during the initial and final eras (Lyth & Stewart 1990b). 

As discussed in Section 2.1, each gauge corresponds to a slicing of spacetime into hy- 
persurfaces, the density perturbation in that gauge being defined by Eq. ( 2.14] ). Choosing 



the averages to be the same, which corresponds to evaluating the perturbations at the same 
coordinate time t, one learns that the synchronous gauge perturbation 5ps is related to the 
comoving perturbation 5p hy 

Sps = 6p + p6t (2.34) 

where 5t is the time displacement going from the comoving to the synchronous hypersurface. 
From Eq. ( ^.18 ) it is given by 



St= -^dt (2.35) 
Jo P + P ^ 

where 6p, like Sp is to be evaluated in the comoving gauge. 

During the initial era, p (x p (x and 5p oc 5 p o(. with a oc t^/^, leading to 
6ps = Thus the synchronous and comoving density perturbations are identical except 
for a scale-independent factor. After matter domination 6p = so that 6t remains con- 
stant, and therefore quickly becomes negligible compared with the timescale t. Thus the 
synchronous and comoving gauges become identical after matter domination. We shall need 
these important facts in what follows. 

2.6 The transfer function 

For scales entering the horizon well after matter domination {k~^ ^ = 80Mpc) the 
initial and final eras overlap so that 7^k(iiiitial) = 7?.k(final). Then, Eqs. ( 2.33| ) and ( |2.32 ) 



determine the final density contrast in terms of the initial one. On smaller scales, there is 
a linear transfer function T{k), which may be defined by 

7^k(final) = T(/c)7^k(initial) (2.36) 

An equivalent, and more usual, definition is 

a~"^<5k|final = AT{k)5k\i^itia.l (2.37) 

where the (time dependent) right hand side is evaluated at an arbitrarily chosen time during 
the initial era, and the constant A is chosen so that T becomes equal to 1 on large scales. 
From the results of the last section, this second definition can be used also to define a 
synchronous gauge transfer function, which is identical with the comoving one. 

To calculate the transfer function, one needs an initial condition specifying the relative 
abundance of the neutrinos, photons, baryons and cold dark matter long before horizon 
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entry. The most natural condition, which we adopt in this subsection, is that the abundances 
of all particle species are uniform on hypersurfaces of constant total energy density. This 
is called the adiabatic condition, for a reason that will be clear in a moment. It is indeed a 
reasonable one, because at very early times a given comoving length scale far exceeds the 
Hubble distance, over which causal processes can operate on the Hubble timescale. On the 
other hand, as we discuss in Section ^.7, it is not absolutely mandatory. 

Adopting the adiabatic initial condition, consider a hypersurface displaced from a hy- 
persurface of constant density by time 6t(x). The density perturbation of any species X is 
given by 5px = px^t, so the perturbations of two species X and Y are related by 

^ - ^ = (2.38) 

PX PY 

To first order in the perturbations, px and py may be taken to have their unperturbed 
values (their averages over a comoving hypersurface). They each therefore satisfy the energy 
conservation condition Eq. ( p.Sp , 

Px = -^H{px+Px) (2.39) 

It follows that each species of radiation has a common density contrast 6r, and each species 
of matter has a common density contrast 5^, with the relation 

Sm -l^r = (2.40) 

This relation holds in both the comoving and synchronous gauges. 

Given the adiabatic initial condition, the transfer function is determined by the physical 
processes occurring between horizon entry and matter domination. If the radiation behaves 
as a perfect fluid, its density perturbation oscillates during this era, with decreasing am- 
plitude. The matter density contrast living in this background does not grow appreciably 
before matter domination because it has negligible self gravity. The transfer function is 
therefore given roughly by 

T{k) =1 (fc < /ccq = SOMpc^i) 

T{k) ={k,Jkf {k>ke^) (2.41) 

An accurate description of the evolution under this perfect fluid assumption is provided 
by Eq. ( ^.25 ), together with Eq. ( 2.52| ) of Section §.7 which describes the evolution of the 



entropy perturbation. Equivalent equations exist in the synchronous gauge, and they were 
integrated numerically by Peebles (1982a) (without including the neutrinos) to give the first 
estimate of the transfer function. 

The perfect fluid description of the radiation is far from correct after horizon entry, 
because roughly half of the radiation consists of neutrinos whose perturbation rapidly dis- 
appears through free-streaming. The photons are also not a perfect fluid because they 
diffuse significantly, except between the Silk scale A;^^ ~ 1 Mpc and the horizon scale at 
decoupling k^^ = 180 Mpc. One might therefore consider the opposite assumption, which 
is that the radiation has zero perturbation after horizon entry. Then the matter density 
perturbation evolves according to Eq. ( |2.25| ), with 6 and p now referring to the matter alone, 

6 + 2H6 - AttGp5 = (2.42) 
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(This equation follows also from Newtonian physics, as is clear from the fact that the fluid 



flow equations Eqs. (2^) and ( |2.6D coincide with the Newtonian ones for p = 0.) The 
growing solution is 

6 = A{1 + 1.5y) (2.43) 



and the decaying solution is 



6 = B{{l + 1.5y) 



In 



' (1 + ^)^/^ + 1 ^ 
(l + y)V2_l 



3(1 + y) 



1/2 



(2.44) 



where A and b are constants and y = a/ueq- Even the growing solution increases only by a 



factor of order 2 between horizon entry (when Eq. ( 2.42 ) first becomes valid) and the epoch 



tteq when matter starts to dominate. The transfer function is therefore again given roughly 



by Eq. (|2^ ). 

A more sophisticated procedure (Blumenthal & Primack 1984; Primack & Blumenthal 
1984) is to keep both modes, matching to the (approximately correct) time dependence 
5 oc {aH)~^ at horizon entry, and one then finds a somewhat bigger transfer function on 
small scales, with Eq. ( p.41| ) replaced by k~^lnk. 

Since the radiation consists of roughly half neutrinos, which free stream, and half photons 
which either form a perfect fluid or only diffuse, neither the perfect fluid nor the free- 
streaming approximation looks very sensible. A proper calculation should presumably take 
quantitative account of the following effects 

• Neutrino free-streaming around the epoch of horizon entry. 

• The diffusion/free-streaming of photons around the time of horizon entry (except on 
scales well above the Silk scale ^ 1 Mpc and well below the horizon scale at 
decoupling = 180 Mpc). 

• The diffusion of the baryons along with the photons. 

• The establishment after matter domination of a common matter density contrast, as 
the baryons fall into the potential wells created by the cold dark matter. 

To first order in the perturbations, all of these effects apply separately to each Fourier 
component, so that a linear transfer function is indeed implied. Excellent accounts of the 
standard (synchronous gauge) formalism for calculating them have been given by Peebles 
(1980) and Efstathiou (1990). 

To date, all numerical calculations of the cold dark matter transfer function have been 
done in the synchronous gauge. They include some, but not necessarily all of the effects 
mentioned above. Unfortunately, none of them is fully documented, and the parametrisa- 
tions which are supposed to fit them do not agree very well with each other. We use one of 
the parametrisations given by Bond and Efstathiou (1984), which has been widely used in 
the literature (Efstathiou 1990) 



T{k) = [l + [ak + {bkf^^ + {ckf 



-Ijv 



(2.45) 
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where v = 1.13 and 



a = 6A{Qrnhr^ h'^Mpc (2.46) 
b = 3 {Qmhy^ h-^Mpc (2.47) 
c = 1.7{Qmhr^li-^Mpc (2.48) 

We postpone until Section |^ a discussion of the dependence of this expression on h and the 
present matter density Qm, taking until then the canonical values = 1 and h = 0.5. 
The transfer function is shown in Figure la, where one sees that the behaviour estimated 
in Eq. (2.41) is roughly correct, though the break at is not at all sharp. The transfer 



function is shown right down to the scale k~ = 10~ Mpc, corresponding to the mass 
M ~ W^Mq enclosed by a dwarf galaxy, but as we shall see it has so far been compared 
with observation only for > 1 Mpc, corresponding to the mass M ~ W^'^Mq enclosed 
by a bright galaxy. 

Another widely quoted parametrisation is that given by Bardeen et al (1986), based on 
a calculation of Bardeen (1986). It is compared in Figure lb with the one we are using, 
along with four other parametrisations. Not all the parametrisations attempt an accurate 
reproduction down to the smallest scales. One sees that even on the presently relevant 
scales > 1 Mpc the fits range more than 10% each side of the one we are using. Adams 
et al (1993) claim that on the scale k^^ ~ 10 Mpc the parametrisation of Bardeen et al 
(1986) is extremely accurate for = 0, and that increasing $7^ to the value — -06 
favoured by nucleosynthesis decreases it by about 15%. One sees from Figure la that this 
would take it about 5% below the parametrisation that we are using. 

In view of the crucial importance of the cold dark matter transfer function a fully 
documented calculation seems to be called for, together with a parametrisation which is 
accurate to 1% or so for all in the range to .1, and all scales of cosmological interest. 

The transfer function, which encodes the solution of linear equations, ceases to be valid 
when the density contrast becomes of order 1. After that, the highly non- linear phenomenon 
of gravitational collapse takes place, as discussed in Section 



2.7 Isocurvature density perturbations? 

Instead of the adiabatic initial condition, one could consider what is called the isocurvature 
initial condition. According to this condition there is no perturbation in the total energy 
density, or in the expansion rates of the separate components (and therefore none in the 
overall expansion rate H). There are however perturbations in the energy densities of two 
or more of the components, which add up to zero. 

Since there is no change in the energy density or the expansion rate the Friedmann equa- 
tion is unaffected, and the comoving hypersurfaces do not receive a curvature perturbation. 
This is the reason for the term 'isocurvature'.^ 

*The specific reference to comoving hypersurfaces is not usually spelled out in the literature, and one 
might form the impression that it is the curvature of space-time which is supposed to be unaffected. Of 
course that is not so, because in general isocurvature perturbations affect the pressure which is related 
through the field equation to the space-time geometry. In particular, the curvature scalar of space-time is 
R= -87rG(p-3p). 
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The most general isocurvature density perturbation may be specified by giving the 
perturbation Sx in the ratio of the the number density of each species X, to (say) the 
number density of the photons. (An equivalent definition is that —Sx is the perturbation 
in the photon entropy per particle.) Except for particular epochs, each species constitutes 
to a good approximation either radiation, p = p/3, or matter, p = 0. Its energy density in 
these cases evolves with time like its number density to the powers 4/3 and 1 respectively, 
so one has 

Sx = l^x - j^'y radiation (2.49) 

Sx = — jS-y matter (2.50) 

Any initial energy density perturbation may be expressed as the sum of an adiabatic per- 
turbation and an isocurvature perturbation. 

We shall focus exclusively on a perturbation of the cold dark matter density. Thus we 
assume that initially there is a common radiation density contrast 6r, a baryon density 
contrast 6b = |(5r and a cold dark matter density contrast specified by 

S = 6c- ^6r (2.51) 

To our knowledge, an isocurvature perturbation of the baryon density has not been discussed 
in the context of cold dark matter, and an isocurvature perturbation of the neutrino density 
has not been discussed at all. In the context of cold dark matter this list is exhaustive, 
since we have agreed to define the perturbations relative to that of the photon density. 

Before seeing how to handle an isocurvature cold dark matter perturbation, we briefly 
discuss its theoretical motivation. In the context of inflation, an adiabatic perturbation of 
some magnitude is mandatory, because it corresponds to a quantum fluctuation of the in- 
flaton field. An isocurvature perturbation, on the other hand, can originate during inflation 
only through the quantum fluctuation of some other field. The resulting inhomogeneity in 
the field has to survive after inflation, until some mechanism converts it into an isocurvature 
density perturbation. 

Until recently it seemed likely that the axion field, if it exists, could easily generate an 
isocurvature perturbation^ through this mechanism.^ Indeed, if a quantum fluctuation of 
the axion field survives to the end of inflation, and if in addition the reheat temperature 
after inflation is too low to restore the Peccei-Quinn symmetry, then the fluctuation will 
survive until the axion acquires a mass at the epoch T ~ Aqcd ~ 100 MeV, to become an 
isocurvature axion density perturbation. The condition that symmetry is not restored after 
inflation is that the reheat temperature is less than the axion decay constant fa ~ lO^'^ to 
10^'^ GeV, which might well be satisfied. However, it has recently been realised (Lyth & 
Stewart 1992b) that the fluctuation may not survive until the end of inflation, essentially 

^ Since the axion inhomogeneity generated through this mechanism is not necessarily small, a significant 
isocurvature cold dark matter perturbation can be generated for any axion density within the cosmologically 
allowed regime 10~* < < 1, even if something else is responsible for the bulk of the cold dark matter. 
The perturbation is however non-Gaussian if the axion inhomogeneity is big. 

®The status of the subject in 1989 is described in the texts of Linde (1990) and Kolb and Turner (1990). 
Developments since then can be traced from the papers of Lyth (1992), Lyth and Stewart (1992b) and Lyth 
(1993a). 
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because symmetry can be restored by the Hawking temperature. The condition for this to 
happen turns out to be Hi > fa, where Hi is the Hubble parameter when the observable 
universe leaves the horizon during inflation. According to Eq. ( 5.46| ) this condition is met 
in many models, so that an isocurvature density perturbation arising from the axion is far 
from inevitable. Nor is there any other well motivated way of generating an isocurvature 
cold dark matter perturbation. 

Despite the lack of motivation for isocurvature initial conditions, they are worth dis- 
cussing if only to see whether they can provide a viable explanation of the observed large 
scale structure, without violating the microwave background anisotropy measurements. Ac- 
cordingly, the evolution of an isocurvature cold dark matter density perturbation is de- 
scribed now, and its effect on the microwave background is described in Section ^. The 
conclusion is that such a perturbation cannot alone be the cause of large scale structure. 



The two fluid formalism 

The evolution of the density perturbation before horizon entry can be described by ignoring 
the baryons, and treating the radiation and cold dark matter as a pair of uncoupled fluids. 
Using the metric perturbation formalism in the synchronous gauge, the resulting rather 
complicated equations have been given by several authors, including Efstathiou and Bond 
(1986) who solved them numerically. We shall work with the fluid flow formalism in the 



comoving gauge, where the equations are simply Bardeen's equation Eq. (2.25), together 
with an equation tracking the evolution of the entropy perturbation which was first given 
by Kodama and Sasaki (1984, 1987). The equation is 



H-^Sk + (2 



3c2)F-25k 



{an) 



-c^'Sk + ^(l + it')"^<5k 



(2.52) 



where = [3 -|- 4(acq/a)] ^. Evaluating dP in terms of 6 and S, Bardeen's equation 
becomes 



H-'^'S], + [2 - 3(2^1- - ci)]H-^6], 

2 



c% + 3 



1 + (a/Oeq) 



l + (4/3)(a/acq)- 



(2.53) 



Eqs. (2.52) and ( p. 53 ) completely describe the evolution of the two fluids. They have 
four independent solutions, which as we now see correspond at earlier times to growing and 
decaying modes, each with either adiabatic or isocurvature initial conditions. 



The initial era 

At an initial epoch, well before horizon entry and well before matter domination, there are 
the following independent modes (Kodama & Sasaki 1984, 1987). (The subscripts k are 
dropped for clarity) 

Adiabatic growing mode 



S/5 = O 



a / a \'' 



(2.54) 
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Adiabatic decaying mode 

6 = Ba-^ S/6 = O 
Isocurvature growing mode 

S = C 6/S = 



a / a Y 
acq \kH ) 



a \ /ay 
acqi \kH ) 



Isocurvature decaying mode 

S = D\n{a/aeq) 6/S = 



{kH) 



(2.55) 



(2.56) 



The growing adiabatic mode corresponds to Eq. ( ^.30| ), the constant A being related to 
the constant initial value of TZ. The decaying adiabatic mode corresponds to Eq. ( 2.31| ), 
with TZ = 0. We dropped it and will do the same with the decaying isocurvature mode. 
Thus, we are interested in the growing isocurvature mode|] 

The large scale radiation density perturbation 

After the initial epoch the evolution of the matter and radiation densities can be followed 
by solving Eqs. (2.52) and ( p. 53 ), provided that diffusion and free streaming are ignored. 
This is certainly permissible before horizon entry, and in fact an analytic solution of the 
equations is available for that case (Starobinsky & Sahni 1984; Kodama & Sasaki 1987). 
It shows that the entropy perturbation S is practically constant, corresponding to the fact 
that particles do not have time to flow on a given scale before horizon entry. It also shows 
that \5/S\ remains small. Substituting 6pc — —6pr into the definition of S, this implies 
that during radiation domination 6r <C Sc and 5 ~ Jo whereas during matter domination 
5c <S Sr and S ~ — |(5r-. Since (xT^, this implies that after matter domination the photon 
temperature has a perturbation 

f = -\S. (2.57) 

We shall discuss the resulting microwave background anisotropy in Section 

The isocurvature transfer function 

Well after matter domination the pressure perturbation is negligible, so that the density 
contrast has the behaviour Eq. ( 2.32| ). One can then define a transfer function. 



(f) 



5p^ 
P 



7iso(fc)5'k(initial) 



(2.58) 



final 



^By considering the separate fluid flow equations, one can show that the decaying modes correspond to 
relative motion between the fluids, which is necessarily absent well before horizon entry (Lyth & Stewart 
1992b). In the case of an isocurvature mode this is a non-trivial result, because it justifies the standard 
assumption that all of the initial S value feeds into the growing mode. It is this initial value which is related 
to the inflationary potential, so for a given potential the predictions would be altered if a decaying mode 
were initially present, even though it would become negligible long before the present. 
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The subscript 'initial' refers to the era well before horizon entry, when S is practically 
constant, and the subscript 'final' refers to the era well after matter domination. This 
transfer function is the same in the comoving and synchronous gauges, because S is the 
same in both gauges and so is the left hand side. 

For scales entering the horizon, free streaming and diffusion are negligible so that the 
transfer function from the analytic solution of Eqs. ( p. 52 ) and ( p. 59 ) which has already 
mentioned. This gives (Kodama & Sasaki 1987) 

Ti,a{k) = 2/15 (fc-i » k-^) (2.59) 

Notice that according to this expression |5k/'S'k| remains small until horizon entry, in accor- 
dance with what was said earlier. 

The situation on smaller scales is similar to the one that we described before for the 
adiabatic case. To estimate the transfer function roughly, one notes that 5m ^ S at horizon 
entry and that there is little subsequent growth of 5m until matter domination. As a result 



the transfer function is again given roughly by Eq. ( 2.41 ). For a more accurate calculation 
one could assume that the radiation free streams away promptly, or go to the other extreme 
of taking it to be a perfect fluid. Neither approximation looks very sensible, since the 
radiation consists of roughly half neutrinos which free stream, and half photons which 
either are a perfect fluid or else only diffuse. Efstathiou and Bond (1986) find however 
that the perfect fluid approximation is valid to 5%. We are not aware of any simple reason 
why the effect of free-streaming and diffusion should be so small; although the radiation 
density contrast relative to the matter density contrast is far smaller, at horizon entry, 
in the isocurvature case than in the adiabatic case (about the same instead of a factor 
Pr/Pc bigger), it is still the thing that determines the shape of the transfer function in 
the vital region around the break at A: = fceq. One would like to know if the perfect fluid 
approximation is also accurate for the adiabatic case. 

2.8 The peculiar velocity field 

After matter domination it is very useful to introduce the concept of a peculiar velocity 
field. Since the comoving and synchronous gauges coincide in this regime, the following 
treatment is valid in both of them. 

In any small region around a comoving observer, the receding fluid defines a three- 
velocity field u^. Its components are to be defined in terms of locally inertial coordinates, 
in which the observer is instantaneously at rest. By a 'small region' is meant one which is 
small compared with the Hubble distance, so that the curvature of space-time is negligible. 
Ignoring the perturbations, 

= Hy' (2.60) 

where the observer is taken to be at the origin of the spatial coordinates y'. Including the 
perturbations, one defines a peculiar velocity field by 

u' = Hy' + v' (2.61) 



where H is the average defined in Eq. (2.16) 
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The three-velocity field cannot be defined over an extended region because of space- 
time curvature, but one might still hope to define a peculiar velocity field, such that 
Eq. ( 2.61| ) is satisfied in the small region around each comoving observer. The field would be 



a four-vector, everywhere orthogonal to the comoving four-velocity and with magnitude 
much less than 1. As discussed in Section ^.3, such a globally defined peculiar velocity field 
does not strictly speaking exist because of gravitational waves. The gravitational waves can 
however be handled separately and it is legitimate to ignore them in what follows. 

Like any perturbation, the peculiar velocity field 'lives' in flat space, described by co- 
moving coordinates related to Cartesian coordinates by = ax^. It can be decomposed 
into a longitudinal (irrotational) part and a transverse (rotational) part, 

V = Vtran + Vlong (2.62) 

where the transverse part satisfies V.vtran = and the longitudinal part is the gradient of 
some velocity potential. For a Fourier mode with k along the z direction, the longitudinal 
part has only a z component and the transverse part has only x and y components. 

The transverse part contributes to the vorticity uJij, and is the quantity that we agreed 
earlier to subtract off so that comoving hypersurfaces can be defined. It decays with time 
like a-^{p + Py and can be ignored (Weinberg 1972; Peebles 1980; Lyth 1993b). 

The longitudinal part of the peculiar velocity is directly related to the density perturba- 
tion. After matter domination, which is the only regime of interest in the present context, 
the relation reduces to a well known Newtonian expression. The same expression follows 
from the present analysis, because the fluid flow equations on which it is based reduce to the 
Newtonian ones during matter domination. To derive it, start with the result V.v = 35H 
which follows from Eq. (|3). From Eqs. (|2T^ ), (|237D , (|]2|) and ( pO| ), it follows that 

V.v = {A7:G6p)t (2.63) 

The solution of this equation is 

V = -tVi) (2.64) 

or 

where 

V'(x,t) = / ^4^^'^' (2-66) 



|x' — x| 

The factor converts coordinate distances into physical distances. Since it is related to 
the density perturbation by the Newtonian expression, ■0 is called the peculiar gravitational 



potential. It is independent of t because 5p oc a^, and from Eq. ( 2.32 ) we see that it is 
related to the spatial curvature perturbation by 

= -^7^(final) (2.67) 





The peculiar velocity on the other hand is proportional to a 



1/2 



27 



Prom Eqs. ( |2.64| ) and ( p.66| ) the Fourier components of v, ^/^ and S are related by 

* = (2.69) 

It follows that on scales much bigger than the horizon they satisfy 

|<5k| < |vk| < iV^kl (aH/k > 1) (2.70) 
For an isocurvature perturbation, it follows from Eqs. ( ^.58| ), (|2.59| ) and ( ^.69| ) that 

V'k = --5k(initial) {aH/k > 1) (2.71) 
5 

3 The spectrum of the density perturbation 

During the linear regime one expects that the perturbations will be random fields with 
simple stochastic properties. For a Gaussian field these are specified completely by the 
spectrum. 

3.1 The spectrum and stochastic properties 

We are interested in the density contrast 6 and the quantity TZ which measures the associated 
spatial curvature perturbation, and after matter domination we are also interested in the 
components Vi of the peculiar velocity field, and the peculiar gravitational potential ■0. Let 
/ denote any one of these perturbations. Its most important statistic is its spectrum, which 
is essentially the smoothed modulus-squared of its Fourier coefficient. To be precise, the 
spectrum may be defined as the quantity 



Lk 



3 



rf ^ j 47r(|/k|^) (3.1) 

where L is the size of the box for the Fourier expansion Eq. ( p.l7| ), and the bracket denotes 
the average over a small region of A;-space. The normalisation is chosen to give a simple 
formula for the dispersion (root mean square) of /, which we shall denote by crj. From 
the Fourier expansion one has = I /k I ' ^^'^ since the possible values of k form a cubic 
lattice with spacing 2tt/L the transition from sum to integral is 

The dispersion (Tj is therefore given by 

a}^{f{^))= frfik)^ (3.3) 



J2 — ^ 47r / k'^dk (3.2) 
k 
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with the brackets now denoting the spatial average. For the density perturbation / = 5 it 
is useful to define the correlation function ^(r) by 

= (/(r + x)/(r)) = fVfikf^^f (3.4) 

The analogous quantity is useful for other perturbations like the peculiar velocity compo- 
nents, though it is not then called the correlation function. For r = it clearly reduces to 
a}. 

If the phases of the Fourier coefficients are random, / is said to be Gaussian, and then 
all of its stochastic properties are determined by its spectrum. In particular the probability 
distribution of /, evaluated at randomly chosen points, has a Gaussian profile. As long as 
the perturbations are evolving linearly this Gaussian property is implied by the inflationary 
models to be discussed in Section |5|, and we take it for granted from now on. 

In linear theory the spectra of all of the perturbations are {k dependent) multiples of that 
of the spectrum Vs of the density contrast. In the literature, our Vs is denoted variously by 
Vp (Salopek, Bond & Bardeen 1989), Pp (Lyth & Stewart 1990b), (Kolb & Turner 1990; 
Taylor & Rowan-Robinson 1992)), 5pk/p (Linde 1990) and dap/dlnk (Bond & Efstathiou 
1991). An alternative definition, used especially in the older literature and denoted by P, 
is equal to Nk~^Vs, where is an author-dependent normalisation factor which is often 
left undefined. 

Prom Eqs. ( 2.32 ) and ( ^.36| ), the spectrum of the density contrast after matter domina- 



tion may be written 

Vs{k) = T\k)6l{k) (3.5) 

The quantity Sh specifies the initial spectrum.^ In fact, from Eq. ( ^.33 ) it is related to the 
spectrum of the initial curvature perturbation TZ by 

= ^^7^(fc) (3.6) 

The subscript H has been used because 5jj is exactly equal to the value of Vs on horizon 
entry on scales k~^ » k~^, and approximately equal to it on smaller scales. As we shall 
see in Section ^.2, this means that dnik) is roughly equal to the mean square value of S at 
horizon entry, for a density enhancement with comoving size of order k~^. 

The standard assumption is that Sfj^ is independent of k, which with the older definition 
of the spectrum corresponds to P (x k. A more general possibility is to consider a spectrum 

4 oc (3.7) 

corresponding to P oc /c". The exponent n is called the spectral index. The standard choice 
of n = 1 was first advocated by Harrison (1970) and Zel'dovich (1970) on the ground that 
it is the only one making the perturbation small on all scales, at the epoch of horizon entry. 
On the other hand, the power-law dependence on n need only be an approximation, valid 
over a limited range of scales. Prom this viewpoint the value n = 1 is not so special, and the 

*For n = 1 it is equal to the r/n of Scaramella and Vittorio (1990) and to in times the en of Abbott and 
Wise (1984c). 
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possibility has recently been explored that n is less than 1, tilting the Harrison-Zel'dovich 
spectrum to give more power on large scales. Of course there is no reason in principle why 
the spectrum should have a power-law dependence at all; the effective value of n could 
change with scale as one goes from the regime ~ 1 Mpc explored by the small scale 
galaxy correlations to the regime ~ 10^ Mpc explored by the large scale cmb anisotropy. 

In Section ^ we shall see how inflationary models of the early universe inevitably give an 
adiabatic density perturbation at some level. For a given inflationary model its spectrum 
can be calculated in terms of the parameters specifying the potential. In typical inflation 
models there is a quite accurate power-law dependence over the above range of scales, the 
spectral index being determined by two small parameters ei and rji, which we shall define 
in Section ^, 

n = l + 2r]i- 6ei (3.8) 

The precise values of ei and rji depend on the infiationary model. In many models they are 
extremely small, making n very close to 1. In others this is not so, and we shall describe 
in Section ^.4 models which span the range < n < 1.2. The prejudice towards negative 
values is a reflection of the fact that ei is positive by deflnition. 

The isocurvature spectrum 

In the case of an isocurvature initial condition, Eq. ( p.5| ) is replaced by 

Vsik) = ^-^y TUk)Vs{k) (3.9) 

At least over a limited range one can hopefully define an isocurvature spectral index by 
Vs{k) oc fe^is^. Since the adiabatic and isocurvature transfer functions have a similar shape, 
a spectrum Vs of given shape can arise either from an adiabatic perturbation with spectral 
index n or an isocurvature one with spectral index niso = n — 1. 

If the isocurvature perturbation originates as a vacuum fiuctuation during infiation its 
spectral index is given by 

niso = -2ei (3.10) 

Many of the usually considered inflationary models require ei to be very small, leading to 
a very flat spectrum, but others allow it to be significant, tilting the spectrum to give more 
large scale power. Since e is positive by definition it is impossible (not just difficult, as in the 
case of an adiabatic perturbation) to go the other way, at least within the almost universally 
accepted slow-roll paradigm. (This restriction applies, of course, only to an isocurvature 
perturbation generated in the manner we have described. An isocurvature perturbation 
produced at a phase transition will on the contrary have most of its power on the current 
horizon scale, and so could form very small scale early structure. An example might be 
the inhomogeneity of the axion density associated with domain walls present at the epoch 
T ~ 1 GeV (Hogan & Rees 1988).) 

3.2 The filtered density contrast 

The density contrast 6{x) will evolve linearly as a <C 1, except in those rare regions where it 
becomes > 1 and gravitational collapse takes place. Following the usual practice, we assume 
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when necessary that the hnear evolution is at least roughly valid right up to the epoch cr ~ 1. 
Soon after that epoch, a large fraction of the matter collapses into gravitationally bound 

objects, and linear evolution becomes completely invalid. For definiteness let us define the 
epoch of non-linearity as precisely the epoch a = 1. Since 6 cx a after matter domination, 
it corresponds to a redshift z^i given by (1 + z^i) = l/o"o, where ao is the linearly evolved 
quantity evaluated at the present time. 

A 'filtered' density contrast 5(i?/,x), which has a smaller dispersion cr(i?/), can be 
constructed by cutting off the Fourier expansion of 5(x) above some minimum wavenumber 
~ or equivalently by smearing it over a region with size — Rf- The filtered quantity 

will evolve linearly until a later epoch u{Rf) = 1. 

A precise definition of the filtered density contrast is made by means of a 'window 
function' W{Rf,r), which is equal to 1 at r = and which falls off rapidly beyond some 
radius Rf (Peebles 1980; Kolb &: Turner 1990). The filtered density contrast is 



(5(i?/,x) = J W{Rf,\x' -:x.\)6{x')d^x' (3.11) 



and its spectrum is 



l2 



where 



and 



The filtered dispersion is 



Vs{Rf,k) = [W{Rf,k)/Vf\ Vs{k) (3.12) 

W{Rf,k) = J e-'^''W{Rf,r)d^x (3.13) 

Vf = J W{Rf,r)d^x (3.14) 

roo |- , -,2 f]]f 

^\Rf) = l [w{Rf,k)/Vf] Ps{k)- (3.15) 

The quantity Vf is the volume 'enclosed' by the filter. It is convenient to define the 
associated mass M = poVf, where po is the present mass density. One normally uses M 
instead of Rf to specify the scale, writing S{A4,x) and a{M). 

The two popular choices are the Gaussian filter 

W{Rf,r) = cxp{-r^/2Rj) (3.16) 

Vf = {2T:f/'^R) (3.17) 

W{Rf,k)/Vf = exp{-kRf) (3.18) 

M = 4.36 X W^^h^iRf/lMpcfMQ (3.19) 

and the top hat filter which smears uniformly over a sphere of radius Rf 

W{Rf,r) = e{r-Rf) (3.20) 

Vf = ATTRj/3 (3.21) 

M = 1.16 X 10^'^h'^{Rf/lMpcfMQ (3.23) 
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The Gaussian filter is the most convenient for theoretical calculations, but the top hat filter 
is widely used to as a means of presenting data. 

Provided that the spectrum is increasing, say like A;™ with m not too different from 1, 
the filtered dispersion will be given roughly by 

a\Rf) ~ V,{k = Rf) (3.24) 

Furthermore, as we shall see in Section ^, the filtered dispersion is roughly the rms density 
contrast of structures of comoving size Rf. Thus the square root of the spectrum of the 
density contrast at wavenumher k is typically a measure of the average density contrast 
within an overdense region of comoving size Rf ~ k~^. This rule holds in the regime where 
the spectrum is increasing reasonably strongly, which corresponds to k~^ > .1 Mpc. 

We have focussed on the density perturbation, but any other perturbation may be 
filtered in the same way. Thus, the entire cosmological perturbation theory can be applied 
to perturbations on any chosen scale, and remains valid for longer the larger that scale. For 
filtering scales Rf 3> 10 Mpc, it is still valid at present. 



3.3 The conventional normalisation 



For a given spectral index, the shape of the spectrum after matter domination is determined 
by the transfer function. It remains to specify its normalisation. In the past this has 
usually been done by comparing the linearly evolved theory with observations of the galaxy 
correlation function on a scale of order 10/i~^ Mpc. Two standard prescriptions exist (Davis 
& Peebles 1983). The first requires that the dispersion crg(r) of galaxy counts in spheres of 
radius r has the observed value unity at r = 8/i~^ Mpc. The second requires that J3{r), the 
integral of the second moment of the galaxy correlation function up to distance r, has the 
observed value 270 Mpc~^ at r = lOh'^ Mpc. 

These normalisation schemes both refer to the statistics of the clustering of optical 
galaxies, not to the clustering of mass as described by the spectrum Vs- It is an important 
ingredient of the CDM cosmogony that these are not necessarily taken to be the same — 
so-called biased galaxy formation (Bardeen et al 1986). To be specific, the mass density 
correlation function ^ given by Eq. ( |3.4| ) is assumed to be related to the correlation function 
^gg of luminous galaxies by 

^99(0 - b'iir) (3.25) 

with a bias parameter b which is roughly independent of the scale r. Including the bias 
parameter, the normalisations implied by the two schemes are 



al{r)/b' 
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V^{k) 



{l/bf 



' sin kr 
(kr)^ 
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cos /crl ^ dk 
(kr)'^ J k 
8h-^ Mpc 



(3.26) 



and 



2r-^Js{r)/b'^ 



2Vo{k) 



{.73/bf 



1 

kr 
at r 
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sin kr 



cos kr 
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Wh-^ Mpc 



dk 
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(3.27) 
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In each case, Vq is the linearly evolved spectrum, evaluated at the present epoch. 

One sees that ag{r)/h is the dispersion of the hnearly evolved density contrast with top 
hat filtering on the scale r, whereas {2r~^ J^{r)Y^'^ /h is its dispersion on the same scale with 
a different choice of filter. In what follows we always adopt the ag normalisation scheme, in 
the sense that normalisation of the spectrum is specified by quoting the quantity b defined 
by Eq. ( 3.26| ), which will be denoted by ftg- In other words, the normalisation is defined by 



giving the quantity 

1/68 ^cfs = aTH(8/i"^ Mpc) (3.28) 

where cjth is the present value of the linearly evolved dispersion, with top hat filtering. 
Only occasionally do we invoke the much stronger assumptions that the galaxy correlation 
function is given by S,gg{r) ~ b'^^{r), with b scale-independent and with ^(r) given by linear 
theory on scales r > 10 Mpc. 

As noted earlier, the dispersion of a filtered quantity evolves linearly only as long as it is 
< 1, and this condition is only marginally satisfied by the quantities encountered in the two 
conventional normalisation schemes. A^-body simulations show that in practice though one 
doesn't do too badly as the nonlinear corrections shrink the comoving size of the fluctuation 
so that it doesn't contribute fully to the filtered variance. Evaluating the integrals, one finds 
that the J3 normalisation is typically 10% lower than the normalisation using ag, across 
the range of n in which we shall be interested. These are shown in Figure 2. As remarked 
above, different choices of transfer function parametrisation can give a correction of around 
10% as weh. 

Figure 3a shows 'Po{k)/a% for a selection of n values, where as always Vo is the linearly 
evolved value of Vs at the present epoch. The spectra for different n typically cross at a 
scale ~ 10 Mpc. Fi gure 3b shows the corresponding dispersion a{M)/as for n = 1 and 
n = 0.7, for both Gaussian and top hat filtering. The mass scale M runs from M = W^Mq, 
the Jeans mass at decouphng (Peebles 1980; Kolb & Turner 1990) to M = lO^^M©, the 
mass of the universe on the scale of the CfA survey. The comoving wavenumber runs 
over the corresponding range of filtering radius Rf. 

For n = 1 the spectrum Vs{k) is practically fiat on small scales. For n < 1 the power on 
small scales is reduced (though even for n = not by as much as in the standard hot dark 
matter model (White, Frenk & Davis 1983) or with string-seeded hot dark matter (Albrecht 
& Stebbins 1992b)). As a result the spectrum starts to fall on scales k~^ < .IMpc. The 
implication of this change for small scale structure is discussed in Section |8[ 



4 The cosmic microwave background anisotropy 

So far we have focussed on the matter perturbation, which can be observed through the 
distribution and motion of the galaxies. In addition there are the perturbations in the cos- 
mological neutrinos and photons, where 'cosmological' denotes particles originating in the 
early universe as opposed to particles emitted later from astrophysical sources. The cosmo- 
logical neutrinos are completely unobservable, but the cosmological photons are observed 
as the cosmic microwave background radiation, variously abbreviated as the cmb, mbr or 
cbr. 
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4.1 The spectrum of the cmb 

Discounting the possibility of early re-ionisation,^ the cmb last scattered at the epoch of 
decoupling. Its surface of last scattering therefore lies practically at the particle horizon 
whose comoving distance is x = 2Hq^ = 1.2 x 10^ Mpc. At this surface an angle 6 degrees 
subtends a comoving distance 

x~ 200 Mpc (4.1) 

The intensity of the cmb as a function of frequency is a very accurate blackbody distri- 
bution (Mather et al 1990), with a temperature which is almost independent of direction, 
T = 2.736 lb .017i^. Here we want to discuss the small anisotropy of the cmb. What one 
will actually calculate and measure is the variation in intensity at fixed frequency, but it 
can be expressed as an equivalent variation AT in the temperature^ 

The temperature as a function of direction e may be expanded in multipoles, 

r(e) =T + v.e + Ar(e) (4.2) 

where 



AT 



oo +/ 

Y: E ^TYrie) (4.3) 



1=2 m=-l 

The dipole term v.e is well measured, and as we shall see v can be ascribed to our motion 
at least if the isocurvature perturbation is negligible. 
The mean square anisotropy over the whole sky is 

(4.4) 

The multipoles, and therefore the mean square anisotropy, depend of course on the position 
of the observer. For a randomly placed observer 
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position 1=2 



-E(2/ + 1)S? (4.5) 



where 



S/ = (|arr) (4.6) 

\ / Dosition 
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position 

As we shall see in a moment, the anisotropy can be calculated in terms of the density 
perturbations of the matter and radiation (together with any cosmological gravitational 



®Very early structure could cause early re-ionisation, which moves the surface of last scattering closer. 
On large angular scales this has little effect on the anisotropy, but it could be significant for scales less than 
a few degrees. As we shall see in Section ^, if the spectrum of the density perturbation is even roughly 
flat the upper limit on its normalisation implied by the large scale anisotropy prevents such early structure 
from forming. Reionisation is more likely in the case where there are nongaussian seeds, which can favour 
sharp overdensities. An example of this is the texture scenario (Turok & Spergel 1990), though it appears 
to be in trouble primarily from the COBE observation, comparison with which is effectively unaffected by 
the possibility of reionisation (Pen, Spergel & Turok 1992). 

^"As we shall see, AT/T is less than 10~^; by specifying such a number one does not imply that the 
blackbody form is so accurate, still less that such accurate measurements of T exist. 
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waves), and as a consequence the mean square multipoles T^f are smooth functions of I. 
One can therefore replace the sum over / by an integral, 




L/;,(2, + l)E?^ (47) 



47r 

position 

Another nice feature of large I is the relation 9^1/1 between the angular size of a feature in 
the sky (in radians) and the order / of the multipoles that dominate it. (This is, of course, 
analogous to the relation r ~ between the linear size of a feature and the wavenumbers 
that dominate its Fourier expansion.) Translating to degrees we have the following relation 
between / and the angular scale 

From Eq. ( |4.7] ), one sees that (21 + l)/$]^/47r gives the contribution to the mean square 
anisotropy per unit logarithmic interval of /, and therefore the expected temperature con- 
trast of structures in the sky with angular size 6 = 60/1 degrees. It is commonly referred to 
as 'the cmb anisotropy on scale 0\ As yet it has been observed only in the COBE observa- 
tion which we shall discuss in Section but the theoretically expected value in the CDM 
model is shown in Figure 4. Before discussing this result in detail for large angular scales, 
let us give a brief overview starting from the smallest scales. 

On scales less than a few arcminutes, the distance Eq. ( [4.1[ ) subtended at the last 
scattering surface is less than the thickness of this surface, which is about 14Mpc (Hogan, 
Kaiser &: Rees 1982). As a result, any initial anisotropy is wiped out, and it turns out that 
no significant additional anisotropy is acquired by the cmb on its journey towards us. 

On bigger scales the cmb has significant anisotropy when it starts out, which can be 
calculated in terms of the cosmological perturbations of the matter and radiation. An scales 
less than a degree or so it acquires negligible additional anisotropy on its journey towards 
usj^ From Eq. (4.1) this arcminute regime explores distance scales of order 10 to 100 Mpc, 



the same as the one explored by large scale galaxy surveys. 

As the scale is increased towards one degree, the cmb begins to acquire significant 
additional anisotropy on its journey towards us, which dominates on scales in excess of 
a degree or so. This anisotropy, called the Sachs- Wolfe effect, is caused by the density 
perturbation, and (if they exist with significant amplitude) by cosmological gravitational 
waves. Eq. ([4.1| ) still holds roughly, so we see that scales in excess of a few degrees explore 
scales in excess of 10^ Mpc. 

The expected anisotropy from these sources is shown in Figure 4. In this report we 
deal only with the large scale anisotropy, corresponding to angular scales 9 in excess of a 
degree or so, and multipoles with / < 60 We have three reasons for making this restriction. 
First, the COBE detection which is the only one reported at the time of writing lies in this 
regime. Second, the theoretical description in terms of the Sachs- Wolfe effect is extremely 
simple. Finally, this is the regime which explores very large linear scales, which cannot be 



^^We here discount the Sunyaev-Zel'dovich effect by which the cmb acquires anisotropy through photon 
up-scattering off dust on its way through a galaxy or cluster. It occurs on scales of order one arcminute and 
is easily identified for nearby clusters. The cumulative effect of distant clusters is expected to be negligible 
(Makino & Suto 1993). 
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explored using galaxy surveys. For a discussion of smaller scales the reader is referred to 
the papers of Bond and Efstathiou (1987), Linder (1988), Bond et al (1991a), Vittorio et al 
(1991) and Dodelson and Jubas (1992). 



4.2 The contribution of the density perturbation 

We now calculate the CMB anisotropy associated with the density perturbation. First 
we calculate the anisotropy of the cmb from this source when it first set out on its journey 
towards us, and then we calculate the Sachs- Wolfe effect describing the additional distortion 
which it experiences on the way. In calculating this latter effect we do not employ the usual 
description in terms of the metric perturbation. Instead, we work with the peculiar velocity 
field, imagining that the cmb passes a sequence of comoving observers. Peebles (1980) has 
given a similar treatment, though he does not explicitly work with the peculiar velocity. 

In the small region around a comoving observer, use a locally inertial coordinate system 
in which he is instantaneously at rest. Denote these coordinates by with = the 
observer's position, and consider light received by the observer which was emitted by a 
comoving source with coordinates dy* = e^dy, with a unit vector. Its wavelength is 
redshifted by an amount dA, given in terms of the recession velocity by 

d\/\ = e'dui (4.9) 

= ^'^dy^ (4-10) 

= e'e^^dy (4.11) 
dy' ^ ^ ^ 



For a perfectly homogeneous and isotropic universe, Eq. ( 2.60 ) gives Bubble's law 



dX/X = H{t)dt = da/a (4.12) 

To calculate the redshift of light from a distant comoving source, we may suppose that it is 
observed by a sequence of nearby comoving observers in its path, and integrate the above 
expression. This leads to the familiar result 

l + z^-^ = ^ (4.13) 

where the subscript denotes the observer's position and the subscript em denotes the 
position of the comoving source which emitted the radiation. 
Including the perturbation, Eq. (|2.6l| ) gives 

dX da ,■ . dvj(-K, t) , / . ., 

— = h e'e^ . ' dx (4.14) 

A a ax* 

Here, as always, the scale factor a is defined with respect to the average Hubble parameter 
on a comoving hypersurface, H = a/ a, and the comoving coordinates are related to 
inertial ones at each point by dy^ = adx^. Since we are dealing with a photon trajectory, 
dt = —dy = —adx. To first order in the perturbations the anisotropy 6T/T = —5X/X is 
given by 

^=(^) -pevWdx (4.15) 

T V T /em Jo dx' ^ ' 
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In this expression, Xem = ^-ffg"^ coordinate distance of the last scattering surface. 

The first term is the anisotropy of the cmb on a comoving hypersurface, evaluated just 
after it leaves the last scattering surface. On scales entering the horizon after decoupling 
it is easily calculated, because the initial adiabatic condition still holds. Since pr tx T^, 
AT{e)/T = ^5r = ^Sc, the last equality following from the adiabatic condition. Since the 
cmb starts out during matter domination, one has 6p = pc6c + Pr^r — P^c- Thus 



Ar(e) 
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1 



5(x) 



(4.16) 



where x = exem is the point of origin of the cmb arriving from the direction e. 

Now consider the second term, which is the Sachs-Wolfe contribution. According to 
Eq. ( |2.64| ) it may be written 
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(4.17) 
(4.18) 



To evaluate the integral, we need the x dependence of (t/a), which using a oc 

^2/3 

is given 

m=f'°'-^ = ^\C-) -'-] (4.1.) 



The other term may be expressed in terms of peculiar velocities using Eq. (HI). The final 
result is 

= i<5(x) + i [V'(x) - Vo] - e. v(x) + e. vo (4.20) 

In this expression, the subscript denotes our position. 

Of these terms, e.vg is the part of the dipole due to our peculiar velocity^ and —ipo/3 
has no angular dependence. Since we are dealing with scales far outside the horizon the 
other terms are related by Eqs. ( |2.68| ) and ( |2.69| ), and one sees that the potential term 
dominates. Dropping the whole of the dipole as in Eq. (|4.3|), the final result is therefore 



AT(e) 



1 



V'(x) 



(4.21) 



This is often called the Sachs- Wolfe effect. 

Inserting the Fourier expansion of iJj and projecting out a multipole leads to the expres- 
sion 



(4.22) 



where fik is the solid angle in k space. 

^^It may be shown (Lyth 1993b) that the total dipole contribution from the above expression is e.(vo — v) 
where the bar denotes the average over aU of space inside the surface of last scattering. This average is 
completely negligible compared with our own peculiar velocity, thus verifying the usual statement that the 
dipole measures our peculiar velocity for an adiabatic perturbation. An isocurvature perturbation can in 
contrast give a significant dipole (Turner 1991a). 



37 



This expression applies to any observer at the origin x = of the coordinates. When 
the position of this observer, and therefore that of the origin, is chosen randomly the 
coefficients Si^ have random phases, and this implies that the multipoles af^ have a Gaussian 
distribution. The variance of the distribution is (Peebles 1982b) 

roD At. 

^f = ^l ^jfm/aH) 5l{k) (4.23) 
where ji is the spherical Bessel function. With 6jj{k) oc k^^^ this becomes 



r((3-n)/2)r(; + (n-l)/2) 
2 ^ ^r((4-n)/2)r(Z + (5-n)/2) 



/(/ + 1) 



(4.24) 



The square bracket is equal to 1 for n = 1. For / ^ 1 and I ^ \n\ it can be replaced by 1, 
if 5h is evaluated on the scale k ~ IHq/2 which dominates the integral. 

The isocurvature contribution 



The Sachs- Wolfe effect Eq. ( 4.21| ) gives the anisotropy acquired by the radiation on its way 



towards us, whether the initial condition is adiabatic or isocurvature. For the isocurvature 
case, Eq. ( |2.71| ) allows one to write it as 

/Ar(e) 



V T 



= -S(x)/15 (4.25) 
sw 



where S is the intial entropy perturbation. In that case, however, one must add the 
anisotropy which the cmb started out with, and according to Eq. (2.57) it is five times 



as big. The total anisotropy is therefore six times bigger than in the adiabatic case, for a 
given density perturbation (Starobinsky & Sahni 1984; Efstathiou & Bond 1986; Kodama 
& Sasaki 1987). 

In the second half of this report we will find, from studies of the distribution and motion 
of the galaxies, that for a spectral index n = 1 the spectrum of the density perturbation 
has more or less the right normalisation to account for the observed cmb anisotropy if 
it is adiabatic. Since the isocurvature and adiabatic transfer functions have more or less 
the same shape, it follows that the cmb anisotropy would be much too big if the density 
perturbation were isocurvature. In other words the isocurvature density perturbation, if it 
exists, can be only be a small fraction of the total density perturbation. Lacking theoretical 
guidance about the size of this fraction, we ignore the isocurvature contribution from now 
on, including its contribution to the cmb anisotropy which could in principle be significant. 

4.3 The contribution of gravitational waves 



Eq. ( [4.15 ) for the anisotropy of the cmb involves the gradient diUj of the four- velocity of 



comoving worldlines. To the extent that gravitational waves are significant, this gradient 
cannot be written as the gradient of a globally defined peculiar velocity. Rather one must 
write (Peebles 1980; Lyth 1993a) 



^ = + a-\t)^-^^ + i^MSll) (4.26) 
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The second term is the one we dealt with aheady, and the third is the contribution of 
the gravitational waves. It involves that part of the metric perturbation in Eq. (^^) which 
is traceless, S^^hij = 0, and transverse, dihij = 0. This means that each Fourier component 
is of the form hij = h+efj + hxe^j, where in a coordinate system where k points along the 
z-axis the nonzero components of the polarisation tensors are defined by e+, = —^yy = 1 
and e^y = e^^ = 1. The spectrum Vg of the gravitational wave amplitude may be defined 



by summing Eq. ( |3.lD over all four components, 

= ^ (S)' + ^l^x(k)P)) (4.27) 



2lT J 

Each Fourier component satisfies the massless wave equation, which in comoving coordinates 
is 

hij + ?,Hhij + {k/afhij = (4.28) 

Well before horizon entry it has constant initial value. For scales entering the horizon after 
matter domination its subsequent evolution is 



hij{t) 



2 ,vr J3/2(a;) 



/ijj (initial) (4.29) 



2 x3/2 

where x = 2k /(aH). Well after horizon entry one has redshifting radiation. 



By substituting this expression into Eq. ( 4.26 ) one can calculate the cmb multipoles 



in terms of the initial amplitude (Fabbri & Pollock 1983; Starobinsky 1985), and hence 
calculate in terms of the spectrum Vg{k) of the initial amplitude. Each gravitational wave 
gives its dominant contribution as it enters the horizon, since its amplitude is practically 
constant before that and redshifts away afterwards. As a result the gravitational wave 
contribution to the cmb anisotropy cuts off below the angular scale ~ l'^, corresponding 
to the smallest linear scale which enters the horizon after decoupling fejg^ = 180 Mpc. The 
corresponding multipole cutoff is at / ~ k^^^Ho/l ~ 70. For / well below this cutoff the 
gravitational wave contribution is given in the case of a flat spectrum {Vg{k) independent 
of k) by (Starobinsky 1985) 

If one ignored the cutoff due to the redshift, the coefficient Ci would become 1 in the limit 
/ » 1. Starobinsky gives the values C2 = 1.118, C3 = 0.878 and C4 = 0.819. Numerical 
calculation, including the effect of the cutoff, shows that a value Q close to 1 is indeed 
achieved for I ~ 10, before the cutoff takes effect (Figure 4). 

For I » 1 the above result is good also \iVg{k) has moderate scale-dependence, provided 
that it is evaluated at the scale k ~ IHq/2 which dominates the lib. multipole. 

Within a given inflationary model one can calculate Vg{k). Defining the spectrum Ug of 
the gravitational waves by Vg oc /c"", one finds 

Ug = -2ei (4.31) 

where ei is the small positive parameter mentioned earlier and defined in Section |5|. 
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Setting the coefficient C/ of Eq. (|4.30| ) equal to 1, tlie ratio of the gravitational and 
density contributions is given in terms of this same parameter by 



T,f (grav 
(density) 



R = . ' , = 12.4ei (4.32) 



In most of the usually considered models of inflation, ei is very small, corresponding to 
a very small gravitational wave contribution with a very flat spectrum. The important 
exceptions are power-law and extended inflation, where ei need not be very small, so that 
the spectral index can be much less than one (more power on large scales) and the ratio 
R can be significantly bigger than 1. It is impossible to have a positive spectral index 
for the gravitational waves, without going outside the almost universally accepted slow-roll 
paradigm. 

As we noted earlier the gravitational wave amplitude is strongly reduced by redshifting 
on scales < 100 Mpc, so that it gives a negligible contribution to the cmb anisotropy on 
arcminute scales. Nevertheless, estimates with a flat spectrum Ug = show that the waves 
could give a significant contribution to the planned LIGO detector (hinder 1988; Sahni 
1990; White 1992; Souradeep &: Sahni 1992). Since the spectrum is tilted, if at all, towards 
more large scale power, these estimates are upper limits. They ensure by many orders of 
magnitude that the gravitational waves do not contribute significantly to the locally defined 
peculiar velocity field, which we observe within 100 Mpc or so of our position. 



5 Inflation 

It is widely supposed that the very early universe experienced an era of inflation. By 
'inflation' one means that the scale factor has positive acceleration, d > 0, corresponding 
to repulsive gravity and 3p < —p. During inflation aH = d is increasing, so that comoving 
scales are leaving the horizon (Hubble distance) rather than entering it, and it is supposed 
that at the beginning of inflation the observable universe was well within the horizon. 

The inflationary hypothesis is attractive because it holds out the possibility of calculating 
cosmological quantities, given the Lagrangian describing the fundamental interactions. The 
Standard Model, describing the interactions up to energies of order 1 TeV, is not viable in 
this context because it does not permit inflation, but this should not be regarded as a serious 
setback because it is universally agreed that the Standard Model will require modification at 
higher energy scales, for reasons that have nothing to do with cosmology. The nature of the 
required extension is not yet known, though it is conceivable that it could become known in 
the forseeable future (Weinberg 1993). But even without a specific model of the interactions 
(ie., a specific Lagrangian), the inflationary hypothesis can still offer guidance about what 
to expect in cosmology. More dramatically, one can turn around the theory-to-observation 
sequence, to rule out otherwise reasonable models. 

In this section we explain how these things work, with particular reference to the density 
perturbations and to cosmological gravitational waves. Both of these could originate during 
inflation as quantum fluctuations, which become classical as they leave the horizon and 
remain so on re-entry. The original quantum fluctuations are of exactly the same type as 
those of the electromagnetic field, which give rise to the experimentally observed Casimir 
effect. 
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5.1 The basic features of inflation 

The implications of an inflationary era were worked out by many people, following the 
seminal paper by Guth (1981), and they have been described in several texts (Kolb & 
Turner 1990; Linde 1990). One of the most dramatic and simple effects is that there is 
no fine-tuning of the initial value of the density parameter Q = Sirp/Sm'piH'^ . From the 
Friedmann equation, is given by 

Its present value fio is certainly within an order of magnitude of 1, and in the absence of an 
inflationary era becomes ever smaller as one goes back in time, implying an initial fine 
tuning. In contrast, if there is an inflationary era beginning when the observable universe is 
within the horizon, Eq. ( ^.1| ) implies that JIq will be of order 1, provided only that the same 
is true of at the beginning of inflation. A value of extremely close to 1 is the most 
natural, though it is not mandatory (Lyth & Stewart 1990a; Ellis, Lyth &; Mijic 1991). |^ 

Another effect of inflation is that it can eliminate particles and topological defects which 
would otherwise be present. Anything produced before inflation is diluted away, and after 
inflation there is a maximum temperature (the 'reheat' temperature) which is not high 
enough to produce all the particles and defects that might otherwise be present. As we 
shall remark later, this mechanism can remove desirable, as well as undesirable, objects. 

The most dramatic effect of inflation is that it may offer a way of understanding the 
homogeneity and isotropy of the universe, or at any rate of significant regions of it. We have 
nothing to say about this complex issue in its full generality, but a more modest version 
of it is our central concern. In this version, one begins the discussion at some early stage 
of inflation, when the universe is supposed already to be approximately homogeneous and 
isotropic (Section 5.1). One then argues that in that case, scales far inside the horizon 
must be absolutely homogeneous and isotropic, except for the effect of vacuum fluctuations 
in the fields. Finally, one shows that after they leave the horizon, such fluctuations can 
become the classical perturbations that one deals with in cosmological perturbation theory. 
This possibility was first pointed out for gravitational waves by Starobinsky (1979) and for 
density perturbations by several people (Guth &; Pi 1982; Hawking 1982; Starobinsky 1982) 
As we shall go to some trouble to demonstrate, the vacuum fluctuations can be evaluated 
unambiguously once an inflationary model is specified. 

To have a density perturbation of the right magnitude, the energy density when the 
observable universe leaves the horizon must typically correspond to an energy scale p^^^ ~ 
10^^ GeV, not far below the Planck scale. One is thus making a big extrapolation above 
the scales < 10^ GeV at which the standard model has been tested (and an even 

bigger one above the scale T ^ 10^^ GeV, corresponding to the epoch a few Hubble times 
before the smallest cosmologically interesting scales enter the horizon, which are all that 
we need to consider in the rest of the present report). Mercifully, a density perturbation 
leaving the horizon during inflation can be followed to the epoch of horizon entry without 

^''An argument has been given for fio very close to 1 on the basis of effects on the cmb anisotropy from 
regions far outside the observable universe (Turner 1991a), but it is not vahd as it stands because it ignores 
spatial curvature. 
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(5.1) 



any knowledge of the intervening era, as we saw in Section |2[ But one would be very lucky 
to find a well motivated model of inflation without understanding that era. At present such 
an understanding requires that one believes in some model of the GUT type which has 
elementary Higgs fields, as opposed to one of the technicolour type where the Higgs fields 
are composite. 



Scalar field inflation 

Two mechanisms for inflation have been proposed. The simplest one (Guth 1981) invokes 
a scalar field, termed the inflaton field. An alternative (Starobinsky 1980) is to invoke a 
modification of Einstein gravity, and combinations of the two mechanisms have also been 
proposed. During inflation however, the proposed modifications of gravity can be abolished 
by redefining the spacetime metric tensor, so that one recovers the scalar field case. We 
focus on it for the moment, but modified gravity models will be included later in our survey 
of specific models. 

In comoving coordinates a homogeneous scalar field (p with minimal coupling to gravity 
has the equation of motion 

4) + 3H<p + V'{(t>) =0 (5.2) 
Its energy density and pressure are 

p = V + ^<p'' (5.3) 
p = -y + ^02 (5,4) 

If such a field dominates p and p, the inflationary condition 3p < p is achieved provided 
that the field rolls sufficiently slowly, 

02 < y (5.5) 

Practically all of the usually considered models of inflation satisfy three conditions. 
First, the motion of the field is overdamped, so that the 'force' V' balances the 'friction 
term' 3i?(/>, 

(/>~-^y' (5.6) 

Second, 

which means that the inflationary requirement (p'^ < V is well satisfied and 

H'^l^V (5.8) 

These two conditions imply that H is slowly varying, and that the scale factor increases 
more or less exponentially, 

a oc e^* (5.9) 
The third condition that is usually satisfied is 

\v\ < 1 (5.10) 
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where 



It can be 'derived' from the other two by differentiating the approximation Eq. ( ^.61) for (j) 
and noting that consistency with the exact expression Eq. ( |5.2| ) requires (j) <^V' is satisfied. 
However there is no logical necessity for the derivative of an approximation to be itself a valid 
approximation, so this third condition is logically independent of the others. Conditions 
involving higher derivatives of V could be 'derived' by further differentiation, with the same 
caveat, but the two that we have given, involving only the first and second derivatives, are 
the ones needed to obtain the usual predictions about inflationary perturbations. The term 
'slow-roll inflation' is generally taken to denote a model in which they are satisfied and we 
are adopting that nomenclature here. Practically all of the usually considered models of 
inflation satisfy the slow-roll conditions more or less well, and we adopt them in this section, 
showing how to handle the general situation in Appendix B. 

It should be noted that the first slow-roll condition is on a quite different footing from 
the other two, being a statement about the solution of the field equation as opposed to a 
statement about the potential that defines this equation. What we are saying is that in the 
usually considered models one can show that the first condition is an attractor solution, in 
a regime typically characterised by the other two conditions, and that moreover reasonable 
initial conditions on cj) will ensure that this solution is achieved well before the observable 
universe leaves the horizon. 



Begining and ending inflation 

Given a parametrisation of the potential, one can calculate the spectrum of the density 
perturbation and of the gravitational waves, and see if they are consistent with observation 
in a theoretically acceptable regime of parameter space. A parametrisation of the potential, 
suggested more or less loosely by some theoretical consideration, is referred to as a model 
(Linde 1990; Olive 1990; Kolb 1991). 

Models of infiation may usefully be classified according to the way in which they enter 
the inflationary era, and the way in which they exit that era. Two modes of entry have been 
discussed, which one may term thermal and chaotic. In the thermal mode, one supposes 
that infiation is preceded by an epoch of radiation domination, with at least some of the 
components in thermal equilibrium. Specific attempts to realise the thermal case, which 
will be discussed in more detail below, are 'old infiation', 'new infiation' and 'extended 
infiation'. In the thermal case, one does not usually discuss what happened before thermal 
equilibrium set in. In the chaotic mode it is supposed that the universe, or at least our 
part of it, contains only smoothly varying fields from the moment of its creation until the 
end of infiation. Subsequently, most of the fields quickly adjust themselves to minimise the 
eff'ective potential, leaving behind the slowly rolling infiaton field and perhaps some other 
fields. 

Two modes of exit are commonly discussed, which one may term first-order and os- 
cillatory. In the case of first-order exit, the universe is supposed to be sitting in a false 
vacuum during inflation. Exit from inflation occurs via a first-order phase transition, in 
which bubbles of the true vacuum form by quantum tunnelling. The original (scalar field) 
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inflationary model, usually called 'old inflation', was of this kind, and was quickly shown 
not to be viable because the bubbles never coalesce. More recently, flrst-order inflation 
models termed 'extended inflation' have been formulated, which rely on a modification of 
gravity to coalesce the bubbles. As we shall show, the COBE data when combined with 
large scale galaxy surveys puts these models in jeopardy too (Liddle & Lyth 1992, 1993). 
Different first-order infiation models discussed below, employing two scalar fields, may still 
be viable. 

In the case of an infiation is supposed to end when the inflaton field moves sufficiently 
close to its minimum that it starts to oscillate, without the intervention of any phase 
transition. For the commonly considered potentials, this epoch is marked by the failure of 
one of the conditions Eqs. ( ^.71) and ( |5.10| ), so denoting it by the subscript 'end' we have 



The oscillations correspond to a gas of 'inflaton' particles with negligible pressure, so that 
there ensues a matter dominated era which ends only when the inflatons decay into radia- 
tion. Of the thermal-entry models mentioned earlier, 'new inflation' invokes an oscillatory 
exit. Chaotic inflation is normally also thought of as invoking an oscillatory exit, though 
this is not mandatory. 

In addition to these two commonly considered cases, there exists the possibility that 
inflation occurs in a false vacuum, which is however exited through a second order phase 
transition as opposed to a flrst order one. That is, the local minimum of the effective 
potential which represents the false vacuum turns smoothly into a local maximum, after 
which the fields quickly adjust themselves to the true vacuum values. By definition, bubble 
production is negligible in a second order phase transition, but as with a first order transition 
there may be production of topological defects such as cosmic strings or domain walls. A 
simple two field model which seems to give rise to a second order exit has been given by 
Linde (1991b); it may be termed two-scale inflation because the fields have very different 
masses m and M, related roughly be mmpi = M^. We will describe the two-scale model in 
Section |5|.4, along with a selection of better known ones. 

In any model, the universe (or at least the constituents dominating the energy density) 
must achieve thermal equilibrium at some point after inflation. Because the first models 
of inflation invoked a thermal entry, this event is usually called 'reheating', even in the 
chaotic scenario where no previous thermalisation has occurred. The corresponding 'reheat 
temperature' is the biggest temperature ever achieved after inflation and plays a crucial role 
in cosmology. Unfortunately, no reliable estimate of it is known at present. Reheating will 
occur promptly if inflation ends in a phase transition (certainly if it is a flrst-order one), but 
it may be long delayed if inflation ends with oscillations of the inflaton fleld. The reason is 
that the couplings of the inflaton fleld to other flelds is typically quite weak. A discussion 
of reheating, with references to earlier work, has recently been given by Kofman, Linde and 
Starobinsky (1993). 

The epoch of horizon exit 

A first requirement for an infiationary model is that the observable universe should be 
within the horizon at the beginning of infiation. We therefore need to know the epoch 



max{eend, \Vend\} - 1 
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when a given scale leaves the horizon, something which is necessary also to calculate the 
inflationary perturbations. Denoting it by a star, this epoch is given by 



k 



(5.13) 



The epoch of horizon exit is therefore related to the present magnitude of the scale in 
Hubble units by 

(5.14) 



k 



Let us denote the end of inflation a subscript 'end' and the epoch of reheating by 'reh', 
assuming a matter domination during the era (if any) between these two epochs. Let us 
also assume that after 'reh' there is radiation domination, until the epoch 'eq' at which 
the dark matter density becomes equal to that of the radiation. Throughout the history 
of the universe the radiation energy density is proportional to a~^, that of the matter is 
proportional a~^, and the total is proportional to H^. It follows that 
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(5.16) 



where pQr = (aeq/ao)yOo is the present radiation energy density and A^* is the number of 
Hubble times between horizon exit and the end of inflation, 



= ln(a^/aend) 

*cnd 

Hdt 



(5.17) 
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It follows that 
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This equation relates the three energy scales V* 
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are related by another expression for A'^* 
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which follows from Eqs. ( 5.18|) , (|5.6|) , (5.7) and 
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The biggest scale that can be explored is roughly the present Hubble distance, ao/k = 
Hq^ = 6000 Mpc. We shall refer to the epoch of horizon exit for this scale as 'the epoch when 
the observable universe leaves the horizon', and denote it by the subscript 1. As displayed in 
Eq. ( ^.461) , the COBE observations require that y/^^ < 10^*^ GeV, with the equality holding 



in most models of inflation. Also, Eq. ( 5.2C| ) gives in most models of inflation V^^'^ ~ 



.1/4 
end ■ 

If reheating is prompt, we learn that in most models the observable universe leaves the 
horizon about 62 e-folds before the end of inflation. If reheating is long delayed A^i could be 
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considerably reduced, being equal to 32 for the most extreme possibility of Vj.J^ ~ 1000 GeV 
(corresponding to reheating just before the electroweak transition which is presumably the 
last opportunity for baryogenesis) . For most purposes, however, one needs only the order 
of magnitude of A'^i. 

The smallest scale on which the primeval perturbation can be probed at present is 
around 1 Mpc, and one sees that this scale leaves the horizon about 9 Hubble times after 
the observable universe. Discounting the case of very late reheating, we conclude that in 
the usually considered inflationary models, scales of cosmological interest leave the horizon 
50 to 60 Hubble times before the end of inflation. 

5.2 Calculating the density perturbation 

Soon after inflation was proposed, it was realised that the vacuum fluctuation of the inflaton 
field would inevitably give rise to an adiabatic density perturbation. Its spectrum, specified 
say by the quantity dnik) that we introduced earlier, can be calculated in terms of the 
inflationary potential. The first calculations (Guth &: Pi 1982; Hawking 1982; Starobinsky 
1982; Bardeen, Steinhardt & Turner 1983) quoted only order of magnitude results, largely 
because they did not specify the normalisation of the spectrum that the result was supposed 
to apply to, but partly also because the calculations themselves were not very precise. The 
precise result, valid subject to the three slow-roll conditions listed above, was first given by 



As always, a star denotes the epoch of horizon exit aH = k. 

We now derive this result, using essentially the original method with a simplification due 
to Sasaki (1986). Here in the text we ignore the back-reaction of the metric perturbation 
on the inflaton field, explaining in Appendix B why it is negligible within the slow-roll 
paradigm. We also explain there how departures from slow-roll, including the back-reaction, 
can easily be evaluated if desired (Stewart & Lyth 1993). A somewhat related formalism 
was developed by Mukhanov (1985, 1989), and yet another by Salopek, Bond and Bardeen 
(1989), but neither of them is as simple as the one advocated here, especially in the context 
of the slow-roll paradigm. 

To define the perturbation 5(f) of the inflaton field, one has to choose a slicing of space- 
time into spacelike hyper surf aces, just as was the case for the perturbations 5p, 5p and 
SH that we studied earlier. Since the inflaton field is supposed to dominate the energy- 
momentum tensor, the momentum density vanishes if its spatial gradients vanish (Bardeen, 
Steinhardt & Turner 1983). In other words, 5(p vanishes if the hypersurfaces are chosen to 
be comoving! 

As we explain in Appendix B, the comoving choice of hypersurfaces becomes singular in 
the limit e — > of exponential inflation. If instead we make a choice which is non-singular, 
it becomes very easy to handle the inflaton field perturbation. Its field equation, given 
the slow-roll conditions, is the one that would naively be obtained by perturbing Eq. ( ^.21) , 
ignoring the perturbation of the metric. 



Lyth (1985), 
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Until a few Hubble times after horizon exit Eq. ( |5.10| ) ensures that V" can be dropped, 
so that the equation becomes 



iScl)ky+3H{5M-+i^-j 50k = (5.23) 

Well before horizon entry is a massless field living in flat space-time, since its wavenum- 
ber k/a is much bigger than the Hubble parameter H. It can be quantised in the usual 
way so that its quantum state is labelled by the number of inflaton particles present with 
momentum k. 

The fundamental assumption leading to Eq. ( |5.21| ) is that on cosmologically relevant 
scales the inflaton field is in the vacuum state, corresponding to no inflaton particles. Let 
us pause to ask about the status of this assumption .p^ The idea behind it is that at any 
epoch during inflation there is a cutoff /cmax in the energy spectrum of inflaton particles. 
Such a cutoff is obviously present in order to avoid infinite energy density, but the point 
at issue is its magnitude, in Hubble units. If inflation is preceded by an epoch of thermal 
equilibrium, the equilibrium occupation number (Kolb & Turner 1990) 2(exp(A;/aT) it 1)~^ 
implies /cmax/o ~ T. In that case inflation begins when the radiation energy ~ T'* falls 
below the inflaton fleld energy ~ V, so the cutoff is independent of the epoch and is given 
by (femax/a)^ ~ ^begin- But since V decreases during inflation, and a increases rapidly, this 
implies that once inflation is underway 

km...V (5.24) 



Irrespective of what happens before inflation, this condition is in any case necessary for 
consistency, just to ensure that the radiation energy density due to the inflaton particles is 
much less than the potential V which is supposed to be dominant. Indeed, since the density 
of states in k space is (27r)~^, one inflaton per state up to A^max implies a radiation energy 
density given by 

1 r^^^^(k\Uk 1 //c^axV 



,\ 1/2 

(5.26) 



Apart from the numerical factor this leads to the desired condition, Eq. (5.24). 

The condition is enough to justify the vacuum assumption. To see this, rewrite it in the 
form 

'^max\ / mp^ y/^ 
aH J \ H ) 

We shall see that typically the right hand side is of order 10^, so assume this value for 
deflniteness. Now suppose that inflation begins some number A'before of Hubble times before 
the observable universe leaves the horizon. Since the scale factor a grows like e^*, with H 
is slowly varying, Eq. ( ^.26|) applied at the beginning of inflation leads to 

^ < 103e-^^^f°- (5.27) 
ki 

^^The following arguments, which we have not seen written down before, were developed in collaboration 
with Ewan Stewart. 
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The vacuum assumption is therefore ensured if inflation begins at least 7 or so Hubble times 
before the observable universe leaves the horizon. Clearly, the result is not very sensitive 
to the assumed value of mpi/H. 



Now let us see how the spectrum Eq. (5.21) is implied by the vacuum assumption. The 
first step is to calculate the vacuum expectation value of If^t/'kP- Well before horizon entry 
we need only the flat space-time version of quantum field theory, rewritten in terms of the 
comoving coordinate and the comoving wave vector k. They are related to the Cartesian 
coordinates and physical wave vectors by the scale factor a{t), whose time dependence is 
negligible on the timescale a/k. Working in the Heisenberg representation, (pi^ is associated 
with an operator 

k{t) = Wk{t)ai, + wUt)al^ (5.28) 
The annihilation operator at satisfies the commutation relation 

["k 1,4 2] = '^ki,k2 (5.29) 
and the field satisfies the commutation relation, 

ikyut), ^^ky2,t)] = i6^{ayi - aya) (5.30) 

= iaHiyi - ya) (5.31) 

As a result the functions Wk are given by 

Wk{t) =a-3/2(2A:/a)-i/2g-i(x+(fet/a)) ^5^32) 

The phase factor x is arbitrary, and like a it must have negligible variation on the timescale 
a/k. The vacuum state is the one annihilated by a^, so the vacuum expectation value of 
the field perturbation is 

(|'5</.k|') = \wk\^ (5.33) 

To extend these results to the epoch of horizon exit and beyond, we have to accept the 
validity of free field theory in curved space-time. It seems to have first been formulated 
for a general Robertson- Walker metric by Parker (1969), and an extended account of the 
subject in this and other space-times is given in the text of Birrell and Davies (1982). 
Hardly any of this general theory is in fact necessary for the present application; all we 
need to assume is that there is a Heisenberg picture, in which operators satisfy the classical 
equations of motion and state vectors are time independent. Then (/)k continues to satisfy 
the field equation Eq. ( |5.22| ), and Eqs. (5.28), (5.29) and ( ^.33 ) still hold, where Wk is the 



solution of the field equation reducing to Eq. ( 5.32 ) well before horizon entry. As one easily 
checks, the required solution is 

Mt) = v:^(^ + :^)e^'/^^ (5.34) 



(2A;3)i/2 V aH 

A few Hubble times after horizon exit, the vacuum expectation value is therefore 

= ^ (5.35) 
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A measurement of the 0k's will yield random phases, and a distribution of moduli whose 
dispersion is given by Eq. ( ^.35 ). Accordingly the spectrum of the inflaton field, defined by 



Eq. (3.1), is given a few Hubble times after horizon exit by (Vilenkin & Ford 1982; Linde 



1982; Starobinsky 1982) 

V]!\k) = ^ (5.36) 

This equation does not generally hold more than a few Hubble times after horizon exit, for 
two separate reasons. First, Eq. ( |5.34| ) fails to be a solution of the massless field equation 
Eq. ( p.23| ) unless H is practically constant, as opposed to just slowly varying. Indeed, after 
horizon exit the third term of Eq. ( 5.23| ) can be dropped, and its solution remains constant 



whether or not H varies. Second, the massless equation will not in general be correct, 
because the term V" in Eq. ( ^.22| ) will become significant. In order to derive our desired 
result however, it will be enough to know that Eq. ( ^.36 ) holds a few Hubble times after 



horizon exit. Since H is slowly varying on the Hubble timescale it can be written 

V'j\k) = ^ (5.37) 

where the star denotes the epoch of horizon exit, and V^p is evaluated a few Hubble times 
after horizon exit. 



Massless field inhomogeneity 



Before continuing the derivation of the spectrum of the density perturbation, we briefly 
digress to consider the case of a massless scalar field ip. Provided that the field has minimal 
coupling to gravity it satisfies Eq. ( ^.23 ), so one arrives at Eq. ( 5.37| ) with (p replaced by 'ip, 



2^ 



(5.38) 



In contrast with Eq. ( ^.37 ) though, this equation remains valid arbitrarily long after horizon 
exit, because as we noted already the massless field equation ensures that ipk remains 
constant. Defining a spectral index by oc k"^^ or equivalently 



din v.. 



dink 



(5.39) 



one concludes from Eqs. (IST 



(|5]|), ( p3D and <^ that (Liddle & Lyth 1992) 

= -2ei (5.40) 

To the extent that the spectral index is well defined e does not vary significantly while inter- 
esting scales leave the horizon, but for definiteness we may evaluate it when the observable 
universe leaves the horizon and this is indicated by the subscript 1. 

If the inhomogeneity of ip is converted later into an isocurvature density perturbation 
S as discussed in Section |^.7, the spectrum Vs has the same shape, so Eq. ( 5.40| ) gives its 
spectral index. 
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The curvature perturbation 

Returning to the derivation of Eq. ( |5.21| ), we have so far deah with the inflaton field pertur- 
bation 5(f), defined on hypersurfaces which are non-singular in the limit e — > 0. It is shown 
in Appendix A that the curvature perturbation IZ of comoving hypersurfaces is given in 
terms of this quantity by 

7^ = ^54) (5.41) 

A relation equivalent to this one was derived by Lyth (1985), but its simple form in terms 
of the IZ variable was first appreciated by Sasaki (1986), who pointed out that it followed 
directly from a formula given by Bardeen (1980). A simple derivation of it is given in 
Appendix A. 

The spectrum of IZ is therefore given by 

= "^Vf (5.42) 

In Section §.4 we learned that 7Z is constant after horizon exit, so Vn remains constant 
even though H/cj) and V^j, might vary separately. It follows that as long as the scale is far 
outside the horizon, 

,1/2 _ 

where the star denotes the epoch of horizon exit. Using Eq. (|5.6| ), this leads to Eq. (5.21), 
concluding our proof of that equation. 

The spectral index 

The spectral index n of the density perturbation is defined by 5"^ oc fc""^, or equivalently 

The power-law dependence on k is not supposed to be exact, so that the latter expression 
should be taken as the proper definition of the spectral index, being useful if n does not 
vary much over a cosmologically interesting range of scales. Differentiating Eq. ( ^.21 ) with 



V^' = -M- (5.43) 



the aid of Eqs. (^, (^, ( pT^ ) and (^ leads to the resuh (Liddle & Lyth 1992) 

n = 1 + 2r?i - 6ei (5.45) 

To the extent that n is scale-independent the right hand side can be evaluated at any epoch, 
which for definiteness we have taken to be the epoch when the observable universe leaves 
the horizon. 

In order to have significant deviations from the fiat n = 1 spectrum one must come 
close to violating at least one of the slow-roll conditions e <C 1 and \r]\ <^ 1. Nevertheless it 
turns out that three of the half dozen or so currently favoured models can naturally give a 
value of n significantly below 1. It is more difficult to have n significantly above 1 because 
e is positive by definition; none of the usually considered models has n > 1 but we shall see 
that small positive values are possible in two-scale models. 
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The inflationary energy scale 

As we see shall see later, the cmb anisotropy observed by COBE determines the value of 
Sh on the scale corresponding roughly to the size of the observable universe, if the density 
perturbation is adiabatic and dominates the anisotropy. The value is = 1.7 x 10^^ with 
a one-sigma uncertainty of about 20%. From Eq. ( ^.21 ) this determines the potential and 



the Hubble parameter in terms of ei, at the epoch when the observable universe leaves the 
horizon 

f//^ = 6.2e}/^ X 10^^ GeV (5.46) 
Hi = 9.1e}/^ X 10^^ GeV (5.47) 

If there are other contributions to the cmb anisotropy these estimates are reduced. One 
therefore has the important conclusion (Lyth 1984) that the inflaton potential when the 
observable universe leaves the horizon satisfies V^^^ < 10^^ GeV. This is also an upper 
bound on the potential Vend at the end of inflation, and in fact with very weak assumptions 
(Lyth 1990) one can say more precisely that V^^^ < 1.6 x 10^^ GeV. 

The classical nature of the predicted perturbation 

To end this long section, we want to explain in what sense the predicted density perturbation 
is a classical quantity. Similar considerations apply to the predicted fluctuation of a massless 
field, and to the predicted gravitational waves to be discussed in a moment. 

In arriving at Eq. ( |5.37 ) for the spectrum well after horizon exit, we supposed that a 



measurement of the field had been made, so that the squared Fourier components |(/>kP took 
on values drawn from the probability distribution predicted by the theory. The expected 
value (It/ikP) is given by the theory, and this is all that is required to evaluate the spectrum. 
Nevertheless, what one observes (through the density perturbation which is linearly related 
to it) is the field itself, as a function of spacetime position. The question is, whether the 
quantum uncertainty is small enough to permit this field to have a sharply defined value 
over an extended period of time. The answer turns out to be yes, but for this positive 
answer it is crucial that the field is measured only well after horizon exit. 

Our fundamental hypothesis is that each perturbation starts out, well before horizon 
entry, as a vacuum fiuctuation of a massless field in fiat spacetime. The state of the system 
is the vacuum state, in which the expectation value of the field vanishes. A measurement 
of the field made well before horizon entry would produce a new state in which it has a 
nonzero expectation value, but of course we have to suppose that no such measurement is 
made or our fundamental hypothesis would be violated. The probability of each possible 
outcome of the measurement would be given by quantum theory; to be precise, the real and 
the imaginary part of each Fourier component is uncorrelated, with a Gaussian probability 



distribution whose variance is given by (half of) Eq. ( 5.33 ). 

The fiat spacetime vacuum fiuctuation is an essentially quantum object because it is not 
possible to put the field into a condition where it remains sharply defined over an extended 
period of time, unless the field strength is far out on the tail of the probability distribution. 
This is clear from the fact that the real and imaginary parts of each Fourier component 
have the same dynamics as a quantum oscillator, the vacuum state being the ground state 
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of the oscillator; it is well known that for such an oscillator it is not possible to make a 
'wave packet' which does not spread significantly, unless its displacement from the origin is 
far in excess of the root mean square value for the ground state. 

The situation is quite different well after horizon exit. Then, one can check that it 
is possible to give the field a sharply defined value, whose strength as measured by the 
modulus squared of the Fourier coefficients is around the quantum expectation value (Lyth 
1985; Guth & Pi 1985; Grishchuk 1993; Albrecht et al 1993). The underlying reason for the 
classical behaviour seems to be the fact that Wk becomes real after horizon exit (Starobinsky 
1982, 1986) so that (/>k has only trivial time dependence 

4>k^Wk{t)[ak + al^^] (5.48) 

This quantum-to-classical behaviour is a tremendous success for the theory, which does 
not seem to have received adequate publicity. If it had failed, the prediction for spectrum 
of the density perturbation would have had nothing to do with reality. 

Of course there remains the usual interpretation problem, arising whenever one talks 
about measurement in quantum theory. Even after horizon exit one is not allowed to say 
that the (/>k's have particular values before they are measured, because in principle one is 
supposed to be able to measure observables whose operators do not commute with the (j)kS. 
A similar situation arises when one observes an electron created in, say, beta decay. Even if 
its observed position is a metre away from the source, one is not allowed to say that it arrived 
at that position on the classical trajectory. The reason is that one might in principle have 
chosen to measure, say, its angular momentum, which involves the full outgoing spherical 
wave function. The problem with the inflaton field is, however, more severe because as we 
shall see it is essentially equal to the density perturbation. Whatever one might think about 
an electron, one is certainly unwilling to say that the density perturbation of the universe is 
not there until it is measured. The problem has been discussed by many authors using the 
concept of decoherence (Sakagami 1988; Halliwell 1989; Padmanabhan 1989), but it seems 
far from a satisfactory solution. 

5.3 Primordial gravitational waves? 

We stated earlier that the gravitational wave amplitude hij satisfies the massless wave 
equation. To derive this equation, one can evaluate the contribution of hij to the Einstein 
Lagrangian 

£ = -^R (5.49) 

where R is the curvature scalar of spacetime. As shown in for example the review of 
Mukhanov, Feldman & Brandenberger (1992) one finds that up to a total derivative each 
of the four non-zero Fourier components discussed in Section ^.3has the same Lagrangian 
as a Fourier component of a massless scalar field ip = (mp;/167r)^/^/ijj. 

This result not only proves the field equation, but also allows us to quantise the gravita- 
tional wave amplitude. Each non-zero Fourier component has the same vacuum fluctuation 
as a massless field, with the above conversion factor. Thus the spectrum of the gravitational 
waves is 
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Putting this expression into Eq. ( 4.3C| ) and dividing it by Eq. ( |4.24| ) gives the ratio Eq. (4.32) 
of the gravitational wave and density contributions to the cmb anisotropy. 

Just as in Eq. ( |5.4[1| ), the spectral index of the gravitational waves is (Liddle Sz Lyth 
1992) 

Hg = -2ei (5.51) 



5.4 A survey of inflationary models 

We now apply the formalism that has been developed to some popular models of inflation, 
as well as to the 'two-scale' model mentioned earlier. The survey is by no means exhaustive, 
and does not consider 'designer' inflation (Kofman & Linde 1987; Kofman & Pogosyan 1988; 
Salopek, Bond & Bardeen 1989; Linde 1992), where one tunes the shape of the potential, 
typically with more than one inflaton field, to obtain more or less any desired spectrum. 
Somewhat similar surveys have been given by Salopek (1992a, 1992b), Davis et al (1992b), 
Liddle & Lyth (1992) and Turner (1993). The outcome of our discussion is summarised in 
Table 2. 



Polynomial Chaotic inflation 

In chaotic inflation (Linde 1983, 1990) the inflaton potential is supposed to emerge, at the 
Planck scale, with a value V ~ "m-pi- If the potential is V oc (p'^ inflation can begin at 
the Planck scale, with an oscillatory exit at the epoch determined by Eq. ( ^.12 ), (f)end — 



a(m|,;/87r)i/2. From Eq. ( ^^201 ) with Ni = 60, (pj ~ 120am|,,/87r. This leads to 1 - n = 
(2 + a)/120 and R = .05a, leading to the following relation between the gravitational wave 
contribution and the spectral index 

i? = 6(l-n)-0.1 (5.52) 

For a = 2, 4, 6 and 10 one has 1 — n = .033, .05, .067 and .1, and R = .10, .20, .30 and 
.50. Shafi (1993) has shown that a potential of this kind follows from a class of superstring- 
inspired gauge theories, with an index a varying between 6 and 10. 



i?^ inflation 

Instead of using a scalar field to drive inflation one can modify gravity with the addition of 
an R^ term, 

^ = -T7P[^ + R'/m^)] + ^nongr (5.53) 

167r 

The first term is the Lagrangian for Einstein gravity, and the third term is the Lagrangian 
of the non-gravitational sector (standard model, GUTS or whatever), and the second term 
modifies Einstein gravity. A more complicated version of this Lagrangian was the first model 
of inflation (Starobinsky 1980). The simple version given here has recently been justified 
in the context of supergravity (Cardoso & Ovrut 1993), and because of its simplicity it has 
been analysed by several authors (Starobinsky 1983; Barrow & Ottewill 1983; Whitt 1984; 
Duruisseau & Kerner 1986; Mijic, Morris & Suen 1986; Starobinsky &; Schmidt 1987; Suen 
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& Anderson 1987; Kofman, Mukhanov & Pogosyan 1987). The term can be eliminated 
by a conformal transformation of the metric (Whitt 1984) to give a new Lagrangian 

C = -^i? + ]-~g^'''d^^d,(t> - V{4>) + £nongr (5.54) 
iovr z 

where 

,2 A^2 



327r 




The new Lagrangian has standard gravity, and a scalar field (j) with canonical kinetic 
terms which did not appear at all in the original Lagrangian, as well as the original non- 
gravitational sector which now appears in modified form. In the regime (j) ^^^PU the po- 
tential V^(0) satisfies the slow-roll conditions so that inflation occurs. The non-gravitational 
sector is irrelevant during inflation, and if the B? term quickly becomes negligible afterwards 
we have a model in which inflation is driven by a scalar field, without any modification of 
gravity. What has been achieved is to motivate the otherwise bizarre form of the inflationary 
potential. 

In the regime </> > mp/ where they are small, the parameters appearing in the slow-roll 
conditions are 



r, = --exp(^-^-— -<Aj (5.56) 
e = -^rf (5.57) 



There is an oscillatory exit to inflation at the epoch given by Eq. ( 5.12| ), (/ignd ~ ("^p;/8vr)^/^. 
From Eq. (|]2§) with Ni = 60, 

~ 5^^ (5.58) 

leading to r/i ~ —.02 and ei ~ 10~^. Thus the spectral index is n = 0.96, but the gravita- 
tional wave contribution is completely negligible. 

Salopek, Bond and Bardeen (1989) have considered different modifications of gravity, 
termed 'induced gravity' and 'variable Planck mass' models, which after a conformal trans- 
formation lead to inflation with the same potential as with B? inflation, in the relevant 
regime of parameter space. These models therefore give the same spectral index, and again 
negligible gravitational waves. 

Power-law inflation 

Power-law inflation a oc (Abbott & Wise 1984b; Lucchin & Matarrese 1985; Barrow 1987; 
Liddle 1989) corresponds to an exponential potential 



With the required value p > \ such a potential does not have any known theoretical mo- 
tivation, except in the context of extended inflation which we consider in a moment. In 
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contrast with the previous two models, the featureless exponential potential cannot provide 
an oscillatory end to inflation, but the calculation of the spectra of the density perturbation 
and the gravitational waves does not depend on what kind of exit occurs. The slow-roll 
conditions are satisfied for all (/> if p ^ 1. The slow-roll parameters are e = = 1/p 
independent of time. The spectral indices are given by — = 1 — n = 2/p, and R = 12/p 
leading to the relation 

R = 6{1- n) (5.60) 
The chaotic result Eq. ( 5.52| ) becomes the same as this result for 1 — n > .1. 



Power-law inflation is of particular interest to researchers in inflation because exact 
analytic solutions exist both for the dynamics of inflation and for the density perturbations 
generated. As a result one can check that the slow-roll approximation remains reasonable, 
even when the condition p S> 1 is not very well satisfied. The exact gravitational wave 
spectrum was calculated by Fabbri, Lucchin and Matarrese (1986), and the exact density 
spectrum by Lyth &; Stewart (1992a). The exact spectral indices are 

-ng = l-n = 2/{p-l) (5.61) 

The normalisations of the spectra do not deviate much from the slow-roll expressions. Prom 
Eq. (|4.24| ) it follows that the same is true of the density contribution to the microwave 



anisotropy. Pull results have not been given for the gravitational contribution, but for £ = 2 
the value of Rp/12.5 can be read from Pigure 2 of Pabbri, Lucchin and Matarrese (1986), 
and is within 10% or so of unity for .6 < n < 1. 



Extended Inflation. 

Extended inflation is based on a modification to the gravitational sector of the Lagrangian, 
which allows a first-order inflationary phase transition to complete satisfactorily (La & 
Steinhardt 1989; Kolb 1991). The modiflcation can be eliminated by a conformal transfor- 
mation of the metric. In general this is not a very useful thing to do because it modifies the 
Lagrangian of the non-gravitational fields, but it is useful during inflation. In fact, after 
the transformation one has power-law inflation, with unmodified gravity (Kolb, Salopek & 
Turner 1990). 

The original model (La & Steinhardt 1989) was based on a Brans-Dicke theory with 
parameter lo, and although this proved insufficient to allow present day tests of general 
relativity to be satisfied, it has remained the paradigm around which more complicated 
working models are based. After the conformal transformation, one has power-law inflation 
with 2p = uj + 3/2. 

A crucial difference between power-law and extended inflation is that extended inflation 
suffers an additional constraint, as one must avoid the large bubbles generated as inflation 
ends from being so profuse as to unacceptably distort the microwave background. This 
constrains w as a function of the inflaton energy scale M as (Liddle &: Wands 1991; Liddle 
& Lyth 1992) 

< 20 + 0.71ogio(M/mp/) (5.62) 

Bounding M using the microwave background limits on the fluctuation amplitude gives 
Lv < 17, corresponding to n < 0.75. We have discussed the signiflcance of this bound in 
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earlier papers (Liddle & Lyth 1992, 1993) and the prognosis for the extended inflation model 
is not good, as we shall recap in the Conclusion. 



Natural Inflation. 

The natural inflation model (Freese, Frieman & Olinto 1990; Adams et al 1993) is based on 
a pseudo-Nambu-Goldstone boson evolving in a potential 

y(,/.) =A4(1±cos(</.//)) (5.63) 

where A and / are mass scales. In this model 

e = r-^ tan^ (5.64) 

r] = -r-i (^1 - tan^ (5.65) 

where r = (167r/^)/mp;. Since e — r] = one must have r 1 to be in the slow-roll regime 
anywhere on this potential. On the other hand, one would like to have / < (?7ip^ 

/87r)i/2 so 

that quantum gravity does not spoil the field theory on which Eq. ( 5.63| ) is based. To this 



extent 'natural' inflation is somewhat unnatural, and one will not wish to have 1/r smaller 
than is necessary to obtain a viable model. 

Inflation ends with an oscillatory exit at the epoch given by Eq. ( 5.12[) , tan(i?i>cnd/2/) — 
which is of order 1. Using this result, Eq. (5.20) with A'^i = 60 gives (pi/2f ~ 
exp(— 60/r), leading to 

= le-i20A (5.66) 
r 

rji ~ —1/r (5.67) 
1 - n ~ 2/r (5.68) 

Thus, natural inflation makes the gravitational waves negligible but tends to give a signifl- 
cantly tilted spectrum. 

A similar analysis is valid for any model which has the inflaton fleld near the top of an 
inverted harmonic oscillator potential, V{(l)) = Vq — ^m?(lP. 

That natural inflation gives the same spectrum as power-law inflation may seem sur- 
prising in the light of claims that one can reconstruct inflationary potentials from a given 
spectrum (Hodges et al 1989; Hodges & Blumenthal 1990). In fact, these two models can be 
regarded as different regimes of an all-encompassing potential (which one can calculate in 
the manner of Hodges and Blumenthal) which is essentially 1 / cosh^ (f) with various factors 
thrown in (Copeland et al 1993). In the ~ region we have the inverted harmonic oscilla- 
tor, while at large (f) we have the exponential region. In a sense, this is the unique potential 
from which all inflationary models giving power-law spectra arise. With this potential, the 
scalar field can roll from the top down to the exponential region, while in a slow-roll ap- 
proximation generating an exact power-law spectrum. However, one must also note that 
as discussed above, the fact that the spectrum has been tilted implies that the slow-roll 
approximations are at best only just satisfied, so there will be corrections to the slow-roll 
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spectrum in all regions of this potential. In the exponential region these corrections have 
long been known to affect only the amplitude and not the slope; Stewart and Lyth (1993) 
have shown that the corrections are also not very significant near the top of the potential. 



Two-scale inflation 

The two-scale model written down by Linde (1991b) corresponds to the potential 

V{i;, (t>) = iA(V'2 - M2)2 + im2</,2 + (5_g9) 

The couplings A and A' are supposed to be somewhat less than unity, and in the numerical 
examples we set A = A' = 0.1. For (j)^ > 0^^^^ = AM^/A', the potential for the ip field has 
a local minimum at ip = 0, in which the field is assumed to sit (except for the quantum 
fluctuation). Inflation can occur as (j) slowly rolls down its potential 

V{(l)) = ^XM^ + ^m2(/.2 (5.70) 

When (j) falls below (pend this local minimum becomes a local maximum, so that there is 
a phase transition. Provided that bubble formation is negligible the transition will be of 
second order, the fields rolling quickly down to their vacuum values ip = M and (p = 0. 
(Otherwise we have a first order exit model akin to the two-field first-order models of for 
instance Adams and Freese (1991) and Linde (1991a), which are briefiy discussed in the 
following subsection.) 

We will suppose that while observable scales leave the horizon the first term of Eq. ( 5.70 ) 
dominates, since in the opposite case we recover the (f^ potential already considered. This 
means that 

<-) 

Of the parameters e and ij which are required to be small, the second is independent of (p, 

r? = (5.72) 

where 

The ratio X must therefore be significantly less than 1. From Eq. ( |5.20| ), 

0^ = ^lx/X'Me^^i (5.74) 
Consistency with Eq. ( 5.71| ) requires roughly rj < .1. The other small quantity ei is given 

by 

and Eq. ( ^.71| ) requires ei <^ rj. It therefore follows that the spectral index n is bigger than 
1 in the two-scale model. 
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Setting A = A' = .1 and imposing the COBE normalisation 5h = ^-7 x 10^^ determines 
all of the parameters in terms of m. For m = 100 GeV one has M = 4 x 10^^ GeV leading 
to r] = 10~^ and ei = 10~'^^. The gravitational waves are absolutely negligible, and the 
spectral index is extremely close to 1. The maximum value of m permitted by Eq. ( [5. 711 ) is 
roughly m = 10^^ GeV, giving M = 2 x 10^^ GeV, r/ = .07 and ei = 10~^. The gravitational 
waves are still negligible, but n = 1.14, significantly bigger than 1. 

The low end of the permitted mass range, m = 10^ to 10^ GeV, is particularly interesting 
because it corresponds to the Higgs field(s) of the standard model. Moreover, the other 
mass scale M is very roughly of order mmpi/(87r)^/^, which is the 'intermediate mass scale' 
invoked in supergravity models (Nilles 1984). 

So far we have considered a single scalar field ip, possessing the discrete symmetry 
ip — > —ip. When inflation ends, domain walls will form along the surfaces in space where ^ 
is exactly zero, so to make the model cosmologically viable one would have to get rid of the 
walls by slightly breaking the symmetry. However, it is clear that one can easily consider 
more fields, with continuous global or gauge symmetries. In particular, if ip is complex one 
can use the same potential with the replacement tp'^ IV'P- This leads to exactly the same 
inflationary model, but now global strings are formed instead of domain walls. A particular 
realisation of this case might be to identify ip with the Peccei-Quinn field. 

Although the two-scale model seems to be quite attractive, one should bear in mind that 
unlike the other models listed it has not yet been the subject of detailed public scrutiny. 

Double field inflation 

The double field inflation model (Adams & Freese 1991; Linde 1991a) is very similar to the 
two-scale model just discussed, whereby one scalar field (say cp) is utilised to trap a second 
(say tp) in a metastable state for a finite amount of time, terminated due to some physical 
eff'ect of the evolution of (p. The difference is that in this model the transition is to be a 
first-order one proceeding through bubble nucleation, so the role of the second scalar field 
is to alter the nucleation probability without allowing the barrier to disappear entirely. As 
the nucleation probability is typically exponentially sensitive to the details of the barrier, 
it is certainly possible for this to be achieved for very modest evolution of (p. 

This model is interesting in that it does appear to permit infiation to end via a first- 
order phase transition, which variants on the extended inflation paradigm seem unable to do 
(Liddle & Lyth 1993) under current constraints (working extended models all rely on ending 
inflation by some dynamical means other than bubble nucleation, or by erasing the effects of 
bubbles by a subsequent inflationary epoch). This is because of the considerable flexibility 
in modifying the nucleation rate permitting one to avoid the danger of any excessively 
large bubbles. However, in practice, the interactions of the two fields must be exceptionally 
finely balanced (Adams & Freese 1991), and so these models should be regarded as rather 
unnatural. 
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Intermediate Inflation 



Intermediate inflation (Barrow 1990; Barrow &: Saich 1990) is a model in which the scale 
factor expands according to 

a = exp (At^) (5.76) 

where A and / are constants with < / < 1. It arises as an exact solution for the rather 
complicated potential 



(/? + 4)2 



(2^/3)1/2 



-/3 



(5.77) 



with /3 = 4(/~^ — 1), and rather more interestingly as a slow-roll solution for a potential 
falling off asymptotically as a power-law V{(j)) oc 

Although rather contrived, this potential has some interesting properties differing from 
those we have already seen (Barrow &: Liddle 1993). Most interestingly, in the case / = 2/3 
the density perturbation spectrum is (in the slow-roll approximation) exactly the Harrison- 
Zel'dovich case. Indeed, because this is brought about by trading off the e and rj slow- 
roll parameters, there is no automatic requirement that the gravitational wave production 
be small despite the flatness of the spectrum. The second interesting aspect is that this 
potential provides a second example of n > 1, which arises with 2/3 < / < 1. 

In practice, a substantial gravitational wave contribution (or n significantly greater than 
one) does prove hard to achieve (Barrow &: Liddle 1993). The drawback of the model is 
that (as with pure power-law inflation) there is no natural exit from the model — one 
needs a mechanism to break the power-law form of the potential. Indeed, in this model 
the slow-roll parameters decrease as the field rolls down the potential, giving smaller and 
smaller deviations from flatness with negligible gravitational waves as it proceeds. Only by 
fine-tuning the timing of the mechanism which ends inflation can one arrange for presently 
observable scales to leave the horizon when the field is placed so as to give large deviations 
from the standard predictions. 



6 Large Scales: Normalisation from the COBE data 

For the rest of this report we are mostly concerned with the comparison between theory 
and observation. We want to know what observation can tell us about the spectrum 5h, 
on various scales, and start here with the largest scales which are explored by the cmb 
anisotropy. We remind the reader that we are expressing our spectral normalisations via 
CTg, using the CDM spectrum to translate between scales. 

6.1 Calculating the observed quantities 

As we write, the only positive detection of anisotropy is the COBE result, announced in 
April 1992 (Smoot et al 1992). f^. However, there are many interesting microwave results 
expected in the near future. Before focussing on the COBE data, we give here some generally 

^^As we were revising for the final version, a Princeton-MIT balloon experiment operating on similar 
angular scales produced a positive result with amplitude similar to the COBE 1-sigma lower limit. 
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applicable formalism, following closely Scaramella and Vittorio (1990). [See also Abbott & 
Wise 1984a; Bond & Efstathiou 1987; Scaramella & Vittorio 1988; Efstathiou 1990, 1991.] 
By analogy with Eq. ( |3.4| ), one can define a temperature correlation function, 

'AT(e) AT(e') 



C{a) 



T 



T 



(6.1) 



Here e and e' specify the directions in which the anisotropy is observed, and the average 
goes over directions separated by an angle a. It is given in terms of the multipoles by 



(6.2) 



where 



Qt 



1 

4tt 



+1 

rn=—l 



m\2 



(6.3) 



Note that our Qf has a different normalisation to that of Scaramella and Vittorio (1990). 

We introduced earlier the notion of a filtered density contrast 6{Rf,x), obtained by 
smearing over a region whose radius is of order Rf,oi equivalently dropping Fourier modes 
with k>RJ^. The same notion is useful for the microwave background anisotropy, so that 
one can smear A{T)/T over a region of angular size 9f (<C 1 radian), which is equivalent 
to dropping multipoles with / > 9J^. For a Gaussian filter with dispersion 6f (equal to 
.425 times the full width at half maximum) the precise statement is that each multipole is 
reduced according to the formula 



m\2 



exp 



-((2/ + 1)0^/2)^ (a[ 



Associated with the smeared quantity is a correlation function 

oo 

C{ef,a) =^exp \-{{2l + l)a/2f] QfPi{cosa) 

1=2 



(6.4) 



(6.5) 



A given experimental setup typically measures something which can be directly related 
to C. The simplest one is a single beam whose resolution can be represented by a Gaussian 
profile. Averaged over the sky the anisotropy measured by such a beam is 

''\=C(0/,O) (6.6) 



AT(0/,e) 



For more complicated setups involving two or three beam switching, still with Gaussian 
profiles, the measured anisotropy is 

' ) =2[C{ef,0)-C{9f,a)] (6.7) 



AT{ef,e) Ar(%,e' 



and 



AT{9f,e) 



l AT{9f,e') 
2 T 



l AT{9f,e") 
2 T 



^Ci9f,0)-2Ci9f,a) + ^Ci9f,2a) 



(6.8) 
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In the second expression, e' and e" lie on opposite sides of e, aligned on a great circle and 
each at an angular distance a. In a typical setup the 'beam throw' a is of the same order 
of magnitude as the 'antenna resolution' 6f. 

If we denote the left hand side of Eq. ( |6.6D , ( |6.7] ) or ( |6.8| ) generically by (AT/T)^, the 
prediction in terms of multipoles may be written 



where Fi is a filter function (Bond & Efstathiou 1987; Efstathiou 1991). The filter function 
is normalised to 1/2, rather than 1, in the limit of no filtering, but we have retained this 
normalisation for compatability with earlier work. 

These expressions involving Gaussian beam profiles are accurate representations of some 
experimental setups, including that of COBE, but require modification in some other 
cases.0 Accurate filter functions for some currently mounted observations are shown in 
Figure 4, which is reproduced from the paper of Crittenden et al (1993). 

The cosmic variance 

The multipoles a™ which appear in the above expressions depend on the position of the ob- 
server. We saw in Section ^ that the probability distribution of each multipole as a function 
of the observer's position is Gaussian with zero mean, and that cosmological perturbation 
theory can predict the variance of this distribution, which is the expected value of 
It thus predicts the expected values of the observed anisotropics defined by Eqs. (|6.3|), ( |6.5[ ), 



In comparing theory with observation, it is important to appreciate that what one can 
actually measure is only the temperature anisotropy at a single observer point, our own. 
Thus in practice one only gets a single realisation of the Gaussian probability distribution. 
This leads to statistical uncertainties, as an individual realisation may not reflect the prop- 
erties of the ensemble averaged system. This effect, normally called the cosmic variance, is 
significant for observables which depend only on a limited number of the a[", and in partic- 
ular for the quadrupole Q2 which only depends on the five a^^ (two of which are redundant 
rotational information). This is the fundamental reason why the quadrupole measurement 
is not particularly useful for constraining theories. On the other hand, the anisotropy on a 
scale of a few degrees depends on a significant number of the and is far less susceptible 
to statistical vagaries. 

6.2 The COBE quadrupole anisotropy 

The COBE observations yield enough data to reconstruct the temperature fluctuation field 
across the entire sky. We focus first on the quadrupole, though we shall find that it is not 
particularly useful for constraining the models, and we shall ultimately drop it in favour of 
the 10^ result. Nevertheless, it is worth examining the predictions of the model to see why 
we come to this conclusion. 

^^We are indebted to George Efstathiou for a clarifying conversation on this issue. 




1=2 



(6.9) 



(U), (P) and {^). 
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The above equations readily allow us to calculate the mean of the predicted distribution 
function for It is given in Figure 5, again as a function of spectral slope, for both 
types of inflationary model. Because the spectral normalisation is uncertain by the bias 
factor, we have plotted the quantity (Qi/'^s)- '^'^^ COBE result of 2.5 x 10~^^ corresponds 
almost exactly to the mean prediction for the flat spectrum when the bias is one, whereas at 
n = 0.6 the predicted mean is nearly 12 times the COBE result even for the natural inflation 
scenario, seemingly requiring a bias of over three. However, that conclusion neglects the 
statistical nature of the quadrupole prediction. 

As discussed above, theory does not predict a unique values for the quadrupoles Qi, 
but rather a probability distribution. Each quadrupole is a sum Eq. ( |6.3|) of the square of 
2/ + 1 quantities with Gaussian probability distributions, and its probability distribution is 
called a X(2i+i) distribution. For a given observer, the prediction is for a given realisation 
from that distribution. Let us discuss the quadrupole predictions in terms of a quantity 
q = 10^° 6^ Qf. Then the mean value q = W^^ 6^ (QDj xK*?) pdf with that mean is 

X'M = „ q'^' exp f-l) (6.10) 



3V2^{q/5f "V 2q 

This distribution has a very broad spread (the variance being 2/5 of the mean squared), 
and thus realisations of it may differ considerably from the mean. The spread is indicated 
in Figure 5 by the vertical bars through the mean predictions. 

Let us for the time being assume the COBE measurement to be perfect {ie not subject 
to observational errors). Then one can only exclude a given theory on the basis of some 
exclusion level, where the experimental result is far along the tail of the distribution. A 
conventional choice would be that 95% of the distribution should predict values above 
(or below) the experimental measurement. Thus one can exclude at 95% confidence only 
those theories which predict a mean sufficiently high that 95% of the distribution is above 
the COBE value of 

'Zexp — 0.25 6^. This one can readily calculate via error functions, 
to discover that one can only say with 95% confidence that q < 1.16^, corresponding to 



(Ql) < 1.05 X 10-^6. 

One can loosen this constraint even further by incorporating the fact that the COBE 
results possess experimental errors. One can model the COBE data with a probability 
distribution function describing the expected values were the experiment to be repeated — 
a sensible choice might be to assume to be gaussian distributed with mean 5 x 10^^ 
and width 1.5 x 10~^ — and use this to construct a pdf PexpiQexp) for the experimental 
measurement ^exp- One then tests each member of the experimental distribution for exclu- 
sion and takes a weighted mean, thus constructing a rejection functional on the predicted 
theoretical pdjh as 

dqexp (6.11) 



= / Pcxp(fep) 1 - 2 / X^{(l)dq 
Jo Jo 



where | signals the modulus. Defined on one side of the probability distribution like this, 
a value TZ > 0.9 signals a 95% rejection of the theoretical pdf in the light of the modelled 
experimental data. With the experimental modelling as suggested above, predicted mean 



values of q up to 2.36^ are allowed, corresponding to 95% confidence that y (Qi) < 1-5 x 
10-^6. 
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Hence one sees that the statistical uncertainties in the theoretical quadrupole prediction 
make it of little use in constraining theories. For example, with this last constraint n = 0.6 
is allowed with a mild bias of 1.15. 



6.3 The variance at 10° 

When one carries out an experiment based on a particular beam configuration, the ex- 
perimental configuration typically involves averaging over fluctuations on scales below the 
beam resolution. This cuts off the contribution from higher multipoles. From Eq. ( |6.9| ), the 
expected value of the anisotropy seen in any particular experiment is 

((^)} = ^E(2^ + 1)S?^^ (6-12) 

For ground-based experiments which typically feature two and three beam configura- 
tions, this filter function can be rather complex, as we saw earlier. For a single beam 
experiment like COBE it is much simpler. The original COBE beam is well approximated 
by a gaussian with Full Width Half Maximum of 7'^. However, before calculating the vari- 
ance they smooth again by convolving with a further 7^^ FWHM gaussian, a procedure 
equivalent to an original smoothing by a 10*^ FWHM gaussian. Such a gaussian has a 
variance a"^ = (4.25'^)^. Thus the appropriate form of the filter function is 

Fi = ^ exp (- (4.25^(/ + l/2)/180)2) (6.13) 



It is now trivial to calculate the predicted lO'^ variance as a function of n, and in Figure 
6 we plot its square root, henceforth denoted ^ multiplied by the bias, for both power- 
law and natural inflation. In the latter case, one sees that the results are essentially exactly 
linear. We have been unable to show why this should be analytically. The appropriate 
fitting function for natural inflation, which has no gravitational waves, is 



AT 



AT 

(n) = exp (2.62(1 -n)) — 

100 



(n = 1) (6.14) 

100 



while that for power-law inflation which does have gravitational waves is 

(n = 1) (6.15) 



AT 



iqo 



15 - 13n , , AT 
{n) = J exp (2.62(1 -n)) — 



iqo 



Remember that we are assuming the same normalisation at 8h^^ Mpc. If instead one makes 
the same normalisation to the COBE lO" result, then the n-dependent factor is transferred 
to the normalisation cig. 

One can also readily calculate the cosmic variance, giving the spread in values about 
these means that would be measured by differently positioned observers. The variance of 
the Qf is 2(Q2)2/(2/ + 1). For the 10° result, this is 10% at n = 1, rising to 12% at 
n=0.6 (this result remains true whether or not there is a gravitational wave contribution 
to the anisotropies) . This is a negligible correction to the larger COBE error in the present 
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observations, but is ultimately a limiting obstacle to the conclusions one can draw on this 
large angular scale. 

Readers may find useful a relation between the lO'^ variance and the expectation of the 
quadrupole as n is varied. For a power-law spectrum one has 

_ r(/ + (n-l)/2) r((9-n)/2) 
S2 r(/ + (5-n)/2) r((3 + n)/2) ^ ' 

For the cases we examine, this is also true of the gravitational wave contribution. The 
factor can now be pulled out of Eq. ( |6.12| ) to give an analytic expression for the ratio as 
a function of n alone, but in a rather awkward form. For the COBE result, a reasonable 
approximation in the range n G [0.5, 1] is just a linear fit 



2.03 + 0.68 (n - 1) (6.17) 



6.4 Microwave background constraints 

In establishing constraints on the normalisation erg = 1/bs at fixed spectral index n, it is the 
10^ data which are of primary interest. Some particular values worthy of note are that for the 

= 1.05x10-^1/68, and that forn = 0.6 of ^ = 5.17x IQ-^l/ftg 
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10" 



flat spectrum of ^ 



lOO 



for power-law inflation and -j^ = 2.99 x 10 l/bg for natural inflation. These of course 

are to be compared with the COBE observations of ^ = (1.1 it 0.2) x 10"^. 

The mean quadrupole prediction is also very similar to the flat spectrum prediction with 
bias one, but the statistical uncertainties make this comparison rather less relevant. It is 
worth recalling that an extrapolation of the quadrupole from the COBE data at smaller 
angles, assuming a power-law spectrum, gives a somewhat larger prediction for the mean 
quadrupole of 16 it corresponding to (6 ± 1.5) x 10^^ (Smoot et al 1992). Note that 

this extrapolation estimates the mean quadrupole, not a specific realisation. Hence this 
again favours values of the bias not much exceeding 1, unless one inserts large-scale power 
to boost the theoretical prediction for the mean. 

Lacking full access to the COBE data, one must make an operational choice as to what 
to take as the observational limits. The COBE team provide fits of their data to power- 
law spectra, where both the slope and amplitude are treated as fitting parameters. The 
two pertinent pieces of information are that the observed (root of the) variance at lO'^ is 
(1.1 lb 0.2) X 10"^, and the allowed slopes in the fit carried out by the COBE team are 
n = 1.1 lb 0.5, where both errors are 1-sigma. These are of course not independent, as lO'^ 
is at the lower end of the fitting range. Nevertheless, if one takes the upper 2-sigma value 
of the extrapolated quadrupole (24// K), then this actually gives a prediction for the lO'' 
variance which is well above the 2-sigma limit on the 10'^ data even for the central n value. 
This conclusion (also implicit in Efstathiou, Bond and White (1992)), argues that the fits 
are not rigid enough to be of much use in constraining theories. Because of the cosmic 
variance and also experimental uncertainties, the quadrupole itself also appears of little use 
in this type of analysis. Consequently, we choose to adopt simply the 10*^ result as it stands, 
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without further incorporation of the COBE data. We relax their error bars to 2-sigma, and 
for many purposes we are only interested in the upper limit thus given of 1.5 x 10^^. This 
is a particularly useful way of utilising the results, as it seems likely that this number can 
only go downwards if it is to avoid conflict with other experiments and so any limits quoted 
on its basis are likely only to become stronger with improved observations. 

The full limits on the bias as a function of n from this criterion, for both styles of 
inflationary model, are plotted in Figure 15 in Section ^ along with other constraints derived 
in later sections. Some sample results are that for n = 1 we have the obvious .7 < bs < 1.5, 
while for n = 0.6 we require 3.4 < bg < 7.4 (power-law inflation) and 2.0 < bs < 4.3 
(natural inflation). The following section will assess the validity of different values for the 
bias parameter. Let us finally recall the uncertainties of normalisation. The use of the 
Bardeen et al (1986) transfer function rather than the one we use would reduce the required 
bias by a further 10%, and a higher baryonic content raise it by a similar amount. 

6.5 Prospects for the near future 

It can be expected that COBE will release improved results once the data from later years 
has been analysed. The ground based Tenerife experiment (Watson et al 1992) operating 
on similarly large scales should also soon be in a position to confirm the COBE detection. 
Much interest has been focussed on the South Pole experiment operating on smaller angular 
scales of around a degree, which has recently announced very strong upper limits (Gaier et al 
1992). This is particularly interesting as those scales are below the horizon size at decoupling 
and so a gravitational contribution to the anisotropics will be suppressed. Indeed, with a 
fiat spectrum the South Pole upper limit is alarming low, and may indeed indicate that 
the COBE result has a significant gravitational wave contribution. However, the situation 
remains to be clarified. In particular, one would like to know whether a tilted spectrum 
without gravitational waves (a la natural inflation) is compatible with COBE and the South 
Pole, or whether the gravitational waves are necessary. Further, the COBE result could be 
towards the 2-sigma minimum of its range, in which case the conflict may disappear. While 
the data remain controversial, we do not feel that this paper is an appropriate place to 
investigate these issues. 

7 Medium scales: Galaxy Clustering and Bulk Velocities 

Now we consider scales which are big enough that the filtered density contrast is still 
evolving linearly, yet small enough that there exist observations of the galaxies and their 
motion with which one might hope to compare the linear theory. Thus we are discussing 
scales from around Wh~^ Mpc up to perhaps 100/i~^ Mpc. 

The cleanest observable to interpret theoretically is the galaxy peculiar velocity field 
smeared on scales > 10 Mpc, which is termed the bulk flow. It is generally assumed that the 
galaxy bulk flow will be the same as that of the underlying matter, which can be calculated 
from linear theory as discussed in Section ^.8. The other quantities which might hope to 
use are the number density perturbations 5N/N of specific classes of object, notably optical 
galaxies, infrared galaxies and galaxy clusters. On the assumption that they are related to 
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the underlying density perturbation 6p by a scale-independent bias factor, 

-J^ -^°bj— (7.1) 

they can be compared with the linear theory on scales r > lOMpc. Otherwise, their 
calculation takes us into the non-linear regime discussed in Section ^ even on scales r > 
lOMpc. 

The peculiar velocity field, if it is available, determines both the normalisation as and 
the spectral index n of the density perturbation. The number density perturbations, on the 
other hand, give information only about n even if the bias is indeed scale-independent. We 
discuss these two types of information now, in reverse order. 

7.1 The galaxy correlation function 

There is ever increasing evidence that the amount of large-scale galaxy clustering is greater 
than can be accommodated in the standard CDM cosmogony. The most striking piece of 
evidence is provided by the measurement of the galaxy angular correlation function w{6) in 
the APM survey (Maddox et al 1990, 1991), based on a sample of over two million galaxies. 
The angular correlation function measures the clustering pattern as seen in projection on the 
sky, and thus does not require redshift data, enabling such large catalogues to be obtained 
by automated measurements on photographic plates. With large redshift samples becoming 
available the further evidence for excess large scale clustering is mounting. Conclusions in 
broad agreement with APM are provided by the 'counts in cells' of the QDOT survey 
(Saunders et al 1991), and more recently by surveys of redshifts of APM galaxies (Dalton 
et al 1992; Loveday et al 1993), and in the power spectrum inferred from the CfA survey 
(Vogeley et al 1992), that inferred from the Southern Sky Redshift Survey (Park, Gott & 
da Costa 1992) and that from the 1.2 Jansky survey (Fisher et al 1993). 

If the bias parameter is indeed scale-independent, all of these data indicate that the 
CDM, n = 1 spectrum normalised at ~ lOMpc lacks power on larger scales. To 
quantify the amount of extra power required, Wright et al (1992) introduced a quantity 
which they called the excess power, which is a functional of the power spectrum defined as 

CT(25/i-iMpc) 

They chose the normalisation so that with their transfer function E = 1 for the CDM n = 1 
model; with our transfer function, E = 0.95 for that model. Another way of defining the 
required extra power is to give the value of n which would produce it (Liddle, Lyth &: 
Sutherland 1992). The connection between these two specifications is 

E[n] = l.U-n/2 ; 0.3 < n < 1 (7.3) 

To estimate the required value of E or n, let us analyse the APM data following Bond 
(1989) and Liddle, Lyth and Sutherland (1992). On the large angular scales where the 
discrepancy arises, the angular correlation function w{9) can readily be calculated in linear 
theory. It is given (Peebles 1980) by Limber's equation as an integral over the galaxy 
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correlation function £,gg{r) = b'^C{r), where the mass correlation function ^(r) is defined by 
Eq. ( |3.4[) . A linear calculation is expected to be reliable on scales above about 2^. The results 
appear in Figure 7, reproduced from Liddle, Lyth and Sutherland (1992). The flat spectrum 
n = 1 falls well below the observational data, but one can see that as the spectrum is tilted, 
the extra power does indeed make itself evident in the clustering statistics. It is suspected 
that the APM survey may contain small residual systematics which bias the observational 
estimates upwards (Will Sutherland, private communication; Fong, Hale-Sutton & Shanks 
1992), so it seems reasonable to regard values of n from 0.3 to 0.6 as good fits to the excess 
clustering data. The corresponding range for E — 1 is 0.29 to 0.14. We postpone for the 
moment the question whether the CDM model becomes viable in all respects, with a value 
of n in this range. 

A different way of specifying the required power, used for instance by Efstathiou, Bond 
and White (1992), is to specify the modification of the transfer function which would gener- 



ate it, through an alteration of the parameter F = ^Ih appearing in, for instance, Eqs. (2.45) 



and (2.48). The equivalent value of n is 



F = - 



1/ 1.88 



2 V 2.88 - n 



) (7.4) 



The range of F these authors consider a reasonable fit to the APM data is 0.15 < F < 0.30, 
corresponding to 0.15 < n < 0.67, which we see is actually rather looser than the range we 
took above. We postpone until Section |9| the question of whether this modification actually 
corresponds to altering the physical value of the present dark matter density 17 (with h 
fixed at say h = .5), and whether such a change would lead to a viable model. 

All of this analysis assumes that linear evolution is valid on scales > lOMpc. The 
assumption has been questioned by Couchman & Carlberg (1992), who claim that a highly 
evolved model with ug = 1 as suggested by the COBE data (though it is important to note 
their analysis preceeded COBE) can explain the APM data. This apparently occurs via 
a combination of non-linear effects altering the shape of the spectrum and the choice of a 
rather unusual method of identifying the galaxies in their simulations. Their analysis, which 
also gives (indeed requires) a significant velocity bias, has been questioned by subsequent 
authors, and is not generally accepted at present. 

Although the galaxy clustering data is usually taken to indicate the need for extra 
large scale power, the alternative explanation of a scale dependent bias factor is sometimes 
mooted. A particularly elegant proposal in this direction is to allow quasars to suppress 
local galaxy formation (Babul & White 1991). There is also 'cooperative galaxy formation' 
(Bower et al 1993), a phenomenological (at present) model in which galaxy formation is 
favoured in the neighbourhood of other galaxies. 



7.2 The bulk flow 

There are many controversial aspects to the measurements of the bulk fiow, but nonetheless 
they provide a useful measure of the absolute magnitude of the power spectrum on inter- 
mediate scales (for a recent review, see Kashlinsky &: Jones 1991). It is well known (Kolb 
& Turner 1990) that typical theories, including standard hot and biased cold dark matter 
models and also models seeded by topological defects, tend to predict bulk velocities rather 
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lower than those observed, particularly for standard CDM with high bias. Here we try to 
quantify this effect for cold dark matter, bearing in mind the uncertainties of theory and 
observation. 

Just as one can filter the density contrast, one can also filter the peculiar velocity field 
to obtain a bulk flow quantity, on any desired scale. As long the density contrast on that 
scale is in the linear regime, the bulk flow can be constructed from the linearly evolved 
Fourier components Vk, given by Eq. ( |2.68| ), and compared with observation. There are 
several different ways in which this comparison can be carried out. 

For our conclusions, the most important is the one which we consider first — an analysis 
of the correlation of the bulk flow with the QDOT density field. This provides strong 
constraints on the bias parameter for galaxies in the IRAS survey, and combined with other 
results from the QDOT survey concerning the clustering of IRAS galaxies leads directly to 
strong limits on the density field power spectrum. We also discuss the reconstruction of the 
full three dimensional density field via POTENT, and the cosmic Mach number. 

Correlating QDOT and peculiar velocities 

One certainly does not at present have a catalogue of galaxy peculiar velocities as large as 
one would like. This is due to the difficulties of requiring independent distance estimators 
in addition to redshifts. The most commonly used catalogue at present is that compiled by 
Burstein (from measurements by a variety of different researchers) containing around 1000 
galaxies, unfortunately distributed very anisotropically across the sky. 

This velocity catalogue has been used to estimate the density contrast on intermediate 
scales by the QDOT collaboration (Kaiser et al 1991). The QDOT survey (Saunders et al 
1991) consists of a redshifting of galaxies sparsely (1 in 6) sampled from the catalogue of 
infrared emitting galaxies compiled by the IRAS satellite. This survey has the advantage of 
even and nearly complete sky coverage, and allows the galaxy density field to be constructed. 

In order to obtain a value for the bias, this density field is correlated with the velocities 
from the Burstein sample. The redshift survey supplies the position of the galaxies in 
redshift space, which amounts to assuming that they have zero peculiar velocity along the 
line of sight. One attempts to invert this into the true galaxy distribution in real space. To 
do this, one constructs a smoothed galaxy number density contrast field, and appeals to a 
constant bias to generate the density field from it. This density field allows one to calculate 
the expected galaxy peculiar velocities, and use this to correct their position. With luck, 
this generates an improved galaxy distribution, on which the whole process can be repeated 
iteratively. This success of this iteration will depend on the accuracy of the choice of bias 
parametei]^, as an incorrect choice will systematically incorrectly estimate the density field 
and hence the peculiar velocities. Thus the correlation with the observed peculiar velocities 
constrains the bias. 

In general, the bias of the infra-red selected (and thus typically young) IRAS galaxies will 
not be the same as that of the optically selected galaxies discussed thus far, and so we denote 
this bias by 6/. It is often stated that IRAS galaxies are somewhat less clustered than their 
optical counterparts; for example Saunders, Rowan-Robinson and Lawrence (1992) suggest 

^^And in general on f2, though the inversion here depends only on the combination OP'^ /h and we are 
anyway presupposing Q. = 1. 
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bi = (0.69 lb 0.09)bg at the 1-sigma level. Kaiser et al (1991) obtain from their analysis 
bf/Q^'^ = 1.16 lb 0.21. This provides the first dynamical evidence which is manifestly 
consistent with Q = 1, as we are assuming here. 

There are other ways of constraining bj. According to Taylor and Rowan-Robinson 
(1992), there are presently three independent dynamical estimates of bj (one being that 
described above), giving respectively 6/ = 1.23 ib 0.23, 1.16 ib 0.21 and 1.2 ib 0.1 where all 
errors are 1-sigma. The first is obtained from the convergence of the velocity dipole, the 
second as discussed above and the third from the covariance tensor of reconstructed flow 
fields. One cannot statistically combine errors from the same data set, so we adopt as a 
1-sigma result bj = 1.2 ib 0.2. 

This estimate of bj can be turned into an estimate of fig (Efstathiou, Bond & White 
1992), because the QDOT survey supplies the dispersion of counts of QDOT galaxies in 
30/i~^ Mpc cubes, given as 0.46 ib 0.07 (Saunders, Rowan-Robinson &: Lawrence 1992) 
(we have doubled the quoted error bars). Because we have a good estimate of the bias, 
we can immediately calculate the mass dispersion in these cubes, and so determine the 
normalisation of the spectrum as defined by as- The result of this calculation is shown in 
Figure 11. 



POTENT reconstruction of the velocity field 

Another simple comparison which can be made is to the mean peculiar velocity of the 
galaxies in a sphere of radius R around us. For R > 10 Mpc this should be equal to the 
linearly evolved filtered peculiar velocity field v(i?, x), with a top hat filter of radius R 
and evaluated at our position. The expected value of its modulus squared is therefore the 
variance <J^{R) of the smeared velocity field, which from Eqs. ( 2.68| ), (3.3) and ( 3.12 ) is 
given by 

a^iR)=Hll^W{R,k)^^ (7.5) 

Here, W{R,k) is the top hat filter W/Vf defined by Eq. ( |3.22| ) Because of the k~'^ factor, 
the infiuence of large scales filters down significantly to small scales. 

Using this approach it is easy to consider the effect of n. Figure 8 illustrates the scaling 
of the dispersion with Rf for a top hat filter. Smaller values of n give significantly higher 
predictions, but the scaling with Rf is roughly the same. 

To compare this result with observation we use the velocity field reconstruction results 
from the POTENT method, pioneered by Bertschinger, Dekel and collaborators (Bert- 
schinger &; Dekel 1989; Dekel, Bertschinger & Faber 1990; Bertschinger et al 1990; Dekel 
1991). This is a very powerful technique, at present rather frustrated by the inadequacies 
of available data, which allows one to take the radial component of the peculiar velocity, 
which is all one can measure, and reconstruct the full three dimensional velocity fiow. 
The procedure is to assume that the velocity field is given as the gradient of a scalar (in 
practice this will be identified as the peculiar gravitational potential of Section §.8). This 
scalar can be calculated as a line-of-sight integral requiring only the radial peculiar velocity 
component; the full velocity field is then obtained by taking the gradient. The power and 
potential of this method has been vividly demonstrated by tests on simulated data, but 
on real data problems of sparseness and uneven sampling present awkward obstacles to a 
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reconstruction of the true velocity field around us. However, although the detailed form 
of the reconstruction is open to question, it certainly provides useful measures such as the 
bulk flow in spheres about us. 

It is not trivial to compare theory with observations, because one predicts only the rms 
velocity averaged over the entire universe and the observed distribution may not be a fair 
sample (phrased another way, the comparison is subject to a large cosmic variance due to 
the limited amount of data available; in the language of Section one is dealing with a xi 
distribution since v'^ = v'^ -\- Vy -\- v'^) . In particular, the possible contaminating effects of 
the great attr actor (if backside infall can be unambiguously identified (Mathewson, Ford & 
Buchhorn 1992)), may distort observations. In the Monte Carlo work analysing bulk flows 
of Tormen et al (1992), special criteria for choosing the observer points are employed before 
statistical comparisons are made. 

For this statistic, it is particularly notable that the influence of large scales filters down 
significantly to small scales, so an unexpectedly large fluctuation (and the great attractor 
would certainly be one in these CDM models) can easily move the whole observed data set 
to well above the mean, as predicted from the spectrum, that would be observed in a fair 
sample. That is, although one might think that there were many independent samples on 
the scale at which one is looking, the dependence on large scales correlates the supposedly 
separate samples leading to a cosmic variance well in excess of that one might have expected. 

The POTENT results are obtained via a two-stage smoothing (Dekel 1991). First, 
the original data is smoothed with a gaussian of radius 12h~^ Mpc to remove short-scale 
sampling difficulties, and then a top hat of radius -R/ is used to provide the published 
data. In Figure 9, we apply this two-stage smoothing to our spectra, which reduces the 
predictions, especially at short scales. These are compared with the POTENT data at 
different normalisations. It is seen that a normalisation ug bigger than .7 is preferred, in 
keeping with the conclusion that we reached earlier. This conclusion holds more or less 
independently of the shape of the spectrum, say in the range .7 < n < 1, because the scale 
probed by the bulk flow data is not all that bigger than the normalisation scale. 

Briefly we note that an alternative comparison framework, using the Burstein sample of 
peculiar velocities and making a direct comparison of power-law models with observations, 
has been made by Tormen et al (1992). The utilise Monte Carlo simulations, and consider 
18 models combining all possible combinations of 63 = l/<78 = li 1-5, 2; n = 1, 0.5, 0; 
O = 0.4, 1. Of particular interest is the Q = 1, n = 0.5, 63 = 1-5 model. They employ a 
maximum likelyhood test, and find that this model is twice as likely as the model with flat 
spectrum and bias 1.5, and is four times more likely than the flat spectrum with bias 1. 
The analysis is however complicated by their requirement that the observer points in their 
simulations must be chosen to match the properties of the local group. 

The cosmic Mach number 

The cosmic Mach number test was devised by Ostriker and Suto (1990) in an attempt to 
probe the shape of the spectrum as opposed to its normalisation. They calculated a cosmic 
Mach number M.os{R) defined by 




(7.6) 
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where v{R, x) is the velocity field smeared on a scale R and () denotes the average over 
all positions x. A large value of M{R) (greater than one) indicates a fairly coherent (cold) 
flow on the given scale, while a small value indicates a fairly random (hot) flow. 

In order to be (marginally) within the linear regime they smeared v on a scale Rg = 
lOMpc before applying the definition, and then used it only for R 3> Rs- This gives, in 
terms of the window function W{R, k) used in Eq. (|7.5D, 



(j^ .p^^2 Jo-'W\R,)WHR)Vs{k)dk/k^ 

{J ios{^)) j;^w^Rs)[l- jW^{R)]V5{k)dk/k^ ^ 

(They also multiplied W'^{R) in the square bracket by (1 + k'^R^ /9) to allow for the fact 
that the spherically symmetric part of the peculiar velocity field around an observer will not 
be detected by local observations, but for simplicity we drop this small correction.) Clearly, 
A^OS depends only on the shape of the spectrum, not on its normalisation. 
For a Gaussian filter 

{a^{Rs)y - [a^iR + Rs)r 
and for any filter one has in the limit Rg <^ R 

Suto and Ostriker calculate A^os, and find that for CDM n = 1 model it is ~ 1 on a scale 
R = 20Mpc falling to ~ .2 at i? = 100 Mpc. Increasing the power on large scales increases 
A^OS) but the effect is not very strong with a reduction to n = .5 increasing Aios by at 
most a factor 1.5. 

From the data that we presented and discussed in Figure 9 it is clear that TWqs is 
difficult to estimate from published observations, especially when one remembers the cosmic 
variance. Such a direct comparison was not attempted by Ostriker and Suto. Instead, they 
essentially considered the observational value of a different quantity, which corresponds to 
a more conventional definition of the Mach number 



M{R) = v{R)/a (7.10) 

where v{R) is the mean, and a is the dispersion, of the peculiar velocity in a sphere of radius 
R around us. From optical galaxy data they estimated (see also Groth, Juszkiewicz and 
Ostriker (1989)) Al(16Mpc) = 2.2 ± .5 and Al(36Mpc) = 1.3 ± .4. If one were prepared 
to regard these as estimates of A^OS) a lot of extra large scale power would be indicated, 
corresponding to a value n ~ 0. Subsequent work (Suto & Fugita 1990; Suto, Gouda & 
Sugiyama 1990; Suto, Cen & Ostriker 1992) confirmed this, and also claimed on the basis 
of statistical analysis that M and Aios can be identified within a factor better than 2 with 
reasonable confidence. 

Our conclusion is that the Mach number test seems to point in the same direction as 
the galaxy clustering data, in that more rather than less large-scale power is indicated, but 
that the latter data are much more decisive at present. 
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8 The non-linear regime 



In the last section we considered the galaxy bulk flow, and the correlation functions of 
galaxies and galaxy clusters, on scales > lOMpc. The bulk flow can be compared directly 
with linear theory, and so can the correlation functions on the assumption that there is 
a scale independent bias factor. Most observational quantities cannot be calculated using 
only linear theory. A partial list of them is the following. 

• The redshift of formation of a given class of objects, such as optical galaxies, infrared 
galaxies or galaxy clusters. 

• The correlation function ^obj(^) of such a class, and therefore the bias factor 6obj 



defined by Eq. (7T) to the extent that it is a valid concept. 



• The number density n{z, > M) of all gravitationally bound systems with mass bigger 
than M, which exist at redshift z. 

• The mean virial velocity of a given class of objects. 

• The number density n{z, > v) of all gravitationally bound systems with virial velocity 
bigger than v, which exist at redshift z. 

• The dispersion cj|| (r) of the relative line-of sight velocity of a pair of galaxies separated 
by distance r. 

Although none of these quantities can be calculated using only linear cosmological per- 
turbation theory, some of them can be estimated if the linear theory is combined with a 
spherical model of gravitational collapse. This quasi-linear approach can be used in different 
ways. At the crudest level its use is relatively uncontroversial; it allows one to show that 
structure forms in a bottom-up manner, with light objects forming first, and it allows one 
to estimate the epoch of formation for objects of given mass. These two applications form 
the subject of the Sections 8.1 and 8.2 below. 

If one is prepared to take the model more seriously, it is possible to make more definite 
statements, by identifying the high peaks of the density contrast in the linear regime with the 
sites of structure formation. This quasi-linear approach has a long history, a very influential 
early paper being that of Press and Schechter (1974). Then, in another influential paper 
Bardeen et al (1986, henceforth BBKS) developed it in considerable detail, at the same 
time applying it to the CDM model. It's defect, as BBKS and many other authors have 
pointed out, is that estimates made from it depend sensitively on at least two parameters 
(the threshold density contrast, and the ratio of the filtering mass to the mass of the objects 
studied which is usually set arbitrarily equal to 1) whose values are quite uncertain; or to 
be more precise, are ill-defined except under very idealised assumptions. 

Because of its widespread use, we explore the quasi-linear approach in detail in Sections 
8.3-8.8. The basic assumptions are explained, then the essential results of Press-Schechter 
and BBKS described, and finally the sensitivity of the results to the unknown parameters 
will be quantified. The conclusion will be that the quasi-linear approach is of somewhat 
limited use, its main value being to exclude very low normalisations of the spectrum on 
scales < IMpc. 
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Finally in Section 8.9 we go on to consider the results of numerical simulations, which 
avoid using the quasi-linear approach or even the spherical collapse model, by modelling 
in various ways the non-linear evolution of the dark matter and baryons. Such simula- 
tions appear to be capable of giving reasonably definite information, and indicate that the 
observations cannot be fitted by the CDM model if it is normalised to fit the COBE data. 



8.1 Gravitational collapse 

The quasi- linear approach makes essential use of the filtered density contrast 6{M,x.). Fol- 
lowing the usual practice for theoretical calculations, we use the Gaussian filter, defined by 
Eqs. (3.16) and ( 3.19| ). The dispersion of the filtered density contrast is thus given by 



c^HM) = J^ Vs{k)eM-R}k^)Y (8.1) 

At this point we need to distinguish between the density contrast of the cold dark matter, 
and that of the baryonic matter. After decoupling the baryons fall into the potential wells 
created by the cold dark matter, and achieve the same density contrast, on scales in excess 
of the Jeans mass of order 10^ Mq (Peebles 1980, 1984a). On smaller scales their pressure 
prevents the baryons from collapsing, so that their density remains practically uniform. 
Thus, luminous galaxies have mass in excess of W^Mq according to the cold dark matter 
model, which is in rough agreement with observation. (Note though that, as discussed in 
for example the review of Blumenthal et al (1984), mass scales of order W^Mq to W^Mq 
are also distinguished by another consideration; they are the smallest mass scales on which 
the baryons can efficiently radiate energy, to become the rather tightly bound structures 
that are observed.) From now on we identify the density contrasts of the cold dark matter 
and the baryons. 

As long as a{M) < 1, the filtered density contrast evolves linearly except in those rare 
regions of space where it exceeds 1 in magnitude. In the 50% of these regions where it is 
positive, gravitational collapse takes place. At least initially, the collapsing regions have 
mass > M, because the filtered density contrast does not 'see' structure on smaller scales. 
In the approach that we are considering, gravitational collapse is modelled by taking the 
collapsing region to be spherically symmetric. One then finds (Peebles 1980, Eqs. 19.50 and 
19.53) that when a given mass shell stops expanding, the mean density inside it is a factor 
97r^/16 times bigger than the mean density of the universe. If the matter inside the shell 
had evolved linearly (density contrast oc (1 + z)~^), its density contrast at that time would 
have been 6 = (3/5)(37r/4)^/^ = 1.06. After it has stopped expanding, the shell collapses. 
If spherical symmetry continues to hold, and one neglects pressure forces, the shell has 
collapsed to a point by the time that the age of the universe has increased by precisely 
a factor 2, corresponding to a mean density contrast within the linearly evolved shell of 
S = (3/5)(37r/2)2/3 = 1.69. Numerical studies (eg. Peebles 1970) indicate that by about 
this time, pressure forces will in fact have virialised the random motion of the constituents 
of the object. After this initial virialisation, the object can lose energy (dissipation) which 
further increases its virial velocity. 

The use of the spherical collapse model in the present context is somewhat problematical. 
First, filtering on a mass scale M will distort the profile of a peak with mass of order M, 
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making it lower and broader. Thus the hnearly evolved filtered density contrast will not 
have precisely the behaviour described above, even for a spherical peak. Second, a peak 
of the filtered density will not be spherically symmetric, and the corresponding peak of 
the unfiltered quantity will be even less so. Finally, any initial departure from sphericity 
is amplified during the collapse (Peebles 1980). Nevertheless, one might hope that two 
features of the spherical collapse model translate into reality. 

• The regions with mass > M which have undergone gravitational collapse can be at 
least approximately identified with the regions where the linearly evolved density 
contrast S{M,:x.) exceeds some threshold 6c which is roughly of order l.P| 

• A collapsing region does not fragment into a large number of separate objects, which 
means that the mass of the resulting gravitationally bound systems is also > M. 

These two assumptions form the basis of the quasi-linear approach, and they will be 
adopted for the rest of the discussion. 

8.2 The epoch of structure formation 

The quasi-linear approach that we are going to discuss is valid only in the regime in which 
the mass density is evolving linearly, which ends when the filtered linearly evolved dispersion 
a{z, M) of the density contrast is of order 1. With the usual n = 1 spectrum it is known that 
this is the epoch of structure formation for objects of mass M, so that one has a bottom-up 
scenario in which low mass objects form first. In this section we derive this result, and ask 
how far it remains valid if n is less than 1. 

For definiteness, let us take the epoch when non-linearity sets in to be the epoch when 
the linearly evolved quantity a{z,M) is precisely equal to 1. Denoting this epoch by z^i 
one has 

1 + Zni(M) = ao{M) (8.2) 

where as always the subscript denotes the present value of the linearly evolved quantity. 
At about this epoch, the regions where gravitational collapse is taking place on scales > M 
become common, whereas before they were quite rare. If a{M) is increasing significantly 
as the mass decreases, then at this same epoch the density contrast filtered on mass scales 
substantially bigger than M is still evolving linearly in most parts of the universe. The result 
is a bottom-up scenario, in which gravitationally bound systems of successively bigger mass 
M form at the successively later epochs Zni(M). After a system has formed, it may become 
part of a bigger gravitationally bound system, remaining a discrete object, but it may also 
merge with other systems. 

Let us ask if a{M) does, in fact, increase significantly as M decreases. For M > IO^Mq 
it has already been plotted in Figure 3b both for the standard choice n = 1 and for some 
smaller values. One sees that it increases quite fast in the regime M > IOMq, but flattens off 
a lower masses and in fact turns over for n significantly less than 1. This is a consequence of 
the small-scale behaviour Vs{k) oc T'^{k)k^~^^ oc fc""^ of the spectrum (up to a logarithmic 
factor). That behaviour is expected to persist down to some very small scale, called the 

^*Note that throughout this section the subscript c denotes the threshold value of 5, whereas in earlier 
sections it was used to distinguish the cold dark matter density contrast from the radiation density contrast. 
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coherence scale, below which the spectrum cuts off sharply and a{M) becomes practically 
equal to the unfiltered quantity a. The coherence scale depends on the nature of the cold 
dark matter. If the cold dark matter consists of subnuclear particles, the coherence scale is 
determined by their interactions. For instance, if it consists of the lightest supersymmetric 
particle the coherence scale is the Hubble scale at the epoch when the particle becomes 
non-relativistic. If it consists of the axion, the coherence scale is the Hubble scale at the 
epoch when the axion acquires mass through QCD effects. 

On the basis of these considerations, we arrive at the following picture of strTicture 
formation, which covers the case n < 1 as well as the familiar case n = 1. The first 
structure forms at some epoch given by 

1 + = cTo (8.3) 

where do is the unfiltered quantity. For n = 1, this epoch is very early and depends on the 
nature of the cold dark matter. Decreasing n makes the epoch later, and to an increasingly 
good approximation allows it to be calculated without knowing the nature of the cold dark 
matter. After this epoch, structure forms according to the bottom-up picture. 

Let us quantify the extent to which the bottom-up behaviour is satisfied, by estimating 
the the maximum mass M^ax of structures forming at the epoch z^i(M); in other words, by 
estimating the width of the band M to Mmax of the masses that are formed at that epoch. 
It is the mass such that a"o(Minax) is significantly less than 1, so that the density contrast 
filtered on the scale Mmax is still evolving linearly almost everywhere in the universe. To 
be definite, let us require that it is still evolving linearly in 90% of the volume of the 
universe. Prom the Gaussian distribution, this corresponds to ao{M^ax) = 1/1-63 = .61. 
For M = W'^Mq this criterion gives Mmax = 10^ M© (10^° M©) for n = 1 (.7). In the regime 
M ~ 10^ Mq, Zni(M) is a slowly varying function of M, so the conclusion is that a broad 
band of masses downwards of IO^^Mq or so collapses at the the epoch z = Zni{W^MQ). For 
M = IO^^Mq, the criterion gives Mmax = 2 x lO^^M©, more or less independent of n. In 
the regime M ~ z^i^M) is a fairly rapidly varying function, so the conclusion is 

that a fairly narrow band of masses M = 1O^^M0 to about lO^^M© collapse around the 
epoch z = Zni(lO^^MQ). 

We conclude that the bottom-up picture is more or less the same for any n in the range 
.7 < n < 1. A broad range of masses below about W^^Mq collapses at about the same 
time, but subsequent collapse takes place in an increasingly narrow mass range. 

Finally we estimate the actual value of the redshift of nonlinearity ZnilM). It can be 
read from Figure 3b, for n = 1 and n = .7 and one sees that it becomes more recent as M 
or n arc reduced. Let tis consider three representative values of M. First take M = IO^Mq. 
Then Zjil(M) is quite early, 1 + Zni(M) = IScrg (Ocrg) for n = I (.7). A broad range of mass 
scales downwards of about W^Mq {W^^Mq) collapse at around this epoch. Second, take 
M = 1(PMq. Then 1 + Zn\{M) = A.5as {S.Jas) forn = 1 (.7). This is the epoch when 
a significant fraction of the mass of the universe collapses into objects with mass of order 
10^2 Mq, the mass of large galaxies. However, the favoured explanation of a bias factor 
^8 > 1 is that luminous galaxies originate from exceptionally high peaks of the evolved 
density contrast, and therefore form well before the epoch z„i(M). We shall discuss later to 
what extent the qTiasi-linear approach is capable of estimating this epoch. Finally consider 
M = IO^^Mq. Then, 1 -|- Zni{M) = .82(78, almost independent of n. This means that a 



75 



significant fraction of the mass of the universe is only now beginning to collapse into large 
galaxy clusters. The prediction is therefore that the observed large clusters originated from 
exceptionally high peaks of the density contrast, which again implies a bias factor for these 
objects (Kaiser 1984). Again, we will discuss later whether the quasi- linear approach can 
quantify this effect. 

A different criterion is sometimes proposed for the epoch of structure formation, namely 
that it is the epoch when the linearly evolved spectrum Vs{k) is equal to 1 on the corre- 
sponding scale k^^ = Rf- This is practically equivalent to the criterion a{M) = 1 if Vsik) 
is significantly increasing, say like /c™" with m ~ 1, because then the filtered spectrum is 
7'5(A;) exp(— /c^i?^) ~ Vs{R'^^)5{\og{kRf)), where 5 denotes the Dirac delta function. This 
condition fails, in particular, when k corresponds to the maximum of Vs which occurs when 
n < 1, so it is not true that in that case the first structure forms at the epoch when \ + z 
is equal to the maximum value of Vs{k), as was incorrectly stated in Liddle, Lyth and 
Sutherland (1992). 

8.3 The statistics of the collapsed regions 

Now we study the stochastic properties of the collapsed regions, defined as regions where 
the linearly evolved density contrast (5(M, x) exceeds a threshold 6c- Since they represent 
exceptionally large fluctuations of a Gaussian random field, there are powerful mathematical 
results concerning their stochastic properties. We shall use some of them in what follows, 
drawing extensively on the work of BBKS. 

The correct choice of the threshold 6c, insofar as it is well defined, is a matter of debate. 
Many authors take the value 1.69 inspired by the spherical collapse model described in Sec- 
tion 8.1. On the other hand, comparison of the quasi-linear estimate of n(> M) (described 
below) with estimates from numerical simulations suggests a smaller value, Carlbcrg and 
Couchman (1989) advocating 6c = 1-44 and Efstathiou and Recs (1988) advocating 6c = 1.33 
(but see also Brainerd and Villumsen (1992) and Katz, Quinn and Gelb (1992)). 

In the collapsed regions, S{M, x) is more than u standard deviations above zero, where 

iy{z,M) = 6c/a{z,M) (8.4) 
= 6c{l + z)/ao{M) (8.5) 

We are working in the linear regime, corresponding io v > 6c > 1. 

The collapsed regions occupy a volume fraction V given by the Gaussian distribution. 



dV 1 



e-'^ (8.7) 



di^ ^/2^^ 
leading to 

V{iy) = crfc(zy/V2)/2 (8.8) 
= (27r)-VV-^e-"'/2(;L - + 0(j.-4)) (8.9) 

For 1/ = 1,2,3,4 the volume fraction is V = .16, .023, .0013, .000031. In practice one is not 
interested in values > 4, because the collapsed regions are then too rare to be physically 



76 



significant. The corresponding mass fraction is about (1 + (5c) = 2 to 3 times bigger than 
the volume fraction. 

To say more one needs to know the shape of the spectrum. We shall list the relevant 
results given by BBKS. They involve only two moments of the spectrum, defined by 



{k\M)) = a-'^{M) k'^exp{-k'^R})Vsik)-r (8.10) 

JO k 

{k\M)) = a-\M) k'^exp{-k'^R})Ps{k)^ (8.11) 

Jo k 



The quantity (A;^) is the mean of the operator, ie., of the quantity 6~^V^6. Similarly, 
(A;^) is the mean of V^. 

A relevant length scale is defined by = 3{k'^)/{k^). For any spectral index n > — 1, 
it is easy to show that in the limit of small filtering scale Rf, 

* (8.12) 



Rf \l + n 

For the case of CDM with .7 < n < 1, the ratio is in the range 1 to 3 for the entire range 

of cosmologically interesting masses. 

Another relevant length scale is (A;^)~^/^. On large filtering scales, such that Vs{k) is 
increasing fairly strongly at k^^ ~ Rf, the ratio is close to 1. As the scale is 

reduced it increases, but is < 10 for M > 10^ Mq. 

Finally, it is convenient to define the dimensionless parameter 

7(M) = {k')/{k^)'/' (8.13) 

It falls from about .7 to about .3 as M decreases from 1O^^M0 to 1O^M0, for .7 < n < 1. 

For sufficiently large u, each collapsed region is a sphere surrounding a single peak of 
S. However, the departure from sphericity is considerable in the cosmologically interesting 
regime. BBKS show that a quantity x"^, which is roughly the fractional departure from 
sphericity, is well approximated by 

x = ju + 9(^,^1/) (8.14) 

where 

) ^ 3(l-7^) + (1.216-.V)exp|-,/2(W2f] 

[3(1 - 7^) + .45 + (7^/2)2]'/^ + 7I//2 

The asphericity is plotted in Figure 10 for M = IO^Mq and IQI^Mq, for both n = 1 and 
n = .7, and is seen to be > .3 even sX v = A and 10^'^ Mq. We emphasise that this is the 
asphericity seen in the linearly evolved, filtered density contrast. The asphericity in the 
true, unfiltered density contrast will be bigger, and will increase during collapse. 

Three useful number densities are given by BBKS. First, the density riy' of the Euler 
number of the surfaces bounding the collapsed regions is 

ln,(., [k^)) = ^^^llf\ '^' - ^y-^'l^ (8.16) 
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Second, the number density of upcrossing points on these surfaces is 



n^^{v, = ^^^^^ X 

572(1 - 572/9)1/2 



e-^ /2 (817) 



An upcrossing point on a surface of constant 5 is defined as one where V5 points along 
some arbitrarily chosen reference direction. 

The third number density is n^eak, the number density of peaks which are more than 
v standard deviations high. BBKS giv6 expressions for 'Hpeak? 

but they point out also that 

in the cosmologically interesting regime it is quite well approximated by riup. We shall use 
this approximation in what follows. It suggests that if a collapsed region contains several 
peaks, they are not buried deep inside it; rather, the boundary of the region is presumably 
corrugated, wrapping itself partially around each peak. 

In the limit 2/7^1, n^eak = 'T^up = \'^x- This is in accordance with the fact that 
in this regime, each surface is a sphere surrounding a single peak. It contributes +1 to 
the number of upcrossing points, and +2 to the Euler number. As v decreases the surface 
becomes deformed, but at first its contributions to nup and are not affected. Eventually 
though, it may become so corrugated that it has more than one upcrossing point, and may 
become a torus so that its Euler number is less than 2 (of course it then has more than one 
upcrossing point whatever its shape). As a result, its contribution to becomes less 

than its contribution to n^p. The ratio 2npeak/'T'x (approximated as 2nup/'n,j(.) is plotted in 
Figure 11. |^ 

One would like to know the number density ncoU of the collapsed regions, but in general 
there does not seem to be an analytic formula for it. It satisfies the inequality < 
J^coU < '"■peak) and so is equal to in the limit 1/7 S> 1. According to Figure 11, it is 

within a factor 2 of ra^ for 1/ > 3. The average number N of peaks per collapsed region 
satisfies 1 < < 2npeak/n^, so according to Figure 11 it is no bigger than 2 for > 3. 



8.4 The mass of the collapsed regions 

As the filtered density contrast does not contain structure on scales much less than the 
filtering scale Rf, we expect the average radius of a peak in it to be > i?/. It can be 
estimated from the number density of peaks with arbitrary height ?T-peak(~oo), which is 
equal to .016i?^^ where the length was defined at the beginning of the last subsection. 
This is also the number density of troughs, so a rough estimate of the average radius of 
a peak or trough at half height is (2npeak(— C)o))~^/^/4 = .8i?*. As we discussed earlier, 
R<^/Rf ~ 1 to 3 for all cases of interest, so the average peak size is roughly of order 
Rf. Since there are few peaks with radius less than Rf, this suggests that the probability 
distribution of peak sizes is fairly narrow, most peaks having a radius around the average. 

We would like to compare the average peak radius Rf with the average radius of a 
collapsed region. Equivalently, we would like to compare the filtering volume Vf with the 

^®We are not interested in the regime u < 1, but for the record is negative there, indicating that 
the surfaces (5(M, x) = vg definitely do not have spherical topology; in fact they percolate leading to a 
sponge-like topology (Melott 1990). 
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average volume of a collapsed region. The latter is equal to the volume fraction occupied 
by the collapsed regions, Eq. (|8.9| ), divided by their number density tIcoh- In general we 
do not know ncoii, but we do know the peak number density ripeak, so we can calculated 
the average volume Vpeak per peak of a collapsed region. It is given by Eq. (|8.9|), and as a 
fraction of the filtering volume Vy becomes 

W/^/ = ^erfc(i//V^)/(npeak't"/) (8.18) 

The corresponding mass fraction is roughly Mp^^i^/M = (1 + (^c)^peak/^/- The volume 
fraction is plotted in Figure 12, and one sees that except for large it is bigger than 1. 

At large v, the ratio may be calculated using the approximation nup = n^f^^ together 
with Eq. (1^, 

W/^/ ^ (27r)3/2((A:2)/3)-3/2^,-3/y^ (g.ig) 

The fact that this ratio is small does not mean that the filtered density contrast is seeing 
structure on a scale much less than the filtering scale.0 Rather, one is probably dealing 
with density enhancements with size is of order Rf, but with mean overdensity only barely 
exceeding the threshold. When one moves just a little distance away from the peak, the 
mean overdensity within a sphere of radius Rf, represented by the filtered density contrast, 
falls below the threshold. 



8.5 The number density n(> M) 

The main application of these results is to estimate the number density n(> M) of gravi- 
tationally bound systems with mass bigger than M, at a given epoch before z^i^M). The 
systems are supposed to be identifiable by looking at the linearly evolved density contrast 
S{M, x). Each collapsed region, defined as one in which S(M, x) > 6c, is supposed to contain 
one or more systems with mass bigger than M. 

If each collapsed region is identified with a single system, then n(> M) = ricoii- In 
general this recipe is useless for lack of an expression for ricon- A different prescription, 
which does lead to a calculable expression, is to identify each peak within a collapsed region 
with a different collapsed object, 

n(> M) = npeak - ?^up (8.21) 

This estimate (usually without the simplifying second equality) is widely used in the liter- 
ature. It is certainly the same as the estimate n(> M) = ricoii for large v, where we know 
that there is just one peak per collapsed region. To what extent the prescriptions are the 
same for lower v is not known, because the number of peaks per collapsed region is not 
known. 

If at some epoch the linearly evolved density contrast does have many peaks within a 
collapsed region, an interesting situation arises. At a somewhat earlier epoch, 6{M,x) was 

^"We are indebted to Peter Thomas for a clarifying conversation about tliis question. 
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smaller, and a separate contour 6{M,x) = 6c was wrapped around each peak. In other 
words, each peak of the linearly evolved density contrast, filtered on scale M, was inside 
a single collapsed region, and presumably represented a separate gravitationally bound 
system. At the later epoch when the collapsed region encompasses many peaks of the 
linearly evolved density contrast, we have a bigger gravitationally bound system. If the 
original systems survive, the identification of each peak with a separate system is correct, 
but it misses the larger system which contains the original systems. Of course, missing this 
one system does not affect the total count much, so if this case is typical of collapsed regions 
containing many peaks the estimate n(> M) = npeak is better than the estimate n(> M) = 
ricow- If, on the other hand, the original systems have merged, that identification is wrong, 
and the whole of the collapsed region should be identified with just one gravitationally 
bound system. If this case is typical, the estimate n(> M) = ricou would be better, if only 
we had a formula for it. Which case is the more likely? A clue is provided by the observation 
made earlier, that if there are several peaks in a collapsed region they typically seem to lie 
near the surface of the region, a part of the surface wrapping itself around each peak. This 
picture would suggest that the estimate n(> M) = npeak is the more reasonable, the peaks 
of the linearly evolved density contrast in a typical collapsed region representing structures 
which have not existed long enough to merge. 

A different prescription was used by Press and Schechter (1974), to derive a widely used 
alternative formula. They worked with the differential number density, 

(in n 

^ = -^n(> M) (8.22) 
dM dM ^ ^ ^ ^ 

At a given epoch, if the filtering mass M is increased by an amount dM then is increased 
by an amount diy, and the volume fraction occupied by the collapsed regions is reduced by 
an amount dV given by Eq. (8^). Press and Schechter suppose that the eliminated volume 
consists of objects with mass between M and M + dM, corresponding to the idealisation 
that filtering the density contrast on any mass scale M cuts out precisely those objects with 
mass less than M while leaving unaffected objects with mass bigger than M. Ignoring the 
overdensity ~ (1 + 6c) of the collapsed regions this implies that the number density dn of 
such objects is given by 



M 



dn 
dM 



M 



d{R}) 



dM 



d{a^{M)) di/ dV dn 



2Rj 



d{R}) 
-a\M){k'^ 



d((T2(M)) du dV 
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27r 
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Att'^R) 



ve 



(8.23) 



(8.24) 



(8.25) 



Press and Schechter multiplied this formula by a factor 2, so that when integrated over 
all masses it would give the total mass density of the universe, rather than just the half 
corresponding to the regions of space where the linearly evolved density contrast is positive. 
They thus arrived at the estimate 



n(> M) ~ n 



ps 



2\l 



M 



Gtt^R', 



/ 



dM' 



(8.26) 
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In this equation, R'^ = Rf{M'), and similarly for (A;^)' and u'. The factor 2 inserted by 
Press and Schechter is not justified by their argument, because the linearly evolved density 
contrast has nothing to do with reality in the non-linear regime a{M) > 1. On the other 
hand, the neglected overdensity gives a factor ~ {l + 5c) =2 to 3. Thus the factor 2 goes in 
the right direction, and the Press-Schechter formula is reasonably well founded theoretically. 
A somewhat different justification for the formula has been given by Bond et al (1991b). 
From Eq. (0), 

The right hand side is independent of ug, Sc and z, so it follows that Ups, like ripeak, depends 
on these quantities through z^, which involves the combination {5c/crs)i^ + z). In Figure 13 
is plotted the ratio of the two alternative estimates n(> M) = rips and n(> M) = ripeak, 
for M = W^^Mq and for M = IO^^Mq. One sees that the estimates agree to better than a 
factor 2 for <2. Presumably, this indicates that in this regime the assumptions underlying 
the two estimates are compatible, in that increasing M by a small amount cuts out portions 
of the collapsed regions which have mass of order M and are centred on peaks with height 
of order v{M). 

For large v, the Press-Schechter estimate falls below ripeak- This can be understood 
analytically, from the expression 

d^pcak dnup _ 9nup du dn^p d(A;^) dn^p d'j 
dM ^ dM ~ du dM d{k^) dM d-f dM ^ ' ' 

The first term dominates for large u, leading to the ratio 

/ o \ 3/2 



dnps/dM ^ 3 
dnpeak/dM [rW) 



i.29) 



Apart from the factor 2, this is just the filter volume divided by the average volume of a 
collapsed region (Eq. ( ^.20|) ). 

As we saw earlier, the smallness of the size of the 'collapsed regions' is an artefact of the 
filtering. The conclusion is that for very rare fiuctuations, npeak is a better estimate that 
nps, the latter being considerably too small (Thomas & Couchman 1992). As we shall see 
however, other sources of error are likely to be more important than the difference between 
"-peak and nps- 

8.6 Bias factors 

Another result given by BBKS is the bias factor at the present epoch, for a given class of 
objects. The factor is present for any class of objects forming when is significantly bigger 
than 1 ie., which originate as exceptionally high peaks of the density contrast. In general 
the bias occurs partly through the excess clustering of these peaks, and partly through 
additional non-linear clustering after the objects have formed. An estimate including both 
effects is 

6o.j = l + ^ (8.30) 
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where 

^ = .-4^ (8.31) 

In these expressions, everything on the right hand side is to be evaluated at the epoch when 
the objects form. 

8.7 The virial velocity in a galaxy or cluster 

The quasi-hnear approach has been pushed further by some authors, to try to predict 
the virial velocity v{M) of structures of mass M. The virial velocity of a gravitationally 
bound system is defined as the rms velocity of its constituents in the centre of mass frame. 
According to the virial theorem it is given by 

2 GM 



V = -fT- (8-32) 



Rg 

where M the mass of the system and Rg is its gravitational radius, defined by the require- 
ment that its potential energy is —GM'^/Rg. The idea is to relate Rg to the comoving size 
Rcom of the object in the early universe, defined by 

M = {A7t/3)RI^po (8.33) 

where po = 3Hq / (SvrG) is the present mass density. To obtain such a relation one can use 
the spherical collapse model of Section ^.1 in the following manner. First, assume that 
the system has virialised, with no energy loss, before the epoch 2;form when in the absence 
of pressure forces the system would have contracted to a point. Next, assume that the 
system collapsed from an initial configuration at rest, with a density profile of the same 
shape as the profile after virialisation and without energy loss (dissipation). It then follows 
from the virial theorem, Eq. ( ^.32| ) , that Rg is equal to one half of the initial gravitational 
radius. Finally, set the initial gravitational radius equal to the initial radius of the edge of 
the object, defined as the sphere containing mass M. Using the results already quoted in 
Section |8|.2 this gives 



Then Eqs. (|]3|) and (|^ give 



Rf = 8^^com(l + ^max)-' (8.34) 

= 32^R-Ul + ^form)-' (8.35) 

= 178i?-?n(l + Zform)~^ (8.36) 

M \ 



126km s-iy Vl0i2M, 



(1 + ^form) (8.37) 



This is the expression quoted by several of the authors mentioned above (Evrard 1989; 
Henry & Arnaud 1991; Evrard &; Henry 1991), and the others presumably used a similar 
expression though they are less explicit. Obviously, the assumptions leading to this relation 
are at best extremely crude approximations. It seems clear, in fact, that any unique relation 
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between the mass, virial velocity and formation epoch of gravitationahy bound systems can 
only be a rough approximation. 

In Figure 14, the arrows on the horizontal axes show Zform{v), calculated from Eq. ( 8.37] ), 



corresponding to estimates of the upper and lower observational limits on observational 
values of v for objects of mass M. These estimates of v are very rough. For galaxies, 
they correspond roughly to those given by Blumenthal et al (1984) without any attempt to 
update that analysis in the light of more recent observations. For galaxy clusters, the upper 
limit corresponds roughly to X-ray observations (Henry & Arnaud 1991); the observational 
lower limit of about 1000 km s~^ would in that case correspond to negative z, a value of 
about 1300 km corresponding to z = 0. 

The most direct way of utilising the relation Zfoy-ai{M,v) is to set ^^form equal to Zni{M). 
This gives the virial velocity v{M) of systems formed when the scale M goes nonlinear. 
Blumenthal et al (1984) identified such systems with luminous galaxies and with galaxy 
clusters, but according to the biased galaxy formation theory that identification is wrong 
in the former case. Alternatively, the relation can be combined with a theoretical estimate 
of the number density dn{M, z) of objects with mass between M and M + dM, which form 
at epochs between z and z + dz, to calculate the number density n(z, > v) of objects with 
virial velocity bigger than v, which exist at redshift z. This is essentially the approach taken 
by the other authors mentioned above. 

8.8 Comparison with observation 



The above estimates of number densities, sometimes combined with the Eq. (8.30) for the 
bias, have been compared with observation by several authors. The statistic n(> M) has 
been considered for galaxy clusters by Kaiser (1984), Dalton et al (1992), Nichol et al (1992), 
Efstathiou, Bond k White (1992), Adams et al (1993) and White, Efstathiou and Frenk 
(1993), and for both clusters and bright galaxies by Bardeen et al (1986), and Bardeen, 
Bond and Efstathiou (1987). Statistics like n(> v) have been considered by Blumenthal 
et al (1984), Bardeen et al (1987), Evrard (1989), Henry and Arnaud (1991). Evrard and 
Henry (1991) and Adams et al (1993). 

We argue in this section that comparisons of this kind, based on the quasi-linear ap- 
proach, are subject to rather large uncertainties which make them of very limited use. 
In doing so we consider only the statistic n(> M), leaving aside the whole new range of 
uncertainties which come into play when one goes to the statistic n(> v). 

Typically, one will wish to compare the theoretical n(> M) with the observed quantity 
in order to determine the normalisation ag. Because of the strong dependence on Ug this can 
indeed be done very accurately if we know accurately the quantities 6c and M. The problem 
is that these quantities are not known accurately, and in fact are not even well defined to 
high accuracy. In Section ^.1 the collapse threshold 6c was estimated to be .69 on the basis 
of a spherical collapse model, but no-one would argue that this result is not accurate, or 
even meaningful, to even one significant figure. As for the filtering mass M, the prescription 
that the filtered density contrast 6{M, x) 'sees' only those structures with mass bigger than 
M is clearly a crude approximation so one should allow considerable uncertainty in this 
quantity. Yet it is clear from Figure 3b that cto(M) depends quite sensitively on M, and 
even on the type of filter. A related point is that the observed quantity n{> M) depends 
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sensitively on the minimum mass M of the objects counted, which again is not accurately 
known. 

The existence of these uncertainties is well known, and has of course been pointed out 
before (for instance by Efstathiou and Rees (1988)), but at the same time it often does not 
seem to be taken seriously. For instance, many of the studies of galaxy clusters that we have 
mentioned consider the astrophysics in considerable detail, while at the same time adopting 
with little reservation a prescription n-peak or nps with a more or less arbitrarily chosen 
threshold like 5c = 1-69 or 5c = 1.33. A similar approach is taken in a recent discussion of 
the scale dependence of the bias parameter (Bower et al 1993), where it is advocated that 
the threshold 5c varies by a few percent according to environment. 

To understand the basic source of the uncertainty, it is enough to consider the high peak 
limit, in which the prediction Eq. ( 8.16] ) for n(> M) depends only on v, which is given by 



Eq. (^ 

v = 5c{l + z)a^\a^/cjQ{M)) (8.38) 

Let us suppose for the moment that the observational quantity n(> M) is known exactly, 
and estimate the uncertainty in A^g due to the uncertainties in 5c and in the theoretical 
value of M, the latter arising because of the identification of M with the filtering mass. 
One has 



A(ln as) = A(ln 5c) - A (ln(ao/a8)) (8.39) 

)(M), 
dM 



A(ln5.)-^^^5^1^A(lnM) (8.40) 



From Eqs. (|1|) and (|8TC| ), 



^^^^4^4«=(^=> (8.41) 

With n = 1, Ff'ik'^) = 1.0 for galaxy clusters (mass M = IO^^Mq) and .76 for galaxies 
(mass M = IO^^Mq). 

We noted in Section ^.3 that values of 5c to be found in the literature span a range 
1-3 5c < 1.7, so let us take 5c = 1-5 it .2. Taking the theoretical and observational 
uncertainties on M to be each a factor of 2 one finds for galaxy clusters 

^^ = ±.13 ±.23 ±.23 (8.42) 
0-8 

Combining these uncertainties in quadrature gives a total uncertainty of 35%. 

The observational uncertainty in n(> M) at fixed M is an additional source of un- 
certainty, and so is the validity of the prescription 7i(> M) — ^peak 

opposed to, say, 

n(> M) = npg (we argued earlier that the Press-Schechter prescription is too low in the 
large v regime, though it is used by most authors). These uncertainties are likely to be 
less important than the ones we discussed, because of the exponential dependence of the 
prediction on v. 

To be more explicit, we now compare theory with observation in some detail, asking 
if the comparison can distinguish between the three parameter choices {n, cig} = {1,1}, 
{1, .62} and {.7, .62} (these values mark the corners of the roughly triangular region of 
Figure 15 allowed by the bulk flow and cmb (COBE) data — this region is introduced 
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properly in Section Using the prescription n(> M) = Upeak, and the formula Eq. (8.30) 
for the biases bg and ^cluster; we compare theory with observation, first using 6c = 1-33 and 
M = 10^^ (rich clusters) and M = 10^^ (bright galaxies), and then looking at the effect 
of changing 6c to 1.69 and changing M by a factor 2. We will find that the effect is so 
big that the comparison cannot make the required distinction. Note that this is without 
asking about the effect of changing to a different prescription for n(> M), such as the 
Press-Schechter one. 

The prediction n(> M) = Upeak is plotted against redshift in Figure 14, for masses 
M = IO^^Mq, IO^^Mq, W^'^Mq and lO^M©, with np^ak given to sufficient accuracy by 



Eq. (8.17). For each case, three curves are given corresponding to the three parameter 



choices. Each curve ends at the epoch z^iiM), when the linear approach ceases to be valid. 
The threshold has been chosen as 6r = 1.33. 



Galaxy clusters 

Let us consider first the case M = 10^^ Mq, which corresponds to very large galaxy clusters. 
Since this mass corresponds to the normalisation scale Rf ~ 10/i~^ Mpc there is little 
dependence on n. 

Since a{W^^ Mq) = .Sag, we assume that the filtered density contrast on this scale is still 
evolving linearly at the present epoch. The quantities of interest are the number density 
n(> M) and the ratio bc/bg of the galaxy cluster bias factor to the galaxy bias factor. We 
will make the identification bg = l/cig. The prediction for n(> M) is 5.0 x 10^^ Mpc~^ if 
o"8 = 1, and 1.0 x 10~^ Mpc~'^ if as = .63. Following for instance BBKS and Bardeen, Bond 
and Efstathiou (1987), we identify galaxy clusters in this mass range with Abell clusters 
of richness class > 1, and hence observed number density 7.5 x 10"'^ Mpc~'^ (Bahcall Sz 
Soniera 1983). Thus, the prediction for n(> M) is about right if cxg = .6, but too big 
if as = 1. Coming to bc/bg, one finds from Eq. ( ^.30D that if as = 1.6, then be = 4.0 
giving bc/bg = 2.5. If, on the other hand, as = 1.0, then be = 1.5 which gives bc/bg = 1.5. 
Recent estimates (Dalton et al 1992; Nichol et al 1992) give icdi") = (ro/r)^-^='=-^ with 
ro = (13 lb 5)/i~^ Mpc. Dividing by the galaxy correlation function $^gg = {5h~^ Mpc/r)^'^ 
gives therefore bc/bg = 2.4 zt .9 as an observational estimate. Again, as = .6 is preferred. 

So far so good, but what about the uncertainty? Unfortunately it is big. First, suppose 
that we take 6c to be equal to 1.69, instead of the 1.33 used in the above estimates. For 
(Tg = 1 (.62) this multiplies n{> M) by a factor .56 (.16), and multiplies be by a factor 
1.4 (1.4). Second, suppose instead that we multiply M by a factor 2, on the ground 
that there may be this amount of uncertainty in the observational value of the masses of 
galaxy clusters (and remembering also that filtering the density contrast on mass scale 
M does not completely eliminate all structure with mass less than M). This multiplies 
n(> M) by a factor .30 (.09) and multiplies be by a factor 1.6 (1.7). Combining these two 
perfectly reasonable changes gives n(> M) = 8 x 10~^ and bc/bg = 3.3 if as = 1, but 
n(> M) = 7 X 10^*^ and be/bg = 6 if = .62. The former choice of as is now strongly 
preferred by both pieces of data! 
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Bright galaxies 



For galaxies, the linear epoch ends before the present. The number density n(> M) at 
the end of the linear epoch cannot be estimated reliably, and it will also evolve with time 
because of merging and other non-linear phenomena. As a result, no reliable prediction is 
possible for the quantity n(> M) at the present epoch. 

However, in the biased theory of galaxy formation, the formation of luminous galaxies is 
supposed to stop during the linear regime, which gives a bias bg given roughly by Eq. (|8.30| ). 
By demanding that bg value reproduces the value bs = l/fTs one obtains the epoch of 
luminous galaxy formation, and hence the observed number density, as a function of as- This 
approach has been implemented for n = 1 by BBKS and by Bardeen, Bond and Efstathiou 
(1987). Here we extend the calculation to n = .7, and comment on the uncertainty. 

The stars in Figure 14b for ug = .62 indicate the epoch of formation of luminous galaxies 
which is needed to reproduce the bias factor bg = l/ag, according to Eq. ( p.3[)| ). This epoch 
is z = 4.7 (2.8) for n = 1 (.7), and the corresponding value of n(> M) is 4.8 x 10~^ Mpc~^ 
(7.9 X lO-^Mpc-3). I gnoring merging etc., this comoving number density should be equal 
to the presently observed number density ^^g(> M) of luminous galaxies for Af = 
An observational estimate of ng{> M) is the Schechter parametrisation, which taking the 
ratio of luminosity L to mass M to be independent of M gives (Ellis et al 1988) 



(This provides a good estimate for the case M = 10^2^0 that we are discussing here. For 
lighter galaxies the assumption of constant M/L is not very good (Ashman, Salucci & Persig 
1993).) Integrating Mdug/dM over all masses gives = 1.6r2gai x W^^h~^MQ, where Ogai 
is the contribution to 0, of luminous galaxies (including the dark halos) and we take Ogai = 
.1. For M = this gives an observational estimate ng{> M) = 1.75 x 10 ^ Mpc , 

which is marked by an arrow on the y-axis of Figure 14a. 

According to this calculation, the predicted number density of luminous galaxies is 
somewhat too small for erg = .62, and in fact one needs a value as — .8 to reproduce 
it. But now set M = 5 x 10^^ in the theoretical calculation, on the ground that the 
observational value of M could be a factor 2 too high. With ug = .62 this gives a somewhat 
earlier epoch of formation z = 5.9 (3.6) for n = 1 (.7). The corresponding number densities 
are Ua (> M) = 5.8 X 10""^ (1.2 X 10^3), which look much more healthy, and the actual value 
of as needed for consistency is now around .7. Even without looking at any other sources 
of uncertainty (such as galaxy merging, which BBKS emphasise), it seems clear that it is 
very difficult to pin down as within the range .6 < erg < 1 from these considerations. 

Small galaxies 

Finally, consider the results for M = 10^^ Mq and W^Mq, as shown in Figures 14c and 
d. The observed number densities ng(> M) according to the Schechter parametrisation 
are again shown by arrows, (though they become increasingly uncertain as the mass is 




(8.43) 



In this formula 



= 1.56 X 10 



2/1^ Mpc 



(8.44) 
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reduced). Within the biased galaxy formation theory one expects luminous galaxies with 
these masses to form before the epoch z^aiM), but the observational bias factor for them 
is not known, and there is no reason why it should be equal to the result 6 ~ l/cs which 
is observed for bright galaxies. However, the epoch of formation would need to be a lot 
earlier than Zni{M) to give agreement with the observed %(> M). It is presumably more 
reasonable to suppose that merging has taken place, reducing ng{> M) from its original 
value. Again, there does not seem to be a useful constraint on the parameters n and (jg. 

The quasar density 

So far wc have focussed on observations at small rcdsliift. One can also ask about observation 
at z > 1. In particular, quasars have now been seen out to a redshift of about 5, and are the 
oldest observed objects in the universe. In order to produce the observed luminosity, some 
galaxies must have evolved to contain a sufficient concentration of mass-energy. One then 
needs to estimate the minimTim mass reqiiircd and the number density, in order to discover if 
the CDM cosmogony can explain the observations. Such a comparison was made for n = 1 
by Efstathiou and Rees (1988). Using both the npeak and Press-Schechter prescriptions 
they calculated n(> M) for various masses as a function of redshift. Our results are far 
higher than theirs, but the difference can be explained if they chose h = 2.5 and defined b 
with respect to J3 normalisation instead of Gg normalisation (they are not explicit about 
these choices, but our calculation essentially agrees with theirs if these changes are made in 
it). As in the case of large galaxy clusters at the present epoch, the extreme sensitivity to 
normalisation is caused by the fact that u is substantially bigger than 1. 

The number density of quasars assumed by Efstathiou and Rees was 10-6 /i^ Mpc-^ out 
to a redshift of around 4. Since then, observations have become more stringent, and to be 
consistent with present data one requires this number density out to a redshift of 5 (Martin 
Rees, private communication). This value also assumes that the quasar lifetimes are not 
too short, in which case one needs multiple generations. One also needs to know the mass 
which is required to be evolving nonlinearly to harbour the quasar, and they estimated at 
the time that lO^^M© was the smallest that one could safely get away with. However, recent 
work (Martin Rees, private communication) has suggested that only IQ^^Mq is required, 
which corresponds to a large loosening of the constraint. As one sees from Figures 14b or 
14c, the upshot of all this is that enough quasars may form at high redshift, even if n = .7 
and cTg = .62. 

Another cosmological requirement on high redshifts is the Gunn- Peterson constraint on 
the amount of neutral hydrogen in the inter-galactic medium. In the CDM cosmogony this 
implies that some structure has formed before 2; = 5, in order to re-ionise the hydrogen 
(Schneider, Schmidt &: Gunn 1989), but it is not clear how much or of what kind. Even 
with n = .7, structure with M < 10^'' Mq forms at the epoch 1 -|- z ~ 9as so there seems to 
be time to form structure with this n and with as — .6. 

Quasi-linear approach plus numerical simulations 

Many papers calculate n{z, > M) in both the quasi-linear approach also with numerical 
simulations. In a limited regime of {z, M} space, usually corresponding to values of ly which 
are not too big, rough agreement is typically found for a choice of threshold Sc somewhere in 
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the range .3 to .7. Where the numerical simulations actually cover the regime of parameter 
space which is required for the purpose at hand, the quasi-linear approach should probably 
be regarded as redundant. On the other hand, if the quasi-linear prediction is pushed far 
beyond the regime where it has been checked, as is often the case, the uncertainties that 
we have discussed still apply. 

Excluding very low normalisations 

Although the quasi-linear approach is subject to rather large uncertainties, it is capable 
of ruling out very low normalisations. In particular, it seems to rule out n < .6 with the 
COBE normalisation (excluding gravitational waves) on the ground that there are too few 
high redshift objects to be consistent with, for example, the abundance of quasars (Cen et 
al 1992; Adams et al 1993; Haehnelt 1993). As we shall recap at the beginning of Section 
P such values are also excluded by the bulk flows, but this independent exclusion from 
different considerations is certainly significant. 

8.9 Numerical simulations 

Let us now turn to numerical simulations. The idea here is to follow the evolution of a gas 
containing a finite number of particles, which is supposed to represent a finite comoving 
region of the universe. The size of the region considered is taken to be of order 10 to 100 Mpc, 
depending on the observational quantity of interest as well as on the available computing 
power. Starting in the linear regime, the particles are given positions and velocities which 
provide a realisation of the required Gaussian density contrast. They are then allowed to 
move under the action of gravity, perhaps with non-gravitational interactions taken into 
account in some way. 

For the standard n = \ CDM model, many simulations have been reported (Davis et al 
1985; Gott et al 1986; White et al 1987; Carlberg & Couchman 1989; Frenk et al 1988; Frenk 
et al 1990; Carlberg, Couchman &; Thomas, 1990; Bertschinger & Gelb 1991; Couchman & 
Carlberg 1992; Bahcall & Cen 1992; Cen & Ostriker 1992a; Cen & Ostriker 1992b; Suto, 
Cen & Ostriker 1992; Katz, Quinn & Gelb 1992; Katz, Hernquist & Weinberg 1992; Gelb 
1992; Gelb & Bertschinger 1992; Gelb, Gradwohl & Frieman 1993; Gelb & Bertschinger 
1993; Klypin et al 1993; Katz & White 1993; White, Efstathiou & Frenk 1993. A useful 
survey of those prior to 1992 is given by Davis et al (1992a). Recently there have also been 
numerical simulations of the CDM model which allow a tilted spectrum, n < 1 (Vogeley et 
al 1992; Park et al 1992; Cen et al 1992; Cen & Ostriker 1992c). 

Numerical simulations can estimate most of the quantities listed at the beginning of this 
section. The most powerful comparison seems to be with the dispersion (J\\{r) of the line- 
of-sight relative velocity between a randomly chosen pair of galaxies separated by distance 
r. It can be measured on scales r ~ 1 Mpc, where its values is 300 zt 50kms~^ (Davis 
& Peebles 1983). For the well-studied n = 1 case, almost all simulations agree that a 
low normalisation .3 < erg < .5 is required to fit this value. The main exception is the one 
reported by Couchman and Carlberg (1992), who claim that = 1 fits the data on account 
of a large velocity bias. This claim is not generally accepted. 

Going on to the case n < 1, one might expect that tilting the spectrum with as fixed 
would reduce v\\ on the scale ~ 1 Mpc. This is certainly true in linear theory, which 



88 



^ives from Eqs. (|2.68|), (|3.3|) and (|3.4D the result 



4{r) ^ ^(|v(r)-v(0)p) (8.45) 
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Compared with Eq. ( |7.5| ), which gives the dispersion of the pecuUar velocity itself, this 
expression contains the square bracket and as a result cr||(r) probes the spectrum mainly 
on scales k~^ ~ r. The linear approach is not however valid on the scale 1 Mpc where 
observations exist, and on the basis of numerical simulations Gelb et al (1993) claim that 
(T|j(l Mpc) is sensitive to scales well in excess of 10 Mpc. As a result, the prediction with ag 
fixed at 0.5 is rather insensitive to tilt, in the range .7 < n < 1. 

Of the other comparisons with observation, we mention only the number density n(> M) 
of galaxy clusters. We have already discussed the considerable uncertainty in the quasi- linear 
estimate of this quantity. These arise primarily from the uncertainty in the threshold 6c, 
the uncertainty in the theoretical mass M (reflecting the uncertainty in the prescription of 
identifying the filtering mass with the observed mass), and the uncertainty of the observed 
mass. Numerical simulation can in principle eliminate the first two uncertainties, but the 
third remains. The estimate that we made of it in Eq. ( 8.42| ) should remain roughly valid 



for the numerical simulations, because they reproduce roughly the quasi-linear results for a 
suitable choice of 6c- 

The most recent comparison with observation of the number density of galaxy clusters 
is due to White, Efstathiou and Frenk (1993). They compare their simulations with two 
different observational estimates, M = 1.8 x W^'^h'^MQ and M = 1.4 x IO^^/i-^Mq of 
the value of M for which n(> M) = 4 x 8 x 10~^/i^ Mpc~^. They deduce normalisations 
respectively dg = .62 and .54. They emphasise, however, that these estimates of M are 
probably too high, by a factor of perhaps 2 to 3, which would lower as- Using the estimate 



Eq. (8.42) a reduction by a factor 2 from the lower value gives cjg = .41, and a reduction 
by a factor 3 gives ag = -32. The conclusion is therefore that a comparison of simulations 
with the observed number density of galaxy clusters gives the result .3 < erg < .6. This is 
about the same range as the one coming from simulations of the pairwise velocity dispersion, 
except that the upper bound on cjg is somewhat weaker. 



9 Alternatives to pure cold dark matter 

In this final section, we assess the extent to which observation requires a modification of 
the CDM model. We then look briefly at the two modifications which at the present time 
have been investigated in some detail. They both invoke the assumption that the cold 
dark matter density is significantly less than the critical value. In the cosmological constant 
model, the difference is made up by a cosmological constant, whereas in the mixed dark 
matter (MDM) model it is made up by hot dark matter. 
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9.1 Comparison of the CDM model with observation 

Figure 15 shows the main constraints which have been assembled in the last three sections. 
Assuming cold dark matter with Q^n = 1, it gives some observational constraints in the 
n — as plane, where n defines the shape of the spectrum of the primeval density contrast 
and fjg = l/bg defines its normalisation. The standard CDM model has spectral index 
n = 1, but we have allowed the possibility of 'tilt' away from this value. 

The significance to be attached to the various results is probably best left to the reader, 
as the observational situation remains fluid. The lines plotted omit the mean value; instead 
they indicate the allowed spread about that unplotted mean. The constraints plotted are 

• The limits from the COBE experiment. The COBE prediction depends on the grav- 
itational wave amplitude produced, and we have plotted the constraints for the case 
where it is negligible (dashed line) as in the case of 'natural' inflation (n < 1) and 
two-scale inflation (n > 1), and also where it is given by power-law or extended in- 
flation (solid line). The lower line corresponds to the observational upper limit from 
COBE. Prejudice from the South Pole null result (Gaier et al 1992) and from the 
MIT-balloon experiment suggests there is little room for the true result to be much 
above the mean. This can be seen as especially troublesome in the power-law case if 
one wants to flt the APM data. 

• The limits from the QDOT survey (dot-dashed), obtained in the manner discussed in 
Section |^.2. 

• The limits on n (independent of ug) obtained from requiring a flt to the APM data 
(dotted). We have taken our constraint to be 0.3 < n < 0.6 (see Figure 7), though 
we note that Efstathiou, Bond and White (1992) allow a range for their F which 
is roughly equivalent to letting n be as large as 0.67, which is a very conservative 
constraint. These lines should be taken very seriously as limits outwith which the 
theory appears to fail to fit the APM data. 

• The limits from pairwise velocities (dot-dot-dot-dashed). According to Gelb, Grad- 
wohl and Frieman (1993) these are roughly independent of n for reasonable values, 
and we have plotted as such, though one should treat the extrapolation much below 
n = 0.7 with skepticism. Note the inherent contradiction with the QDOT lines in 
the CDM cosmogony. Again, it would appear that one really has to lie between these 
lines in order to satisfy the present observational data. 

The simplest conclusion that one could draw is that standard CDM cannot fit all these 
data. Perhaps one might hope to stretch the QDOT and APM lines towards n = 0.7, 
(78= 0.5 if the gravitational waves are negligible. However, if one takes the QDOT data 
seriously at say the 2-sigma level this possibility appears to be excluded even for natural 
infiation. In any event, the situation appears hopeless for power-law infiation due to the 
large gravitational wave contribution. 

9.2 The cosmological constant model 

Now let us consider the effect of making the present matter density less than the critical 
value 1, making up the difference with a cosmological constant contribution (Peebles 
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1984b; Kofman & Starobinsky 1985; Turner 1991b; Gorski et al 1992; Lilje 1992; Wright et 
al 1992; Efstathiou et al 1992; Kofman, Gnedin & Bahcall 1992; White, Efstathiou & Frenk 
1993). Fohowing most of these authors we set the spectral index n equal to 1. 

The cosmological constant may be thought of as a contribution with homogeneous, time- 
independent energy density and pressure related by PA+PA = 0. Since pm oc oc (1 + z)^, 
the matter density is equal to pA at the epoch 

(1 + Za) = (17A/f)™)l/3 (9.1) 

Until just before this epoch the effect of the cosmological constant is negligible, being < 10% 
of the matter density before the epoch {1 + z) = 2(1 + za)- The theory developed in Section 
§ is then an adequate approximation. The spectrum of the matter density is specified by 
the transfer function, but its shape is affected because of the fact that the epoch of matter- 
radiation equality is earlier, 1 + Zeq being increased by a factor 1/0^- As long as Silk 
damping is neglected the shape of the transfer function is unchanged when written as a 
function of k/keq, or equivalently (assuming Zeq » ^a) as a function of k/Qm, as exhibited 
in Eq. (|1|). 

Around the epoch za the effect of the cosmological constant starts to become important. 
The growth of the density contrast slows down, and ceases altogether well after za- As a 
result the normalisation of the transfer function begins to fall, but its shape is retained. 
According to the discussion in Section 0.1, this means that the slope of the galaxy correlation 
function is fitted with a value of in the range 0.3 to 0.6. 

To calculate the effect of the cosmological constant on other observable quantities one has 
to follow the evolution of the transfer function to the present, and calculate the associated 
peculiar velocity and the Sachs- Wolfe effect (Kofman & Starobinsky 1985). The authors 
cited earlier find find at least marginal consistency with the data, for — -4, corresponding 
to r^A — -6. The most serious potential problem is with the bulk flows, for which the 
prediction seems to be too low. 

The cosmological constant model with il^ ^ ^A may be viewed as unnatural, in that 
we are living at the epoch when the cosmological constant is starting to be significant. 
It does not seem possible to offer even an anthropic justification for this fact, since the 
cosmological constant is negligible until z ~ 1, by which time galaxy and star formation 
are well under way. The most effective way to constrain the scenario appears to be via 
the abundance of lensed quasars (Fukugita et al 1992), which measures the volume at 
high redshift, expected to be much larger in the cosmological constant model than in the 
standard case. New constraints of this type should shortly be available from the Hubble 
space telescope snapshot survey. 

One could, of course, consider the possibility that 

is less than 1 without making up 
the difference at all, < 1. In this 'open' model (so called because the spatial sections are 
open, as opposed to the compact case O > 1) the growth of the density contrast is reduced 
more than in the cosmological constant case, because one has to go to higher redshift before 
the matter dominated O achieves its early time value ~ 1. As a result the model 
predicts a bigger cmb anisotropy for a given normalisation ug, which requires higher values 
of (Blumenthal, Dekel & Primack 1987). 

This is a good place to mention another possibility, which is to reduce the Hubble 
parameter. The effect of such a reduction is again to alter the scale A^max, so that the 
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transfer function remains the same when written in terms of the variable k/hi^ as exhibited 
in Eq. ( |2.45| ). The galaxy correlation function, which is known as a function of /i~^r is 
therefore fitted with a value h ~ .2. The problem, of course, is that such a value is ruled 
out by direct observation, the lowest reported value being around .4. 

9.3 The mixed dark matter (MDM) model 

Many authors, especially since the advent of the COBE data, have considered the possibility 
that the dark matter has a hot component (Shafi & Stecker 1984; Fang, Li &; Xiang 1984; 
Valdarnini & Bonometto 1985; Achilli, Occhionero & Scaramella 1985; Ikeuchi, Norman & 
Zhan 1988; Schaefer, Shafi & Stecker 1989; Holtzman 1989; Schaefer 1991; Wright et al 
1992; van Dalen &; Schaefer 1992; Schaefer & Shafi 1992; Taylor & Rowan-Robinson 1992; 
Davis, Summers & Schlegel 1992; Klypin et al 1992; Shaefer & Shafi 1993; Holtzman &: 
Primack 1993). We discuss this possibility now, again following the majority of authors by 
setting the spectral index n equal to 1. 

Loosely speaking, the term hot dark matter denotes dark matter whose primeval per- 
turbation free streams away on scales up 100 Mpc or so, as opposed to warm dark matter 
where the maximum scale is appreciably lower and cold dark matter where it is too low to 
be cosmologically significant. In practice the term is more specific, denoting matter which 
consists of a massive stable particle species which falls out of thermal equilibrium while still 
relativistic. The only known candidates are the three neutrino species fg, and Vt- 

We briefly summarise some of the main features of such 'neutrino' hot dark matter, 
following Davis, Summers & Spergel (1992). Denoting its mass by rriy, the present mass 
density of the hot dark matter is 

It becomes non-relativistic at the epoch 

For > .1, this is sufficiently early that the epoch Zeq = -62 x lO'* of matter-radiation 
equality is the same as for pure CDM. In that case the scale entering the horizon at Znr is 
given by 

k-^ = {^Y^^ ' ' = ^Q^l" Mpc (9.4) 

On scales smaller than this the perturbation of the HDM free-streams away upon horizon 
entry, causing the growth of the CDM perturbation to slow down. 

This free-streaming makes pure HDM completely unviable. However, if CDM is also 
present, then on becoming nonrelativistic the HDM can fall into the potential wells already 
created by the CDM so that its density contrast grows to match that of the CDM. At any 
epoch, this can occur on scales below an effective Jeans length 

A:7^ = .11(1 + z)^/2j^^ Mpc (9.5) 

By the present time, the CDM and HDM have a common density contrast on scales A;^^ > 
.IMpc, which are the only ones of interest. 
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Because the HDM continues to have significant random motion even after matter dom- 
ination, the transfer function for MDM continues to evolve right up to the present. At a 
given epoch it has less power on short scales than the pure CDM transfer function, in the 
regime 1 to 100 Mpc. As a result, it can generate the 'excess power' on large scales relative 
to short scales discussed in Section ^.1, if Q,y ~ .3. To demonstrate this we have plotted in 
Figure 16 o"o(M) normalised to the pure CDM value, using the transfer function given by 
Klypin et al for fij^ = 0.3. (For illustration, we also show two tilted models, one with and 
one without gravitational waves.) With the COBE normalisation, as = 1/1.7. [Klypin et 
al quote as 1/1.5, but this appears to be because they have normalised to a high value 
of the rms quadrupole rather than directly to the lO'^ result as we do.] Thus one has less 
need of galaxy bias than in the pure CDM model. 

This same value of Qi, is claimed by the above authors to fit various other observations. 
One would particularly like to know the prediction for the pairwise velocity (T|| (r). Numerical 
simulations with MDM are more complicated than for pure CDM, because one has to take 
into account the random motion of the neutrinos, and so far one set of simulations has 
been done with sufficient resolution to look at this quantity (Klypin et al 1992). A fit is 
claimed, but the same simulations run for pure CDM give a lower pairwise velocity than 
most other authors. Specifically, at r = IMpc and with the normalisation as = 1/1-5, the 
simulations give cr|| ~ 450kmsec~^ whereas those of Gelb, Gradwohl and Frieman (1993) 
give o"|| ~ 600 to 700 Mpc (interpolating their results for as = .5 and .75). According to 
the simulations of Klypin et al, going from pure CDM to MDM with as fixed at 1/1.5 only 
reduces the pairwise velocity by around lOOkmsec"^. [Note though our comments above 
that our COBE normalisation gives a lower as for MDM than does theirs.] 

In summary, while it is clear that MDM gives a much lower pairwise velocity than pure 
CDM, when normalised to COBE, it is not yet clear that the reduction is sufficient to fit 
the observed pairwise velocity. 

In Figure 17 the MDM and CDM predictions for o"o(M) are plotted directly. One sees 
that the bottom-up structure formation will proceed more or less as with pure CDM, as 
discussed in Section 8.2. The reduction of small scale power means, though, that there is 
a lack of high redshift objects, and in fact it is clear that a value ily significantly bigger 
than .3 is fatal (Klypin et al 1992; Haehnelt 1993). Whether fi^ ~ .3 is viable remains to 
be seen. 



Fine tuning in the MDM model 

The fine tuning required in the MDM model is of a different kind from that required in the 
cosmological constant model, in that it involves parameters which one might hope eventually 
to calculate from a fundamental theory. One of these is clearly the neutrino mass, which has 
to be around 7eV to make fijy = .3. In addition, there are whatever parameters determine 
the cold dark matter density ~ -7. 

As a purely illustrative example, suppose that the cold dark matter consists of axions, 
and take as an estimate of its density the commonly quoted result (Kolb & Turner 1990) 
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where /„ is the axion decay constant. Then ^ly/^lc is of order 1, provided that /„ ~ 10^^ GeV 
and that rrii, ~ 10 eV. It is natural to suppose that the massive neutrino is the r neutrino, 
in which case one suggestion for its mass is the see-saw formula (Ellis, Lopez k. Nanopoulos 
1992) 

m, = mljM (9.7) 

where M is some very high mass scale. With mtop ~ 100 GeV as suggested by collider 
results, one obtains ~ lOeF if M ~ 10^^ GeV. With these estimates, the question of 
whether MDM is 'natural' boils down to the question of whether the existence of mass scales 
vn-T ~ 100 GcV and ~ M ~ 10^^ GcV is 'natural'. Whether natural or not, such scales 
have certainly been invoked by particle theorists and experimenters, for reasons that have 
nothing to do with cosmology or dark matter (Ellis, Lopez h Nanopoulos 1992). 

Leaving aside specific formulas for and Oo which at the present time cannot to be 
taken seriously, the point we are making is that one hopes eventually to relate the ratio 
of these quantities to masses and couplings which appear in a Lagrangian describing the 
fundamental theory of nature. Then the problem of 'explaining' this ratio will be the same 
as the problem of 'explaining' these masses and couplings, and thus belongs firmly in the 
realm of particle theory or, some would hope, of superstring theory. 

Finally, we comment on the prospect of verifying the MDM model by direct measurement 
of the neutrino properties. The model requires a neutrino, presumably z^r, and a mass of 
order 10 eV and a lifetime > 10-*^° years. The lifetime is obviously inaccessable to any 
expected astronomical observations. What about the mass? It is far too small to measure 
directly through energy and momentum conservation, the present upper limit from such 
techniques being about 35 MeV. However, if they have mass the neutrinos are expected to 
mix, and unless the mixing angle of with i/^ is extremely small it should be observable in 
the forseeable future. At present the limit on the mixing angle if v^- has the required mass 
is about .03, and this can be improved by a factor of 3 or so in the proposed CHORUS and 
NOMAD experiments at CERN and P803 (Ellis, Lopez & Nanopoulos 1992). Thus one has 
a real prospect of verifying the model, or ruling it out, in the forseeable future. 

9.4 Topological defect models 

Finally, let us briefly summarise the present status of the chief rivals to the inflationary 
model for providing initial density perturbations — models based on topological defects 
where inhomogeneities are created in an initially homogeneous background via the Kibble 
mechanism (Kibble 1976). Possible types of defect are gauge strings (Albrecht & Stebbins 
1992; Perivolaropoulos 1993; Perivolaropoulos and Vachaspati 1993), global strings or tex- 
tures (Pen, Spergel & Turok 1992), global monopoles (Bennett & Rhie 1992; Pen, Spergel 
& Turok 1992) and domain walls forming at a very late phase transition (Xiao-chun h 
Schramm 1992; Jaffe, Stebbins & Frieman 1992). 

In these models, the seed (defect) field evolves not under gravity but instead in accord 
with the scalar field dynamics appropriate to the defect field in question. At lowest order, 
one can assume in fact that these defects evolve in a Friedmann background driven by the 
matter content of the universe. The defects act on the matter to provide seed inhomo- 
geneities, but at lowest order these inhomogeneities do not act back on the defect evolution 
(Veeraraghavan h Stebbins 1990). Once gravity is not the only force acting, the evolution 
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of inhomogeneities is vastly complicated, and because of this topological defect theories 
have considerable ground to catch up relative to the inflation-based models in terms of 
theoretical development. These complications are outside the scope of the present paper. 

At the present time, textures seem to be more or less ruled out, because the COBE data 
require a normalisation erg ~ .25 (Pen, Spergel & Turok 1992), too low to explain other types 
of data, given that the perturbation spectrum induced by textures is very similar in shape 
to that of standard CDM. The other cases have not been ruled out, but neither have they 
been confronted with data to anything like the same extent as the CDM model. The most 
promising case is that of gauge strings with hot dark matter (Albrecht &: Stebbins 1992b). 
As already discussed, most theories with hot dark matter are in deep trouble as the dark 
matter free-streaming removes too much small-scale structure for them to be viable at the 
COBE normalisation. Cosmic strings provide a possible exception, as small-scale structure 
on the strings regenerates small-scale inhomogeneities after (or during) free-streaming. On 
the other hand, this process appears to generate rather too much small-scale power in a 
model with strings and cold dark matter (Albrecht & Stebbins 1992a). 

A possible problem for gauge strings is the fact that the temperature of the universe 
cannot be appreciably higher than 10^^ GeV after inflation, as we discussed after Eq. ( |5.46| ). 
This bound implies that cosmic strings with enough energy to form structure probably 
cannot form by the usual Kibble mechanism of a thermal phase transition. (Pollock 1987). 
To see this, one has to compare two results. First, in order to form structure, gauge 
strings would need to have an energy per unit length of at least 2 x (10^*^ GeV)^ (Albrecht 
&; Stebbins 1992; Perivolaropoulos 1993). Second, the energy per unit length of a string 
formed in a thermal phase transition at temperature T is typically ~ T^. Combined, 
these results require a temperature after inflation T > 2 x 10^^ GeV, or an energy density 
p^/^ > {tt"^ / IbY / ^ X 10^^ GeV where > 100 is the effective number of particle species. 
Thus, structure forming gauge strings probably cannot be created after inflation. 

Forming gauge strings, or anything else, before the observable universe leaves the hori- 
zon during inflation is useless because inflation dilutes their density exponentially. There 
remains the possibility of forming the strings after horizon exit, but before the end of infla- 
tion. This can be done, if the string forming field has a suitable coupling to the inflaton field 
or to gravity (Shafi &: Vilenkin 1984; Vishniac, Olive &: Seckel 1987; Lyth 1987; Yokoyama 
1988; Shafi 1988; Yokoyama 1989; Lyth 1990; Hodges & Primack 1990; Basu, Guth & 
Vilenkin 1991; Nagasawa & Yokoyama 1992; Shafi 1993; see also the two-scale model of 
Section |5[4). 

The situation is different for global strings because they are surrounded by a Goldstone 
field. As a result, their effective energy per unit length grows like \'n.H~^ after formation, 
and can increase by a factor of order 100 by the time of structure formation, allowing 
them to form in a thermal phase transition after inflation. A structure formation scenario 
with global strings has yet to be properly investigated, though some tentative, and rather 
unfavourable, conclusions have been drawn by Pen, Spergel and Turok (1992). 
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10 Discussion and Conclusions 



In this paper we have provided a detailed review of the machinery which allows one to 
make calculations in the Cold Dark Matter model and its variants, including considerable 
material which is applicable in more general surroundings. While it is certainly possible to 
have 'designer' models of inflation which are capable of adapting to almost any conceivable 
observational setback, the typical outcomes of presently available models are fairly simple. 
Typically, inflation provides a spectrum of adiabatic gaussian perturbations which can be 
well described by a power-law spectrum, tilted from the Harrison-Zel'dovich spectrum, 
normally tilted so as to provide extra large scale power. The magnitude of the tilt may 
be modest or pronounced. Long wavelength gravitational waves are also produced, in 
most models with amplitude increasing with the degree of tilt in the density perturbation 
spectrum. 

The details of structure formation are very sensitive to the matter content of the uni- 
verse. It appears that if cold dark matter is the main constituent of the universe, present 
observations require that the initial perturbations be adiabatic — isocurvature perturba- 
tions generate excessively large cmb anisotropies for the same final density perturbation. 
Adiabatic perturbations are exactly what inflation provides. In CDM models, the only 
remaining alternative would appear to be texture seeded models, which have been placed 
in jeopardy by a combination of microwave anisotropy and velocity data, though the death 
blow apparently remains to be struck (Pen, Spergel & Turok 1992). 

To end this review, we summarise the present observational situation with regard to 
realistic inflation-seeded CDM models. Our final constraints on models possessing only 
CDM are those illustrated in Figure 15. It is clear that nothing works on all data; adding 
tilt does not appear to offer any immediate advantages over standard CDM. The most 
striking problem appears to be the reconciliation of the pairwise velocity data with the bulk 
flow data (though the former could also be substituted by the galaxy cluster abundance 
(White, Efstathiou h Frenk 1993)). The excess clustering of galaxies on large scales as 
reflected in the slope of the galaxy correlation function is also clearly problematical if one 
takes the bulk flows seriously, but as the hardest piece of data to interpret {ie the easiest 
to dismiss as the consequence of complex astrophysics (see eg Bower et al 1993) ) one could 
argue that it be ignored. 

With nowhere providing clear agreement, it is not too clear how one should assess the 
constraints on n. If one ignores completely the clustering and pairwise data, but decides 
that the bulk flows should be believed at some level not much more than two sigma (and 
assuming at least that the COBE two sigma upper limit cannot be violated), then for 
natural inflation one gets something like n > 0.70. For power-law and extended inflation, 
the large gravitational wave contribution to COBE pushes the limit up to around 0.85. 
(Liddle & Lyth 1992, 1993). This clearly excludes the possibility that these models can 
fit the clustering data, and indeed we have already remarked (Liddle & Lyth 1992) that 
this limit rules out simple forms of extended inflation, as they give a big bubble constraint 
(Liddle & Wands 1991) n < 0.75. (Changing from a cold dark matter cosmogony to a 
different choice also does not seem to salvage the situation (Liddle & Lyth 1992).) 

With a lot of work, one could aim for agreement around n ~ 0.6 - 0.7 and a% ~ 0.5 
- 0.6 in the natural inflation model. Note though that this requires pushing all data to 
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extremes, and in particular relies on the true 10 variance on the mwb being well above the 
COBE mean, which is strongly disfavoured by other microwave results. Such a model does 

however seem to deal adequately with the epoch of formation of structure (contrary to our 
more pessimistic assessment in Liddle, Lyth and Sutherland (1992)). 

To conclude, we have found strong constraints on the slope of the primeval spectrum 
when generated by various inflationary models. Most extended inflation models appear 
to be ruled out completely. Power-law inflation is perhaps viable, but only for values of 
n > 0.85, too high to allow an explanation of the clustering data. Natural inflation (and 
related) models are the most promising candidates for generating useful power-law spectra 
— provided the true level of fluctuations is close to the top of the COBE range they seem 
marginally able to explain the excess large scale clustering as now seen in many optical 
surveys. In most aspects, such a model does at least as well as an unbaised standard CDM 
model, with the advantage of a more plausible galaxy cluster abundance as well as helping 
with galaxy clustering statistics. 

We completed our survey by examining variants on the CDM model which may be better 
suited to explaining the observational data. The standard technique is to utilise additional 
matter (be it a component of hot dark matter or of a cosmological constant) to remove 
short-scale power from the CDM spectrum. Hot dark matter does this by free-streaming, 
a cosmological constant by delaying matter-radiation equality. Because this power can be 
removed over a much shorter range of scales than with tilt, it enables an explanation of the 
observed deficit of short-scale power relative to intermediate scale power in the spectrum. 

We have not attempted a detailed investigation of these models; however, this can be 
pursued using the techniques we have aimed to illustrate in this report. We note though 
that provided modifications to the spectrum remove short-scale power, then the limits on 
the deviation of the spectral slope from unity will tighten. 

The main moral of this paper is that if one is to take inflation seriously as the origi- 
nator of the adiabatic density perturbation, then one must deal with the range of possible 
predictions that inflationary models may make. With options both of tilt and gravitational 
waves, the uncertainty in the inflationary prediction (which can easily halve the short-scale 
power, for example) appears greater than the uncertainty from any other source. In gen- 
eral, the predictions for tilt and gravitational wave contribution to COBE from inflation 
are independent, corresponding to choices for the e and rj slow-roll parameters, thus feeding 
two new parameters into large scale structure studies. 

Indeed, given two new parameters, it is in some ways remarkable that this does not 
appear to provide enough freedom to fix up the CDM model, as Figure 15 illustrates. 
MDM adds yet another new parameter, roughly speaking an ability to remove short-scale 
power from the spectrum while leaving large scales untouched, and may be necessary should 
all present observations stand up. It appears likely that MDM will however need an initial 
spectrum close to n = 1 with no gravitational waves if it is to succeed. Should MDM be 
vindicated, this appears to raise the possibility that a wide range of inflation models not 
providing these conditions with sufficient accuracy may be excluded. 
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Appendix A 



In the text we found it convenient to regard a choice of the time coordinate t as a sHcing 
of spacetime into hypersurfaces labelled by t. Each such slicing was referred to as a gauge, 
and we dealt explicitly with two choices; the comoving gauge where the hypersurfaces 
are orthogonal to the comoving worldlines, and the synchronous gauge where they are 
orthogonal to geodesies which start out comoving in the limit of very early times. Here 
we derive some results relating different gauges, the main object being to prove a couple of 
results which were used in the text. 



Proper time versus coordinate time 

Consider a slicing of spacetime into hypersurfaces labelled by a parameter t. The gradient 
d^t is some multiple of the unit normal vector n^. 



frif, = df,t 



(A.l) 



It follows that du{fu^) is antisymmetric. 







(A.2) 



First take the hypersurfaces to be orthogonal to the comoving worldlines, so that 
is the fluid 4- velocity = dx^/dr, where r is proper time along a comoving worldline. 
Taking the inner product of both sides of Eq. (A.l) with gives / = dt/dr. Then, in a 
basis such that = (1,0,0,0), Eq. (A^) becomes 



^^f-f 



= 
diP 
P + P 







(A.3) 
(A.4) 



The first equation is the necessary and sufficient condition that comoving hypersurfaces 
exist. The second determines /. To first order in the perturbations. 



/ = constant 



1 + 



5p 
P + P 



(A.5) 



Setting the constant equal to 1 corresponds to choosing t to be the average of r over each 
hypersurface. 



The relation STZ = -HSt 

The immediate purpose of the following discussion is to prove the relation TZ = {H/4>)54>, 
which was used in deriving Eq. ( 5.21| ). It was first derived by Sasaki (1986) using the metric 
perturbation formalism of Bardeen (1980), but it is easy also to derive it using the more 
geometrical viewpoint adopted here. 

The quantity on the left hand side of the relation is related to the curvature of comoving 
hypersurfaces by Eq. ( 2.29| ). The infiaton field perturbation vanishes on such hypersurfaces, 
because its spatial gradient corresponds to momentum fiow and in the text we said that 
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5(f) on the right hand side of the equation is to be defined using any gauge (choice of 
hypersurfaces) which is non-singular in the slow-roll limit. Before justifying this statement 
in general, we show that it holds in the gauge where the hypersurfaces have zero intrinsic 
curvature scalar. This will turn out to be true irrespective of the slow-roll limit, as long as 
(j) dominates the momentum density, and will allow us in the next Appendix to extend the 
calculation to that case. 

The relation we want is in fact a special case of a more general relation 5TZ = —H6t, 
giving the change in the intrinsic curvature of a spacelike hypersurface when it is displaced 
by a proper time interval 6t. In this more general relation H and IZ do not refer to a hyper- 
surface orthogonal to comoving worldlines. Rather, they refer to an arbitrary hypersurface. 
Thus, for each hypersurface we define H and TZ as in Section 2, Eqs. (2/7) and ( ^.29|) , except 



that is taken to denote the unit normal to the hypersurface rather than the 4-velocity 
of the fluid. Using a locally inertial coordinate system such that = (1, v^) with of first 
order, 3H = div\ and from Eq. (p39| ) = -AV'^IZ where = did"- and i?^^) is the 
curvature scalar of the hypersurface. 

The proof of the required relation is very simple, if we take as a starting point the 
following relation (Hawking &: Ellis 1973) 

3 6 

In this relation, G^^ = B^^ — i^g^^R is the component of the Einstein tensor along the 
unit normal to the hypersurface, where R is the curvature scalar of spacetime, and R^^^ 
is the curvature scalar of the hypersurface. The quantity H is essentially the extrinsic 
curvature scalar of the hypersurface, so the relation we are using connects the intrinsic and 
extrinsic curvature scalars of a hypersurface, given the curvature of the spacetime in which 
it is embedded. The relation is valid for any choice of hypersurfaces, to first order in the 
perturbations. 

For a comoving hypersurface, Einstein's field equation says that the left hand side is 
equal to SirGp/S, where p is the energy density in a comoving frame. It is therefore not 
necessary that the unit normal is orthogonal to the hypersurface, it need only be of the 
form n'^ = (l,u*) where f* is of first order. The energy density p and therefore G^^ is the 
same for any such choice, to first order. 

Now consider two different choices of hypersurface, and consider the one from each 
choice which goes through a given point in spacetime. The left hand side is the same for 
both choices so the differences 6H and 5R^'^^ are related by 

2H6H = -6R^^^ (A.7) 
6 

We will now calculate 6H, and hence deduce the required change 5TZ. 



To calculate the change in H, start with Eq. ( A.l ), which is valid for any gauge. In the 
limit where the perturbations vanish all gauges have to become equivalent, so / = 1 + x 
where x is of first order in the perturbations. Now consider two different hypersurface 
choices, 

(1 + x)n^ = d^t (A.8) 
(1 + x')n; = d^t' (A.9) 



99 



In any basis such that = (1,0,0,0) up to first order corrections, the relative velocity 
Vi = n'^ — of the normals to the hypersurfaces is given in terms of the time displacement 
6t = t' — t between them by 

Vi = diSt (A. 10) 

The required change in H is therefore given by 36H = diV^ = V'^6t 
Prom Eq. (^]^) it follows that 

HV'^St = -V^7^ (A.ll) 

As we are working with Fourier series (periodic boundary conditions), the unique solution 
of this equation is the advertised relation. 

Now consider inflation. The field perturbation 5(j) vanishes on comoving hypersurfaces, 
so it follows that their curvature perturbation is given as advertised by 

n = {H/(j))6()) (A.12) 

where Sep is the field perturbation defined on those hypersurfaces where the curvature per- 
turbation vanishes. 

In the slow-roll limit the comoving hypersurfaces become infinitely distorted, as wit- 
nessed by the divergence of their displacement 6t = 6(j)/4> from the undistorted zero cur- 
vature hypersurface, as well as the divergence of their curvature TZ. It follows that the 
displacement between gauges remaining finite in the slow-roll limit becomes negligible com- 
pared with 5t, so that Eq. ( A.12| ) becomes valid in all of them. 



Appendix B 

We first explain why the effect of the metric perturbation on the field equation Eq. ( |5.22D 
for (5(/> is of the same order as other departures from slow-roll (Lyth &: Stewart 1990a). 
Including both metric and field perturbations, the field equation is □(/> + V' = 0, where 

□ = (-gy^/^d^(-g)y^g-bd, (B.l) 

To first order, the metric perturbation adds to the left hand side of Eq. ( |5.22| ) a term [dO)cj). 
The unperturbed field (p is position independent so that the spatial derivatives in 60 have no 
effect on it. If it were time-independent as well the metric perturbation would have no effect 
at all. Thus the effect of the metric perturbation vanishes in the limit where the slow-roll 
approximation becomes exact, for any coordinate choice which remains non-singular in that 
limit. (This argument also shows explicitly what one usually takes for granted, that there 
is no back-reaction for a massless (non-inflaton) scalar field; there too, the background field 
has no time dependence.) 

The exact field equation for 5(j), including the effect of the metric perturbation and 
making no slow-roll assumption, has been derived by Mukhanov (1985) in a coordinate 
system with fiat spatial hypersurfaces. In terms of u = a6(f) the equation is 

K + (^^ - z'/z)uk = (B.2) 
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where z = a(j)/H and the prime denotes differentiation with respect to conformal time r/, 
drj = dt/a. Mukhanov derived this equation using the metric perturbation approach, but 
it follows straightforwardly in the fluid flow approach that we are using, from Eqs. ( |5.41| ), 
( 2.25| ), ( p.28 ) and ( 2.29| ). (In deriving it one has to note that the last two equations are 



valid even when IZ is not constant, and that on comoving hypersurfaces 5p = 6P, because 
V{(j)) is unperturbed.) 

One easily checks that the equation reduces to Eq. ( ^.22 ) in the slow-roll approximation. 



in accordance with the general argument above. For the general case, all of the formalism 
in Section H goes through, with awk the solution of the above equation satisfying Eq. ( 5.32 ) 



well before horizon entry. It has been used to calculate the exact spectrum for power-law 
inflation (Lyth & Stewart 1992a), and to estimate the error in the slow-roll approximation 
for an arbitrary inflaton potential, in the slow-roll regime (Stewart & Lyth 1993). 
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Table 1: Current models of structure formation 



non-baryonic 


Structure origin 


cold dark matter 
30% hot dark matter, 60% cold dark matter 
60% cosmological constant, 40% cdm 
hot dark matter 


adiabatic density perturbation 

adiabatic density perturbation 
adiabatic density perturbation 
cosmic strings 



Table 2: Models of inflation 



Model 


n 


R 


f/^Vi GeV 




= l + 2r]-6e 


= 12e 


= 6eV4 X 10^6 


Power-law 


1 - (2/p) 


6(1 - n) 




J9 = 10 


.80 


1.2 


3.4 X 10^^ 


'Natural' 


1 - (2/r) 






r = 10 


.80 


7 X 10-^ 


1.2 X 10^^ 




.96 


4 X 10"^ 


8 X 10^^ 




1 - .0083(a + 2) 


6(1 -n)- 0.1 




a = 2 


.966 


.10 


1.9 X 10^^ 


Two-scale 








m = 10^ GeV 


1.0001 


10~20 


4 X 10^1 


m = 10^3 GeV 


1.14 


.007 


9 X 10^5 



Figure Captions 



Figure la 

The transfer function used to obtain results in this paper, following the parametrisation 
of Efstathiou (1990) based on 4 fitting parameters. This is obtained by scaling the values 
obtained by Bond and Efstathiou (1984) for h=0.75. 

Figure lb 

A selection of transfer functions plotted by fractional departure from that of Efstathiou 
(1990), labelled Tg, as used in this paper. We show Bardeen et al (1986), Davis et al 
(1985), Holtzman (1989), Starobinsky and Sahni (1984) and an alternative from Bond and 
Efstathiou (1984) based on simulations with h = 0.5. Each given as appropriate to CDM 
universes with Q = 1, h = 0.5, low baryonic content and three neutrino species. They 
feature 4, 5, 3, 2 and 4 fitting parameters respectively. Not all paramctrisations have 
attempted to fit down to the smallest scales. As we shall see later, spectra normalised to 
galaxies typically cross at k = 0.1 Mpc"^. Consequently, the relative normalisation of the 
spectra is well defined by the value of the transfer function at this point. The different 
parametrisations show sizeable discrepencies. 

Figure 2 

The normalisation b^duik = IMpc^^) as a function of n. The solid line is from the 
normalisation (used in this paper), while the dashed line shows the rival J3 normalisation. 
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Figure 3a 

The present day spectra, as calculated in the linear approximation, for a selection of values 
of n. One sees the additional large-scale power and the deficit on short scales when one 
compares n < 1 to the standard CDM spectrum. Note k is in Mpc~^, without a factor of 
h. 

Figure 3b 

The dispersion b^a{M) as a function of mass, for n = 1 and n = 0.7, and with both 
choices of filter. The top hat filtered spectra are unity at M ~ IO^^Mq, as required by the 
normalisation. The gaussian filtering gives significantly smaller answers than does the top 
hat, as its smearing gives a higher contribution to the larger scales. 

Figure 4 

This figure is reproduced from the paper of Crittenden et al 1993, which cites references 
for the various experiments. The top Figure shows the spectrum of the cmb anisotropy for 
a slightly tilted standard CDM model with spectral index n = 0.85 and = 0.05. As 
discussed in the text, the case n = 1 is recovered for / » 1 by multiplying the curves by a 
factor = Z'"''^. The contribution of an adiabatic density perturbation is the middle 

line, labelled 'S', and the contribution of gravitational waves is the bottom line, labelled 
'T'. The light dashed line is the density contribution with ^Ib = 0.01. For each curve, the 
quantity plotted is /(/ + 1)E^, normalised to 1 for the quadrupole, / = 2. If the gravitational 
wave contribution to the quadrupole is equal to that of the density perturbation, as is 
roughly the case for power-law inflation with n = 0.85, the top curve indicates the total. On 
the other hand, it could well be that the gravitational contribution is negligible. The bottom 
Figure shows the filters 2Fi for various experiments, as defined in Section ^1. SP89 denotes 
the South Pole experiment of Meinhold and Lubin (the UCB/UCSB MAX experiment has 
a similar filter). Ten denotes the Tenerife experiment, MIT denotes the balloon experiment 
of Page et al and SP91 denotes the Gaier et al and Schuster et al experiment. OVR022 is 
a planned experiment to be undertaken using the 5 meter dish at Owens Valley. OVRO is 
the original Readhead experiment. 

Figure 5 

The triangles indicate the prediction for the mean quadrupole as a function of n for power- 
law inflation, and the stars for natural inflation. The dotted lines plot the COBE observation 
for two choices of bias (note the bias implicit in the y-axis). The CDM prediction for the 
mean at bias 1 is very close to the COBE result. The vertical bars on the starred data 
indicate the spread of the pdf for the quadrupole (for clarity the bars have been omitted 
for the triangles — they are the same size); 95% of the pdf is above the bottom of the 
bars, while 95% is below the top. The xi distribution is not symmetric, so the bars are 
skewed to higher values (somewhat concealed by the log plot). A value of n is allowed at 
95% exclusion if the observations cut through its vertical bar. Modelling the observational 
errors (see text) gives an even looser criterion. 

Figure 6 

The predictions for fog as defined in text, along with the COBE limits (at 1-sigma) 

for bias 1 and 2. The upper line represents the power-law inflation predictions, the lower 
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those for natural inflation. As for the quadrupole, the n = 1 prediction at bias 1 is very 
dose to the COBE result. 

Figure 7 

The predicted angular correlation functions for a choice of n are plotted alongside the 

observational data from the APM survey (Maddox et al 1990). With the anticipated resid- 
ual systematics, values of n between about 0.3 and 0.6 provide reasonable fits, while the 
standard CDM curve falls well below the data. 

Figure 8 

The predicted rms velocity flows, when smoothed with a top hat of radius Rf, for different 
choices of n. 

Figure 9 

The predicted rms velocity flows in a conflguration mimicking the POTENT observational 
data. The velocity field is first smoothed with a 12h~^ Mpc gaussian, reducing the short 
scale power, and then smoothed with top hat filters of radius Rf, giving predictions smaller 
than in figure 9. The solid lines indicate the predictions for n = 1 and n = 0.7. The stars 
indicate the POTENT observational data at bias 1 (read from figure 4 of Dekel (1991)), and 
the triangles the same at bias 1.6. The error bars on the data (the last ones just overlap) 
are 1-sigma. Finally, we emphasise that the theoretical curves are averages over all observer 
points, whereas the observations are a single realisation, with correlated errors due to long 
wavelength domination of bulk fiows. 

Figure 10 

The asphericity parameter plotted against v. 
Figure 11 

The ratio In-a^jn^ plotted against v. 
Figure 12 

The ratio T^eak/^/ plotted against u. 
Figure 13 

The ratio Uap/n-ps plotted versus u. The ratio is not very sensitive to n, and only the case 
n = 1 is plotted. 

Figure 14 

The number density n(> M) as a function of redshift. It is calculated theoretically by 
equating it with the number density of peaks of the linearly evolved density contrast, filtered 
on the mass scale M. Figures 16a-d refer respectively to M = 10^^ M©, 10^° Mq 

and 10^ Mq. For each case there arc three curves. They correspond to the three choices 
{n, b} = {.7, 1.6}, {1, 1.6} and {1, 1} which mark the corners of the more or less triangular 
region allowed by the QDOT and COBE data (figure 15). The arrows on the vertical 
axes give the observed galaxy number densities at the present epoch. The arrows on the 
horizontal axes indicate very roughly the 'observed' range of formation epochs, deduced 
from the indicated range of virial velocities. Each curve ends at the epoch when a{M) = 1, 
signalling the end of linear evolution. According to the theory of biased galaxy formation. 
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luminous galaxies form significantly before that epoch. In the case of bright galaxies, the 
epoch can be calculated by demanding that the resulting bias factor is equal to bg, and it 
is indicated by a star in figure 16 b. The theoretical and observational uncertainties in are 
discussed in the text. 

Figure 15 

The limits on the tilted CDM models in the n-a^ plane, where fog = l/cs- COBE limits 
are shown for two different inflation models, and extrapolated to n > 1 in the case with 
no gravitational waves. The COBE and QDOT limits should be considered as 1-sigma 
bands about an unplotted mean, while the APM and pairwise velocities limits have greater 
significance. See the text. Section 0.3, for a detailed discussion. 

Figure 16 

The dispersion of the density field is plotted for both tilted models {n = 0.7 both with and 
without gravitational waves) and for mixed dark matter. The dispersion is normalised to 
the standard CDM value. The ability of mixed dark matter to fit the clustering data is due 
to the sharp fall-off of the relative power on scales 20-50 Mpc. 

Figure 17 

The dispersion of the density contrast, when normalised to COBE, as a function of mass for 
both standard CDM and for mixed dark matter, utilising the transfer function of Klypin et 
al (1992). 
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